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SUMMARY 

This  report  presents  results  obtained  In  the  final  six  month  period,  July 
1985  to  March  2,  1986  and  summarizes  the  principal  accomplishments  of  the 
entire  period  of  performance  starting  In  November  1982. 

Changes  In  doping  density  near  the  passivated  Interface  of  HgQ^yCdQ 
have  been  measured  using  a  newly  developed  algorithm  (reference  3  of  Sec¬ 
tion  1)  which  allows  accurate  determination  right  up  to  the  Interface 
Itself.  Doping  density  Is  significantly  higher  within  250  to  500A  of  the 
Interface  than  In  the  bulk,  at  distances  beyond  lOOOA  from  the  Interface. 

This  Increase  Is  Interpreted  as  Inhomogeneous ly  distributed  surface  damage  to 
the  Te  lattice.  As  such,  this  Is  another  manifestation  of  the  same  effects 
seen  as:  Cd  loss  from  the  processed  surface,  by  photoelectron  spectroscopy; 
and  low  X  value  regions  at  the  PH0T0X"'-S102  Interface  measured  by  low  tempera¬ 
ture  electroreflectance,  and  reported  In  the  fifth  Interim  technical  report. 

A  significantly  more  detailed  measurement  than  any  previously  reported  of 
the  doping  profiles  resulting  from  attempted  channeled  Implants  of  various 
Ions  Into  the  (110)  surface  of  HgQ^^CdQ^^Te  are  presented.  The  Implanted 
wafer  was  etched  at  a  low  angle,  with  many  MIS  devices  then  made  on  the  new 
surface.  Measurements  of  average  doping  density  across  the  depletion  region 
of  these  devices  produced  a  high  resolution  depth  profile.  Despite  the  use  of 
surface  processing  designed  to  minimize  mechanical  damage  to  the  surface  (and 
hence  maximize  channeling),  the  profiles  are  still  dominated  by  damage 
mechanisms  and  show  an  essentially  gausslan  doping  profile  as  expected  for 
random  Implantation.  Actual  channeling  Is  achieved  only  for  a  small  fraction 
of  the  Incident  Ions  In  each  sample. 

We  discuss  the  role  of  the  weak  bonding  of  Hg  In  the  Hgj_j^Cd^Te  lattice 
on  Interface  morphology  when  metal  Is  deposited  on  the  surface.  The  Ag/ 
HgCdTe,  Cu/HgCdTe,  Al/HgCdTe  and  Cu/CdTe  Interfaces  are  studied  In  ultra-high 
vacuum  with  X-ray  and  ultraviolet  photoemlsslon  and  with  low  energy  electron 
diffraction.  The  three  metals  are  found  to  yield  distinctly  different  Inter¬ 
face  morphologies.  For  Ag/HgCdTe,  the  Hg  Instability  results  In  vacancy  dlf- 
fusion  of  Ag  10  -  10  A  Into  the  semiconductor.  For  Cu  overlayers,  the 

Instability  greatly  enhances  Te  and  Cd  Intermixing  Into  the  Cu  for  HgCdTe 
compared  to  CdTe.  The  more  reactive  A1  over layers  deplete  the  weakly  bound  Hg 
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from  the  HgCdTe  surface.  We  compare  the  metal/HgCdTe  Interfaces  with  the 
corresponding  metal/CdTe  and  metal/GaAs  Interfaces.  For  the  three  overlayer  "  C 

metals  on  HgCdTe  the  Instability  of  the  Hg  bond  results  in  Interfaces  with  a 
much  greater  degree  of  intermixing  than  for  the  corresponding  metal/blnary 
semiconductor  interfaces. 

As  an  aid  in  completing  the  analysis  of  the  angle-resolved  data,  calcula¬ 
tion  of  the  band  structure  of  CdTe  using  the  nonlocal  pseudopotential  of 
Chellkowsky  and  Cohen  (reference  1  of  Section  3.2)  was  performed.  The  spec¬ 
trum  of  computed  energies  and  wave  functions  spanned  the  valence  bands  and  the 
high  lying  conduction  bands  up  to  30  eV  above  the  valence  band  maximum.  The 
results  of  this  calculation  shed  Insight  into  the  role  of  the  actual  final 
states  Involved  in  the  normal  emission  from  CdTe  (and  by  analogy,  from 
Hgj_jjCdj^Te)  in  shaping  the  features  observed  in  the  spectra.  While  a  simple 
model  of  approximating  the  final  states  by  free  electron  parabolas,  success¬ 
fully  accounts  for  much  of  the  observed  structure,  certain  peaks  in  the  data 
reflect  Instead  the  mixing  of  the  parabolic  bands  that  takes  place  in  the 
crystal  potential. 

Angle-resolved  photoemlsslon  spectroscopy  (ARPES)  study  of  the  band 
structure  of  along  the  {110}  direction  using  normal  emission  method 

has  been  published  by  Sllberman  et  al  (reference  1  of  Section  3.3).  This  work 
is  an  extension  of  the  band  mapping  along  other  major  symmetry  directions 
({111}  and  {100})  using  off-normal  emission  method.  By  assuming  a  direct 
transition  model  (OTM)  and  the  free  electron  final  state,  the  band  mapping 
along  {100}  direction  shows  very  good  agreement  with  the  band  calculations  by 
Chellkowsky  et  al,  for  CdTe  using  nonlocal  pseudopotential  method,  and  Hass  et 
al,  (References  2  and  3  of  Section  3.3)  for  HgTe  and  Hg,  Cd  Te  using 
empirical  tight  binding  method. 

The  dislocation  energies  and  hardness  of  semiconductors  are  calculated  by 

an  extension  of  Harrison's  method.  It  is  demonstrated  in  agreement  with  * 

experiment  that  dislocation  energies  per  unit  length  are  proportional  to 
-3  -9 

d  -  d  ,  where  d  is  the  bond  length  and  hardness  is  proportional  to 

d”^  -  d”^^.  The  hardness  is  related  to  the  Interaction  energies  among  dislo- 

c 

cations.  It  is  argued  that  dislocation  densities  of  semiconductors  will  be 
reduced  if  they  are  alloyed  with  a  second  constituent  that  has  a  shorter  bond 
length.  Experimental  evidence  supporting  this  strategy  is  noted. 
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The  bond  length  and  energy  changes  of  the  constituents  of  alloys  relative 
to  their  pure  crystal  values  are  calculated  from  an  extension  of  Harrison's 
method.  It  Is  demonstrated  that  the  already  weak  HgTe  bonds  are  destabilized 
by  adjacent  CdTe,  HgS,  or  HgSe,  but  are  stabilized  by  ZnTe.  It  is  also  argued 
that  dislocation  energies  and  the  hardness  of  semiconductors  vary  as  a  high 
Inverse  power  of  the  bond  length  of  the  constituents.  Hence,  the  shorter 
ZnTe  bond  as  an  additive  should  improve  the  structural  properties  of  HgTe  and 
CdTe.  Experiments  that  support  these  predictions  are  noted.  The  electronic 
transport  properties  of  0.1  eV  band  gap  HgZnTe  are  about  the  same  as  those  of 
HgCdTe,  and  the  structural  properties  of  the  Zn  compound  are  superior;  thus, 
we  conclude  that  HgZnTe  is  likely  to  be  the  better  material  for  IR  devices. 

The  bond-length  relaxation  in  pseudobinary  alloys  can  be  predicted  by  a 
simple  radial  force  model.  In  tetrahedral  structure  alloys  the  bond-length 
deviation  of  the  solute  In  the  dilute  solution  Is  a  quarter  of  the  bond-length 
difference  between  the  two  components.  This  result  agrees  with  the  experi¬ 
mental  work  done  on  the  Gaj_jjIn^Aa  system  performed  by  Mlkkelsen  and  Boyce 
(Reference  1,  Section  4.3). 

A  calculation  is  presented  to  explain  the  anomalous  experimental  behavior 
of  the  Si  2p  core-exciton  binding  energy  and  llnewldth  in  Sl^^Gej^j^  alloys. 

The  observed  minimum  In  the  llnewldth  near  x  ■>  0.13  can  be  explained  as  the 
result  of  a  competition  between  Intrinsic  broadening  due  to  screening  and 
extrinsic  alloy  broadening.  For  pure  Si,  the  binding  energy  Is  estimated  to 
be  0.15  ±0.05  eV  and  the  width  Is  shown  to  be  smaller  than  that  observed  at 
X  -  0.15. 

Correlations  between  atoms  of  Cd  and  Zn  in  Cd^_^Zn^Te  and  Hg  and  Zn  in 
Hg2_^Zn^Te  are  calculated  within  an  extended  quaslchemlcal  approximation.  A 
tetrahedron  of  alloy  atoms  Is  used  for  the  microcluster  unit.  The  statistical 

mechanics  of  the  correlations  is  outlined  for  the  present  problem  and  an  ap-  _ _ 

proximate  method  of  solution  is  Indicated.  Microcluster  energies  are  obtained  U 

from  an  approximation  of  a  Keating  model  due  to  Martin,  plus  a  correction  due  q 

to  Harrison  to  account  for  the  coupling  between  the  tetrahedral  bonds.  It  Is  □ 

found  that  the  correlations  are  significant  at  typical  growth  temperatures  - . . 

(300  to  600K).  predicted  to  undergo  splnoldal  decomposition  be- 

tween  300K  and  600K.  .  . 
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PRINCIPAL  ACCOMPLISHMENTS  OF  THE  HgCdTe  SURFACE  AND  DEFECT  STUDY  PROGRAM 


The  purpose  of  the  study  was  to  examine  the  fundamental  physics  and  chem¬ 
istry  of  interface  formation  on  HgCdTe  to  the  eventual  goal  of  Improving 
Interface  technology,  particularly  surface  passivation. 

•  Passivated  interfaces  of  HgCdTe  are  dominated  by  surface  process¬ 
ing  Induced  damage  to  the  HgCdTe  lattice.  The  damage  extends 
over  an  effective  depth  of  between  200  and  500A  from  the  Inter¬ 
face.  This  Includes  both  grown  (e.g.,  anodic  oxide)  and  depos¬ 
ited  (e.g.,  FH0T0X”'-S102)  passivation  technologies. 

-  Storage  time  and  hence  minority  carrier  lifetime  are  shortened 
(measured  by  capacitance  transient  analysis) 

-  Local  doping  at  the  Interface  Is  Increased  by  at  least  an 
order  of  magnitude  (measured  by  capacitance-voltage  (C-V) 
depth  profiling  using  a  new  data  reduction  algorithm) 

-  Band  structure.  Including  energy  bandgap  Is  altered  (measured 
by  admittance  spectroscopy,  electroreflectance  and  C-V) 

-  Lattice  alteration  Is  associated  with  loss  of  Cd  during 
processing;  Interface  composition  Is  distorted  from 
stoichiometry  (measured  by  photoelectron  spectroscopy  and 
electroreflectance) 

e  These  technology  effects  are  a  direct  result  of  the  underl3'lng 
susceptibility  of  the  HgCdTe  to  be  easily  damaged  due  to  the  weak 
HgTe  bond  which  Is  further  destabilized  by  the  presence  of  Cd. 

e  The  underlying  causes  of  the  bandgap  variation  of 

alloys  and  their  fragile  structural  properties  have  been 
Identified  with  the  unique  nature  of  the  Hg  6s-valence  state 
energy. 

e  Understanding  of  this  effect  has  allowed  us 

-  To  explain  the  breakdown  of  the  virtual  crystal  approximation 
observed  In  the  valence  band  by  PES. 

-  To  Identify  the  reason  for  Hg  and  Cd  core  d-state  energy 
shifts  as  being  caused  by  charge  transfer  from  CdTe  to  HgTe 
bonds. 

-  To  predict  that  the  already  weak  HgTe  bond  would  be  destabil¬ 
ized  by  an  adjacent  CdTe  bond,  while  ZnTe  would  cause  a  less 
deleterious  effect. 

e  Our  study  of  the  strain  coefficients  led  us  to  appreciate 

Harrison's  prediction  that  strain  coefficients  vary  roughly  as 
d~^,  where  d  Is  the  bond  length,  means  that  substitution  of  a 
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short  bond  length  constituent,  e.g.  ZnTe  in  a  longer  bond  length 
host  will  strengthen  the  host. 

•  Based  on  these  concepts,  we  have  suggested  that  Hgt.^Zn^Te  was  a 

good  candidate  material  to  replace  long  wave 

length  Infrared  detectors.  Preliminary  experimental  results 
strongly  support  these  contentions. 

•  Also  based  on  this  theory  suggested  that  JnAs  substituted  Into 
GaAs  Is  a  preferred  way  to  minimize  dislocations  In  GaAs.  An 
explanation  of  the  mechanism  through  which  the  longer  bonded  InAs 
functions  to  decrease  dislocations  densities  was  also  proposed. 

•  Derived  an  entirely  new  hardness  of  semiconductors  that  agrees 
with  experimental  results  for  many  systems. 

•  Properly  predicted  the  observations  by  Mikkelson  and  Boyce  of  the 
bond  length  variations  of  GaAs  and  InAs  in  the  Ga^^j^In^^As  alloy, 
and  predicted  equivalent  results  for  all  group  IV,  and  II-V  and 
II-VI  compound  alloys. 

•  With  no  adjustable  parameters  In  the  theory,  mixing  enthalpies, 
which  are  In  general  agreement  with  experiment,  are  calculated 
for  all  semiconductor  alloys. 

•  Devised  a  new  correlated  statistical  theory  of  order-disorder 
transitions  which  demonstrates  that  as  a  rule,  all  alloys  are 
partially  ordered  (not  random),  and  the  degree  of  ordering  de¬ 
pends  on  growth  conditions. 

•  Appreciated  the  reasons  that  calculated  deep  defects  states  are 
sensitive  to  the  quality  of  the  host  band  structure.  In  particu¬ 
lar,  It  was  demonstrated  that  no  tight  binding  few-neighbor 
calculation  currently  used  Is  capable  of  producing  accurate 
results,  and  most  will  not  even  3d.eld  proper  trends. 

While  we  have  made  considerable  progress  toward  understanding,  the  bulk 
and  Interface  properties  of  hgCdTe  our  most  recent  work  has  raised  several 
questions  that  remain  unanswered.  A  collection  of  the  more  Important  ones 
are: 


e  What  Is  the  Influence  of  strains  that  are  present  during  growth, 
either  because  of  effects  at  the  melt-solid  Interface  In  solid 
state  recrystalline  growth  or  because  of  lattice  mismatches 
between  the  substrate  and  the  alloy  In  epitaxial  growth,  on  the 
cluster  populations? 

•  Do  dipolar  energies  due  to  Madelung  terms  and  electron-electron 
Interactions  play  a  significant  roll  In  cluster  populations? 

(Our  preliminary  answer  Is  yes.) 
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•  What  is  the  relation  between  cluster  populations  and  vacancy 
formation  energies  and  diffusion? 

•  What  Influence  do  the  clusters  have  on  electrical  transport 
properties? 

•  What  drives  Interface  and  metal-contact  structural  behavior? 

•  What  effects  can  one  expect  In  superlattices  as  a  consequence  of 
these  cluster  Ideas? 

•  Make  the  relation  between  alloy  composition,  growth  method, 
strain  coefficients,  dislocation  energies,  and  hardness  explicit. 

There  are  many  other  questions  that  could  be  added  to  this  list,  but 
these  are  a  set  of  pressing  ones  that  should  be  pursued. 

Finally,  even  though  our  Investigation  of  the  feasibility  of  avoiding  the 
damage  due  to  Ion  Implantation  by  channeling  Into  the  {110}  face  was  not 
entirely  successful,  there  are  several  Important  technological  conclusions 
that  can  be  reached: 

•  Channeling  in  HgCdTe  does  occur,  though  it  does  not  dominate. 

The  lattice  damage,  mentioned  above,  which  affects  passivated 
Interfaces  also  clogs  the  channels  at  the  surface  so  that  most  of 
the  Incident  Ions  are  scattered  Into  a  random  Implant  range 
profile. 

•  Damage  Induced  donor  level  follows  the  'random'  portion  of  the 
Implant  profile  but  Increased  In  concentration  by  a  factor  of  10 
(except  In  H  and  CL). 

•  Hydrogen  shows  a  damage  level  below  the  peak  atomic  concentration 
inferring  that  the  H  takes  most  of  the  recoil  resulting  In  little 
lattice  distortion. 

•  Chlorine  Is  the  only  Ion  evaluated  so  far  which  shows  a  damage  to 
Ion  concentration  ratio  of  100.  Chlorine  also  shows  a  shoulder 
In  the  damage  profile  which  roughly  follows  the  channeled  portion 
of  the  Ion  profile. 


vlii 


RPT60411 


SANT*  SARSANA  Of  SEARCH  CENTER 


CONTENTS 


Section 


Page 


Summary  and  Principal  Accomplishments  of  the  HgCdTe  Surface  and 
Defect  Study  Program . . . 

Introduction . 

INTERFACE  STUDY  TASK . 

Introduction. . . . 

Sample  Preparation  and  Measurement.. . . . 

Results . 

Conclusion . . . . 

Acknowledgments . . . . 

References . . . . . 

CHANNELED  IMPLANT  TASK . 

2.1  Channeled  Ion  Implantation  Range  Studies  In  HgCdTe . 

Introduction. . . . . . . . . . . 

Experimental  Techniques  and  Accuracies . 

Experiments . . . 

SIMS  Measurement  Technology . . . 

Discussion . . . 

Conclusions . . . . . . 

References . . . 

2.2  Doping  Profiles  of  Channeled  Ion  Implantation  In  HgCdTe.... 

Implant  Profiles . 

Introduction . . . 

Sample  Prop,  and  Experiment . . . 

Results . . . . . 

Conclusion. . . . . . . 

SURFACE  STUDY  AND  OXIDATION  CHEMISTRY . 

Introduction. . . 

Summary . . 

References . . . . . . . 

3.1  ROLE  OF  Hg  BONDING  IN  METAL/Hgj.j^Cd^Te  INTERFACE  FORMATION 
References . . . . 


RPT60A11 


HUGHES 


SANT*  BAASAA*  Rf  SCAACM  CINTER 


CONTENTS  (Cont) 

-  € 

Section  Page 

3.2  CALCULATION  OF  CdTe  BAND  STRUCTURE .  70 

References .  71 


3.3  ANGLE-RESOLVED  PHOTOEMISSION  STUDY  OF  THE  CdTe,  HgTe, 
AND  Hg,  Cd^Te  BAND  DISPERSIONS  ALONG  MAJOR  SYMMETRY 


DIRECTIONS,  AND  THE  ALLOYING  EFFECT  IN  Hgj.j^Cd^^Te .  72 

References . . . . .  74 


4  ALLOY  BONDING .  77 

4.1  Dislocation  Energies  and  Hardness  of  Semiconductors .  78 

4.2  Effects  Influencing  the  Structural  Integrity  of 

Semiconductors  and  Their  Alloys..... . . . .  81 

4.3  Bond  Length  Relaxation  in  Pseudobinary  Alloys . . .  88 

4.4  Binding  Energy  and  Spectral  Width  of  SI  zp  Core  Exltons  In 

Si^Ge^_^  Alloys .  90 

4.5  Correlations  In  Pseudobinary  Alloys.. .  94 

5  LIST  OF  PAPERS,  PRESENTATIONS  AND  THESES  SUPPORTED  IN  PART 

OR  IN  WHOLE  BY  THIS  CONTRACT .  116 


X 


RPT60411 


N'.^  i’-i' 


J.i.' 


HUGHES 


MAMTA  tAKBAAA  HtUAKCM  CINTIA 


ILLUSTRATIONS 


Figure 

1-1 

1-2 

1-3 

1-4 

1- 5 
1-6 

2- 1 


2-2 

2-3 

2-4 

2-5 


Page 

Measured  Doping  Density  versus  .Depth  for  Distances  Up  to 

lOOOA  From  the  Insulator/HgQ  yCdQ  Interface .  10 

Equivalent  Circuit  of  Measured  MIS  Capacitor  Showing 
Contributions  Due  to  Insulator  ),  Depletion  Layer  (C^), 

Interface  Traps  and  Inversion  Layer  .  12 

l/Cj  versus  Surface  Potential  Excursion  for  MIS  Devices  Made 

With  PH0T0X"-S102  on  n-Type  NgQ^yCd^^^Te .  13 

y 

l/Cj^  versus  Surface  Potential  Excursion  for  MIS  Devices  Made 

with  PH0T0X"-S102  on  p^-Type  HgQ^^CdQ^jTe .  14 

1/Cj^  versus  Surface  Potential  Excursion  for  MIS  Devices  Made 

with  Anodic  Oxlcd  Grown  on  n-Type  HgQ^yCdp  ^Te .  15 

1  MHz  Capacitance  Voltage  Curves  for  the  Samples  Used  In 

This  Study . . . . .  17 

Partial  Periodic  Table  of  the  Elements,  Showing  the  First  Four 
Rows  (Through  Kr)  Plus  In  and  Sb  and  the  Three  Substrate 
Components,  Hg,  Cd,  and  Te.  The  potential  dopants  in 
columns  1,  3,  5  and  7,  which  were  also  the  Ions  Implanted  and 
profiled  In  this  work,  are  outlined  In  boxes .  25 

Channeled  Depth  Distributions  for  Selected  Ions  Implanted  In 

the  {110}  direction  of  HgCdTe  at  150  keV . .  25 

Random  Depth  Distributions  for  Selected  Ions  Implanted  In 

{110}  HgCdTe  at  150  keV .  26 

{111}  Random  and  {110}  Channeled  Depth  Distributions  for 

Ions  Implanted  In  HgCdTe  at  200  keV... . .  26 

Channeled  Depth  Distributions  for  Selected  Ions  Implanted  In 

the  {110}  Direction  of  HgCdTe  at  200  keV  (to  5  pm) .  27 


2-6  {111}  Random  and  {110}  Channeled  Depth  Distributions  for  LI 

Implanted  In  HgCdTe  at  200  keV . . .  27 

2-7  {111}  Random  and  {110}  Channeled  Depth  Distributions  for  F 

Implanted  In  HgCdTe  at  200  keV .  28 

2-8  {111}  Random  and  {110}  Channeled  Depth  Distributions  for  Na 

Implanted  In  HgCdTe  at  200  keV . .  28 

2-9  {111}  Random  and  {110}  Channeled  Depth  Distributions  for  A1 

Implanted  In  HgCdTe  at  200  keV . . .  29 

2-10  {111}  Random  and  {110}  Channeled  Depth  Distributions  for  Cl 

Implanted  In  HgCdTe  at  200  keV . . .  29 

2-11  {111}  Random  and  {110}  Channeled  Depth  Distributions  for  K 

Implanted  In  HgCdTe  at  200  keV . . .  30 

2-12  {111}  Random  and  {110}  Channeled  Depth  Distributions  for  Br 

Implanted  In  HgCdTe  at  200  keV .  30 


RPT604 1 1 


xl 


Figure 

2-13 


HUGHES 


tAMTA  AAMAAKA  Af  Sf  ARCH  Cf  NTf  A 


ILLUSTRATIONS  (Cent) 


{111}  Random  and  {110}  Channeled  Depth  Distributions  for  In 
Implanted  in  HgCdTe  at  200  keV . . . 


Channeled  Depth  Distributions  for  200  keV  B  Implanted  Into 
Different  HgCdTe  Substrates . 


Valence  Band  Spectra  at  hv  *  21.2  eV  of  the  Clean 
Surface  as  Well  as  for  the  Surface  Covered  by  5  ml  of  Each 
Overlayer  Metal.  The  Cu/Hgj_jjCd^Te  spectrum  Is  for  the 
composition  x  »  0.25;  the  other  three  spectra  are  for  x  ■  0.23. 
A  spectrum  of  clean  x  »  0.25  Hgj  Cd^Te  Is  not  shown  as  it  is 
very  similar  to  the  clean  x  »  0.23  spectrum . 


Attenuation  of  the  XPS  Peak  Areas  I  of  the 

Levels  With  Increasing  Cu  Coverage.  The  peak  areas  are 

normalized  to  their  zero-coverage  values  1^ . 

Attenuation  of  the  Semiconductor  Core  Level  Peak  Areas  as  in 
(a),  but  for  Cu/CdTe . . . 


Attenuation  of  the  XPS  Peak  Areas  I  of  the  Hgj_^Cd^Te  Core 
Levels  With  Increasing  A1  Coverage.  The  peak  areas  are 
normalized  to  their  zero-coverage  values  I_ . 


Calculated  CdTe  Conduction  Bands  of  Symmetry  and 
Identification  of  States  With  Greater  Than  30X  (  )  or  Between 
%  and  30%  (A)  [220]  Plane  Wave  Character.  Also  shown  are 
portions  of  the  experimental  final  states  for  peaks  No.  1, 

No.  3,  and  No.  7  (0) . 


The  Band  Mappings  of  CdTe  (•)  and  HgTe  (o)  Along  [100] 
Direction  using  Off-normal  Emission  ARPES.  The  theoretical 
calculations  for  CdTe  from  Ref.  2  (solid  line)  and  HgTe  from 
Ref.  3  (dash  line)  are  also  shown  for  comparison.  Only  the 
lower  bands  are  presented...... . . . 


Comparison  of  the  Theoretical  Calculation  of  (E)  for 
HgQ  yCdQ  jTe  and  the  Experimental  Result  of  the  ARPES  Peak 


Shift  From  the  VCA  Value  of  <100>  Lower  Band. ............... . 

Comparison  of  21"  (E)  for  Hgg  yCdg  and  the  Experimental 
Result  of  ARPES  Peak  Broadening.  T^e  broadening  data  was 
obtained  by  comparing  ARPES  data  of  HgTe  and  Hgn^yCdQ^jTe  of 
the  same  photon  energy  hv,  and  the  analyzer  angle  6^.1 . 


RPT60411 


HUGHES 


Figure 


ILLUSTRATIONS  (Cont) 


Density  of  States  Projected  on  Cation  s  Orbitals  for  HgTe 
(solid  curve)  and  CdTe  (dashed  curve)  [cf.  Eq.  (17)] . 


TABLES 


Heats  of  Te  Fomatlon  and  Heats  of  Solution  AHgQj^(C;M) 

of  Cation  C  in  Metal  M  In  the  Solid  Phase,  for  Various  Metals 
on  Hgj.^Cd^Te . 


RPT60411 


^  y  ' 


ipi  , 


tANTA  tAAAAAA  Af  Sf AACM  CCNTCA 

«  »ubs«i>«<v 

INTRODUCTION 


The  bandgap  tunablllty  of  HgCdTe  between  the  absorption  energies  of 
0.0  eV  (HgTe)  to  1.5  eV  (CdTe)  makes  the  HgCdTe  alloy  system  the  premier 
candidate  for  use  In  Infrared  detection.  The  potential  electronic  properties 
of  this  semiconductor  system  had  been  realized  with  very  limited  success  In 
the  time  period  leading  up  to  the  Inception  of  this  program.  Large  single 
element  photo-detector  devices  had  been  manufactured  with  some  success,  but  as 
more  and  more  sophisticated  device  demands  were  made  on  this  system,  for 
example,  the  manufacture  of  large  scale  photovoltaic  focal  plane  arrays. 
Inherent  problems  In  the  passivation  of  HgCdTe  emerged  which  greatly  reduced 
the  success  of  these  devices.  The  extreme  desirability  of  obtaining  working 
Infrared  Imaging  systems  caused  many  research  and  development  efforts  to 
search  for  quick  solutions  to  the  problems  plaguing  HgCdTe  technology.  It  was 
often  felt  that  merely  extending  Si  and  III-V  semiconductor  process  technology 
to  the  processing  of  II-VI  semiconductor  systems  with  a  few  alterations  would 
produce  successful  yields  for  these  II-VI  devices.  These  empirical  approaches 
have  shown  limited  success  at  reducing  some  of  the  key  problems  hindering 
successful  device  fabrication.  It  was  felt  that  an  Investigation  aimed  at 
examining  the  nature  of  the  bonding  In  the  alloy  and  the  electrical  and  me¬ 
chanical  properties  of  this  alloy  system  and  their  ultimate  effect  on  the 
device  physics  was  a  necessary  first  step  in  solving  the  fundamental  problems 
plaguing  HgCdTe  technology. 


Our  most  Important  results  have  been: 

•  establishing  the  breakdown  of  the  Virtual  Crystal  Approximation 
(VCA), 

•  recognizing  the  close  relationship  of  this  with  the  relative 
Instability  of  Hg  in  the  lattice  as  well  as  tunablllty  of  the 
bandgap,  and 

•  establishing  that  Hg|_jjZn^Te  would  have  important  advantages  over 
HgCdTe  because  of  the  Increased  Hg  stability  In  the  lattice 
(based  strongly  on  the  theoretical  work  of  Arden  Sher  In  the  SRI 
portion  of  this  contract). 

Further  work  on  oxidation  and  metallization  showed  that  the  weak  bondin'^  on  Hg 
Is  a  key  consideration  In  HgCdTe  interfaces  and  makes  the  Interfaces  of  this 
material  very  different  from  those  of  other  zlncblend  materials  such  as  Si, 
GaAs,  or  even  CdTe. 
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In  the  time  period  leading  up  to  this  joint  contract  between  SBRC,  Stan¬ 
ford  and  SRI,  there  were  some  key  observations  and  questions  concerning  HgCdTe 
that  were  not  well-understood.  One  of  the  most  critical  parameters  Involving 
the  optimization  of  device  fabrication,  operation,  and  reliability  is  the 
Interface  found  between  the  HgCdTe  and  any  passivating  or  metal  overlayer. 

For  HgCdTe  the  sticking  coefficient  of  oxygen  had  been  observed  to  be  very 
small  by  Bennlnghoven ,  and  Eblna  found  the  sticking  coefficient  of  oxygen  on 
a  clean  cleaved  CdTe  {110}  surface  to  be  less  than  10~^^.  The  native  oxide  of 
HgCdTe  and  the  simpler  CdTe  were  not  well-understood,  and  the  critical  chemis¬ 
try  between  the  HgCdTe  and  the  native  oxide  layer  was  unknown.  Another  un¬ 
usual  observation  by  Bennlnghoven  was  the  Hg  out-dlf fusion  from  the  bulk  even 
at  room  temperature.  The  difficulty  of  maintaining  structural  Integrity  and 
stoichiometry  during  the  rigors  of  device  processing  were  unknown.  These 
observations  provided  a  possible  link  to  the  behavior  observed  In  practical 
devices.  Surface  leakage  being  a  persistent  problem  In  device  structures 
could  be  explained  by  a  build-up  of  out-diffused  Hg  in  the  Interface  region 
between  the  passlvant  and  the  HgCdTe.  The  difficulty  In  binding  oxygen  to 
stoichiometric  HgCdTe  could  affect  any  passivation  process  on  the  HgCdTe 
surface.  Another  Important  observation  by  Bennlnghoven  was  that  HgCdTe  was 
particularly  sensitive  to  electron  beam-induced  desorption.  Therefore  the 
ability  to  characterize  HgCdTe  using  electron  spectroscopies  without  introduc¬ 
ing  artifacts  or  anamolous  stoichiometry  values  was  unknown.  Given  the  high 
thermod)njamlc  activity  of  Hg  In  the  alloy,  resulting  In  a  very  high  equilib¬ 
rium  vapor  pressure,  it  was  not  even  known  If  HgCdTe  would  be  stable  in  the 
ultra-high  vacuum  (UHV)  environment  employed  during  photoemlsslon  spectroscopy 
experiments  (PES).  There  were  certainly  many  more  unusual  observations  of  the 
HgCdTe  alloy  system  In  addition  to  those  mentioned  above,  but  these  exemplify 
the  unusual  behavior  of  an  alloy  that  appeared  to  have  an  electronic  structure 
that  was  strongly  Influenced  by  the  mechanical  defect  structure,  possibly 
related  to  the  peculiar  nature  of  the  bonding  of  HgCdTe  as  compared  to  GaAs 
and  other  III-V  semiconductors. 

Because  of  our  success  In  employing  PES  to  gain  great  Insights  Into  the 
physics  and  chemistry  of  surfaces  and  Interfaces  of  SI  and  III-V  semiconduc¬ 
tors,  we  felt  that  we  could  draw  upon  our  scientific  experience  with  these 
systems  and  provide  for  a  more  fundamental  framework  of  understanding  Into  the 
HgCdTe  alloy  system.  Since  many  of  the  critical  device  parameters  are  dlffl- 
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cult  to  control  in  this  material,  the  complexity  of  these  parameters  demanded 
a  more  scientific  approach  to  be  employed  In  solving  and  understanding  In  a 
more  fundamental  way  the  characteristics  of  this  material  that  directly  lead 
to  reduced  reliability  and  lower  yields.  The  success  of  this  approach  Is 
manifested  by  the  many  contributions  to  the  present  HgCdTe  technology  from 
SBRC,  Stanford,  and  SRI  while  funded  by  this  contract.  It  would  be  Impossible 
to  do  proper  justice  to  all  of  the  work  that  we  have  conducted  under  this 
contract  In  this  report,  therefore  the  following  paragraphs  are  Intended  only 
as  a  summary  of  the  key  work.  The  complete  list  of  publications  made  over  the 
duration  of  this  contract  that  Is  presented  In  Section  4  gives  testimony  to 
the  successes  and  achievements  of  our  HgCdTe  research. 

Discussions  of  specific  accomplishments  over  the  entire  program  period 
and  results  obtained  In  the  period  July  198S  to  March  1986  are  Included  In  the 
following  four  sections.  Each  section  deals  primarily  with  one  of  the  per¬ 
forming  organizations. 
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DAMA(X  ENHANCED  DOPING  AT  THE  HgCdTe  SURFACE 
J.A.  Wilson  and  V.A.  Cotton 


INTRODUCTION 

The  passivation  of  HgCdTe  surfaces  Involves  the  use  of  certain  chemical 
treatments  which  have  recently  been  shown  to  cause  distortions  In  the  surface 
composition  from  that  of  a  stoichiometric  alloy. ^  This  Is  a  form  of  surface 
damage  which  Is  expected  to  have  a  direct  effect  on  the  electronic  structure 
of  any  Interface  which  Is  later  formed.  Some  compositional  changes  have  also 
been  recently  Identified  at  the  PH0T0X"-S102  Interface,  by  electroreflectance 
(ER),  which  may  be  due  to  the  actual  deposition  of  the  passivation  layer 
Itself.  Both  of  these  effects  are  very  local  to  the  surface,  a  few  tens  to 
hundreds  of  angstroms.  Measurement  of  the  doping  profile  through  the  deple¬ 
tion  region  of  an  MIS  device  formed  on  HgCdTe  Is  probably  one  of  the  most  sen¬ 
sitive  probes  of  the  electrical  consequences  of  these  forms  of  near  surface 
damage.  Composition  distortions  of  up  to  a  few  atomic  percent  are  seen  within 

20  to  3QA  of  the  surface  with  Photoelectron  Spectroscopy  (PES)  while  changes 

14  -3 

In  doping  concentration  on  the  order  of  10  cm  can  be  detected.  This  Is 
nine  orders  of  magnitude  lower  than  the  conventional  levels  detected  by  lES 
and  ER  and  will  allow  detection  of  the  deeper  tall  of  this  surface  damage. 

He  have  measured  the  doping  profile  through  the  depletion  region  of  MIS 
devices  made  on  x  ■  0. 3  using  PH0T0X*-S102  or  anodic  oxide  as  the 

Insulator.  These  profiles  are  calculated  from  high  frequency  and  quaslstatlc 
capacitance  voltage  curves  using  an  algorithm  which  avoids  the  divergence  of 
the  computed  profile  within  an  extrinsic  Debye  length  of  the  Interface,  common 
to  previous  measurements.  A  sharply  higher  doping  concentration  Is  seen  at 
the  Interface  In  all  the  HgCdTe  samples  measured.  The  concentration  drops  to 
the  bulk  value  generally  on  a  scale  of  200  to  500A,  though  some  Increased 
doping  can  sometimes  be  seen  at  lOOOA.  Comparison  of  the  behavior  of  n  and 
p-type  samples  Indicate  the  resulting  surface  structure  is  more  complex  than 
expected  for  the  case  of  one  type  of  uniformly  Introduced  defect.  The  results 
suggest  a  heterogeneous  distribution  of  composition  and  defects  across  the 
surface  under  the  gate  with  neighboring  regions  contributing  n  or  p-type 
behavior. 
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SAMPLE  PREPARATION  AND  MEASUREMENT 

The  MIS  structures  were  made  on  n  and  p-type  RgCdTe  with  a  bulk  composi¬ 
tion  (mole  fraction,  x)  of  about  0.30.  Ingots  were  grown  by  the  horizontal 
zone  melt  technique,  cut  Into  wafers,  polished  to  remove  saw  damage  and  then 
Isothermally  annealed  In  Hg  vapor  to  reduce  the  cation  vacancies  left  from  the 
high  temperature  (bOCC)  growth.  The  anneal  lowers  the  native  defect  acceptor 
concentration  below  that  of  the  residual  donor  Impurities  to  produce  n-type, 
or  not  quite  to  that  concentration  to  produce  p-type.  n  or  p-type  Is  selected 
by  setting  the  anneal  temperature.  Only  the  wafers  outermost  100  to  125  gm 
are  affected  In  this  way  during  the  four  week  anneal. 

The  wafers  were  chemomechanlcally  polished  In  a  solution  of  Br2  in  OMF 
(dlmethylformamlde)  and  then  dip  etched  In  Br2  In  ethylene  glycol  to  prepare 
the  surface  for  passivation.  Passivation  was  performed  by  deposition  of  be¬ 
tween  1500  and  4000A  of  PH0T0X"-S102  or  growth  of  600A  of  anodic  oxide.  The 
anodic  oxide  was  grown  using  a  voltage  source  stepped  In  two  volt  Increments 
to  a  final  voltage  of  15V  during  the  growth.  In  this  process  the  current  den¬ 
sity  Is  not  held  constant.  Gates  consisting  of  400A  of  T1  followed  by  4000A 
of  Au  were  applied  by  thermal  evaporation  through  a  shadow  mask.  No  photo¬ 
resist  techniques  were  employed.  In  addition  a  C-V  standard  made  with  SIO2 
grown  on  SI  was  used  as  a  control  to  verify  the  data  reduction  technique. 

Capacitance  measurements  were  made  at  1  MHz  using  a  PAR  410  capacitance 
meter  and  at  quaslstatlc  using  an  electrometer  (ramp  rate  Is  between  10  and 
50  mV/sec).  Typical  current  levels  are  10~^V.  Measurements  were  made  In  a 
dewar  cooled  with  12^2  with  zero  field  of  view. 

The  C-V  data  was  reduced  to  give  surface  potential  versus  applied  gate 

n 

bias,  l/Cjj  versus  surface  potential  (^g)  and  doping  versus  depth,  using  an 
algorithm  based  on  the  analysis  presented  In  reference  3. 

Integration  of  the  quaslstatlc  capacitance  curve  gives  the  dependence  of 
surface  potential  (i^.)  versus  applied  bias,  to  within  an  additive  constant.^ 

9 

The  doping  density,  N,  as  a  function  of  depth,  W,  Is  given  by 


N(W)  - 
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where  g  Is  defined  by 


W 

X 


(g  -  In  g  - 


The  value  of  g  used  in  determining  the  doping  Is  obtained  from  the  slope  of 
2 

l/Cj  versus  by  a  search  routine  fitting  the  expression 


,1^C2  ‘I 


q  -d  d4.^  1-g  „2 

^  72 


where  the  extrinsic  Debye  length  Is  given  by 

2k  Te  ,  /- 

»-  It — s-)‘". 

The  value  of  l/C^  la  determined  from  the  high  frequency  capacitance 
voltage  data  by  the  expression 


1/c/  -  (^-1)^ 

where  Is  the  Insulator  capacitance,  C|||.  Is  the  measured  high  frequency 

capacitance,  k  Is  Boltzmans  constant,  q  Is  the  electron  charge,  e_  Is  the 
dielectric  constant  for  HgCdTe  and  T  Is  the  temperature.  The  negative  sign  Is 
for  n~type  and  the  positive  sign  for  p-type  HgCdTe.  The  additive  constant  In 

resulting  from  the  Berglund  Integration  can  be  determined  from  the  Intei^ 

2 

cept  of  the  tangent  to  the  portion  of  the  1/C^  versus  curve,  for  (^s  In 

2  6 

depletion,  with  the  horizontal  axis  (1/C^  *  0)  .  It  will  become  apparent 

from  the  data,  that  In  many  cases,  the  region  of  the  l/C^  curve  through 
depletion  Is  not  closely  approximated  by  a  single  slope.  This  Indicates  that 
there  Is  a  no  single  appropriate  offset  constant  and  therefore  no  single 
flatband  condition  for  the  bands  near  Interfaces. 


RESULTS 

Doping  profiles  for  samples  made  on  both  n  and  p-type  HgCdTe  using  PHO- 
T0X”-S102  as  the  Insulator,  and  on  n-type  using  anodic  oxide  show  measured 
doping  levels  using  this  technique,  which  agree  with  those  expected  for  the 
bulk  material  at  depths  far  from  the  Interface;  0.1  um  or  deeper  and  agree 
with  Hall  data  with  a  factor  of  2.  Near  the  Interface,  however,  the  measured 
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doping  Is  seen  to  depart  significantly  from  the  bulk  values  for  all'  the  I]gCdTe 
samples  as  shown  in  Table  1-1.  Table  1-1  compares  peak  surface  region  doping 
to  that  In  the  bulk,  at  depths  between  half  and  one  ym.  The  measured  bulk 
doping  Is  uniform  at  these  depths.  The  Increased  apparent  surface  dopings  are 
at  least  an  order  of  magnitude.  Figure  1-1  shows  the  portion  of  the  doping 
profiles  within  lOOOA  of  the  Interface.  Also  shown  Is  a  doping  profile  of  a 
SI  MIS  device,  derived  using  the  same  algorithm,  which  shows  no  significant 
change  In  Its  bulk  doping  near  the  Interface.  The  interface  value  of  the 
doping,  for  the  SI  device.  Is  within  the  range  of  fluctuation  seen  over  the 
entire  depletion  region,  about  2  x  lO^  cm~  .  This  high  degree  of  fluctuation 
Is  due  to  the  hand  digitizing  required  of  this  particular  data  set.  The 
HgCdTe  data  was  recorded  directly  by  computer  and  consequently  shows  much  less 
fluctuation  over  the  entire  depletion  region  for  each  sample. 

Table  1-1. 


Sample  I  MIS 


PHOTOX 

PHOTOX" 

PHOTOX" 

PHOTOX" 

PHOTOX" 

Anodic 

Anodic 


Doping 

Surface 

Bulk 

cm“^ 

.nn-3 

rm 

2x10*^ 

2.2x10^^ 

Ixio^^ 

5x10^^ 

3x10^^ 

6x10^^ 

4x10^^ 

2.5x10^^ 

IxIO^^ 

1.8x10^^ 

1x10^^ 

2x10^^ 

1x10^^ 

2x10^^ 

2x10^^ 

2x10^^ 

Surface  Potential 
Range 
eV 


Figure  1-1  shows  doping  profiles  for:  three  n-type  devices,  numbers 
7161,  2026  and  2069;  two  p-tjrpe  devices,  3001  and  3002,  all  made  using 
PH0T0X"-S102 ;  and  two  n-type  devices  made  with  Anodic  oxide,  005A  and  0058. 
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In  all  cases  the  surface  doping  Is  higher  than  the  bulk  by  over  an  order 
of  magnitude  while  In  the  case  of  the  anodic  oxide  the  Increase  Is  nearly  five 

orders  of  magnitude.  It  can  also  be  seen  In  Figure  1-1,  that  this  region  of 

higher  doping  Is  almost  entirely  confined  to  the  250A  nearest  the  Insulator 
Interface.  We  also  note  that  some  of  the  curves  are  Incomplete  (broken  lines) 
In  that  some  data  Is  not  plotted  between  the  surface  and  the  bulk  In  the  re¬ 
gion  where  the  doping  Is  expected  to  be  changing  the  most  rapidly.  These 

points  are  rejected  during  a  step  In  the  calculation  of  the  doping  density  due 

to  lack  of  convergence  of  the  search  routine.  Physically  these  are  regions 
where  the  variation  of  doping  Is  too  rapid  and  the  algorithm  used  breaks  down 
(the  value  of  1/C^  varies  very  rapidly  with  .  Each  point  Is  calculated 
Independently  of  neighboring  points  so  that  omission  of  any  point  has  no 

effect  on  the  accuracy  of  the  others.  The  calculated  doping  profiles  are 

2 

reduced  from  the  slope  of  the  dependence  of  1/C^^  on  The  calculation  of 

1/Cj  from  the  data  Is  straight  forward  and  does  not  Involve  the  same  search 
routine  so  that  the  points  missing  from  the  doping  profiles  are  contained  In 
these  curves. 

Depletion  layer  capacitance  was  determined  from  the  high  frequency  capac¬ 
itance  voltage  data  based  on  an  MIS  equivalent  circuit  of  the  form  shown  In 
Figure  1-2.^ 

High  frequency  measurements  were  made  at  1  MHz  which  Is  high  enough  that 
no  known  HgCdTe  Interface  state  can  contribute  significantly.^  The  curves  of 
l/Cjj  shown  In  Figures  1-3  through  1-5  are  accurate  for  surface  potentials 
from  accxmulatlon  (C^j  Is  Infinite)  through  depletion  to  the  onset  of  Inver¬ 
sion.  Since  the  doping  Information  displayed  In  Figure  1-1  Is  contained  as 
the  slope  of  these  curves.  It  Is  significant  that  the  slope  is  not  constant 
over  the  surface  potential  range  representing  depletion,  to  the  same  degree  as 
the  SI  device.  This  Indicates  that  the  flatband  condition  is  in  general, 
depth  dependent  In  I^CdTe.  A  single  flatband  bias  is  a  useful  concept  In  the 
case  of  SI  MIS  devices  and  is  a  useful  approximation  for  most  MIS  devices  In¬ 
cluding  HgCdTe,  but  It  Is  not  strictly  correct  at  these  Interfaces. 

2 

The  value  of  1/Cj  versus  for  SI,  shown  for  comparison  as  the  inset  in 
Figure  1-3,  Is  essentially  a  straight  line  (the  small  fluctuations  are  again 
an  artifact  of  the  hand  digitization  used  on  this  data)  while  those  of  the 
HgCdTe  samples  are  of  nonuniform  slope.  From  the  plots  of  l/C^  versus  (fg  it 
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Figure  1-2.  Equivalent  Circuit  of  Measured  MIS  Device  Showing 

Contributions  Due  to  Insulator  Depletion  Layer  (Cj), 

Interface  Traps  and  Inversion  Layer  (C^nv*%nv^* 

is  clear  that,  with  the  exception  of  0054,  there  are  no  clear  Indications  of 
general  type  reversal  and  therefore  of  any  p-n  junctions  In  series  with  the 
depletion  layer  capacitance.  Type  reversal  would  be  Indicated  by  a  region 
where  the  slope  changes  sign.  In  sample  0054,  Figure  1-4,  there  Is  a  change 
of  sign  of  the  slope  between  ■  -0.2  and  -0.25  eV.  The  substrate  of  this 

anodic  oxide  MIS  Is  n-type  so  chat  In  this  region  It  appears  there  Is  local 
type  conversion.  The  equivalent  depth  of  this  feature  In  Figure  1-1  Is  the 
region  beginning  around  250A  and  extending  another  several  hundred  angstroms 
Into  the  bulk.  In  general,  however,  the  increased  surface  doping  Is  of  the 
same  type  as  the  bulk  material  and  there  Is  no  clear  evidence  of  a  junction 
layer  at  the  Interface  so  that  the  damaged  region  Is  probably  In homogeneously 
distributed  across  the  surface  with  neighboring  regions  being  either  n  or 
p-type.  This  allows  Iro-type  paths  to  short  and  heterotype  barriers  In 
general. 

Since  the  Increased  effective  surface  doping  Is  of  the  same  type  as  the 
bulk.  It  Is  unlikely  that  It  Is  caused  directly  by  a  damage  mechanism  such  as 
that  which  always  results  In  donor  ac  '.Ivity  following  Ion  Implant  into 
HgCdTe.  The  increased  doping  most  likely  results  Indirectly  from  damage  to 
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the  Te  lattice  resulting  In  net  alteration  of  the  band  structure  In  such  a  way 
that  the  ferml  energy  Is  relocated  closer  to  the  nearest  band  edge  than  In  the 
undamaged  bulk.  We  note  that  the  effect  Is  slightly  more  pronounced  In  p  than 
n-type  HgCdTe.  However  this  Increased  effect  cannot  be  taken  by  Itself  to 
Indicate  that  the  valence  band  edge  Is  more  strongly  affected  since  the  p-type 
sample  were  about  an  order  of  magnitude  higher  to  start  with.  The  effect  Is 
significantly  more  pronounced  at  the  anodic  oxide  Interface  Indicating  that 
the  applied  PHOTOX”  SIO2  passivation  is  Inherently  a  more  benign  process. 

As  Indicated  In  Table  1-1  the  total  excursion  of  the  surface  potential 
for  these  Interfaces  does  not  correspond  closely  to  .  nor  vary  simply  from 
that  expected  from  the  bulk  bandgap  at  77  for  x  -  0. 3  HgCdTe,  of  0.26  eV.  The 
PHOTOX"  Interfaces  on  n-type  wafers  all  show  a  more  restricted  excursion  while 
the  other  cases,  particularly  the  anodic  interfaces,  show  a  wider  excursion. 
The  surface  potential  excursion  Is  the  value  obtained  directly  from  Berglund 
Integration  of  the  quasi-static  capacitance  voltage  curves.  For  MIS  devices 
which  exhibit  a  measurably  long  storage  time,  greater  than  about  2  ms,  there 
is  the  possibility  that  the  MIS  device  does  not  closely  approximate  equllib- 
rlun  while  the  gate  bias  Is  ramped  from  accumulation  to  Inversion.  In  this 
case  the  low  frequency  curve  does  not  accurately  represent  the  true  quasl- 
statuc  condition  the  derived  surface  potential  excursion  appears  wider  than 
the  bulk  bandgap  and  any  surface  state  densities  derived  from  this  data  would 
appear  artificially  low.  The  parts  used  In  this  study  were  selected  to  avoid 
bulk  storage  times  long  enough  to  cause  this  effect.  Also  the  results  ob¬ 
tained  for  bias  sweeps  from  accumulation  to  Inversion  were  compared  to  those 
for  sweeps  from  Inversion  to  accumulation,  to  avoid  any  nonequilibrium  condi¬ 
tions.  The  nonequivalence  of  the  surface  potential  excursion  and  the  bulk 
energy  gap  Is  taken  as  a  real  effect  of  these  Interfaces  Indicating  a  major 
alteration  of  the  interface  crystal  structure. 

The  part  of  the  C-V  curve  which  gives  rise  to  these  apparent  Increases  in 
surface  doping  are  seen  as  a  very  slight  flattening  of  the  curve  at  the  tran¬ 
sition  between  accumulation  and  depletion.  This  region  Is  Indicated  by  the 
arrows  In  Figure  6.  One  of  the  samples  (2026)  was  chosen  because  the  flatten¬ 
ing  is  especially  pronounced.  Similar  types  of  distortions  of  the  C-V  curves 
can  occur  due  to  the  presence  of  discrete  or  very  high  densities  of  interface 
traps.  The  densities  of  Interface  traps  for  the  devices  of  this  study  are 
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typical  of  l^CdTe  with  nldgap  nlnlna  of  about  1  to  5  x  10^^  eV~^  cm''^  and  are 
free  of  unusual  densities  of  discrete  traps.  Also  the  degree  of  flattening  Is 
similar  In  both  high  frequency  and  quasl-statlc  curves  so  that  any  interface 
state  associated  with  this  feature  would  have  a  response  frequency  well  above 
1  MHz.  To  a  detection  limit  of  about  5  x  10^  eV“^  cm“^  there  are  no  high  fre¬ 
quency  states  present  at  these  Interfaces.^  This  type  of  distortion  of  the 
C-V  curve  cannot  be  due  to  any  fundamental  property  of  the  HgCdTe  lattice, 
such  as  non  parabolic  It  y  of  the  conduction  band,  since  the  effect  varies  trldely 
In  degree  from  device  to  device  for  material  with  the  same  composition.  It  Is 
also  not  a  peculiarity  of  the  algorithm  because  the  Increased  doping  corre¬ 
sponds  In  degree  to  a  real  feature  of  the  capacitance  voltage  data. 


GATE  WAS  (V) 


Figure  1-6.  1  MHz  Capacitance  Voltage  Curves  for  the  Samples  Used  in 

This  Study.  Arrow  indicates  the  bias  region  where  nonldeal 
structure  Interpreted  as  Increased  near  Interface  doping. 


CONCLUSION 

We  have  examined  the  dependence  of  doping  on  depth  at  the  passivated 
Interface  of  HgCdTe  formed  with  anodic  oxidation  or  deposited  PHOTOX"  S102< 
In  all  cases  the  doping  density  within  250  to  500A  from  the  Interface  Is 
Increased  sharply  over  that  of  the  bulk  material.  Because  the  Increased 
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doping  Is  generally  of  the  same  type  as  the  bulk,  this  Increase  Is  taken  to  be 
the  Indirect  result  of  a  fundamental  Inhomogeneous  alteration  of  the  Te  lat¬ 
tice  and  hence  Its  electronic  band  structure  by  the  processes  used  to  form  the 
MIS  device.  As  the  result  of  this  lattice  distortion  the  bands  are  curved 
near  the  Interface,  to  some  extent,  at  any  applied  bias  and  the  energy  gap  Is 
changed  from  the  btilk  value.  The  Increase  In  effective  doping  and  band  dis¬ 
tortion  are,  hence,  consequences  of  the  weakened  HgTe  bond  In  the  presence  of 
Cd  which  Is  the  fundamental  cause  of  the  ease  with  which  this  alloy  Is  dam- 

O 

aged.  In  general,  passivated  Interfaces  on  HgCdTe  are  broader  and  more  rich 
In  structure  than  the  corresponding  situation  on  other  semconductors . 
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2.1  CHANNELED  ION  IMPLANTATIION  RANGE  STUDIES  IN  HgCdTe 
INTRODUCTION 

We  have  measured  the  depth  distributions  for  elements  of  Interest  as 
potential  dopants  Implanted  in  the  {111}  random  and  the  {110}  channeling 
orientations  of  HgCdTe  using  secondary  ion  mass  spectrometry  (SIMS).  The 
HgCdTe  material  used  was  preannealed  solid  state  regrown  (SSR).  The  elements 
studied  were  all  of  odd  Z,  taken  from  columns  1,  3,  5,  and  7  of  the  periodic 
table,  namely  H,  Li,  Na,  K,  Cu,  B,  Al,  Ga,  In,  N,  P,  V,  As,  Sb,  F,  Cl,  Mn,  and 

Br,  plus  Be.  Ihe  implantation  energies  were  150  or  200  keV.  The  implant 

1  ^  —2  15  —2 

fluences  were  mostly  3  *  10  cm  ,  with  extrapolation  to  1  x  10  cm  be¬ 
cause  of  SIMS  measurement  technology.  Alignment  for  the  channeled  Implants 
was  performed  by  using  Rutherford  backscatterlng  of  protons.  We  have  deter¬ 
mined  the  range  parameters,  Rp,  ARp,  maximum  channeling  range  R^. ,  and  the 
shape  factors  or  third  and  fourth  moments,  and  02>  performing  pearson  IV 
fitting  routines  to  the  experimental  depth  distributions.  Sufficient  detec¬ 
tion  sensitivity  using  SIMS  was  achieved  to  determine  many  {110}  channeling 
ranges,  which  vary  substantially  and  reflect  the  Z^-dependent  electronic 
stopping  S_.  Some  channeling  ranges  in  the  large  {110}  channel  of  HgCdTe  are 
from  6  to  50  pm  for  200  keV  ion  energies. 

When  ions  are  channeled  in  a  crystal,  they  lose  energy  by  electron  inter¬ 
actions,  not  by  nuclear  collisions.  Therefore  they  do  not  cause  displacement 
damage  or  vacancy-interstitial  pairs,  until  near  the  end  of  their  range  when 
they  have  lost  all  but  the  last  few  keV  of  energy.  Therefore,  channeled  ion 
implantation  does  less  damage  in  the  deeper  channeled  region  of  the  implanted 
depth  distribution.  Damage  is  still  done  in  the  surface  region  by  ions  that 
do  not  enter  channels  that  contribute  to  the  random  component  of  a  channeled 
depth  distribution. 

The  motivation  for  this  work  was  to  determine  implantation  range  param¬ 
eters  and  depth  distributions  for  selected  elements  implanted  in  both  the 
random  and  the  open  {110}  channeled  orientations  of  HgCdTe  crystals  that  wsuld 
be  used  to  fabricate  devices.  In  particular,  we  desired  to  learn  how  deep 
these  selected  ions  channel  and  what  atom  densities  can  be  achieved  in  the 
deep  channeled  regions,  with  the  idea  that  the  damaged  surface  region  could  be 
etched  away  to  leave  a  less  damaged  region  that  might  be  chemically  doped 
rather  than  defect  doped.  The  random  range  parameters  and  profiles  are 
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Important  for  understanding  standard  Implantation  technology.  Establishing  a 
valid  SIMS  technology  for  profiling  In  HgCdTe  was  another  goal  and  result  of 
this  work. 

Implantation  of  the  lighter  elements  Into  HgCdTe  was  emphasized  because 
Implant  depths  are  deeper  and  less  Implant  damage  Is  caused  In  the  HgCdTe 
compared  with  higher  atomic  nunber  iZ^)  elements.  Therefore,  elements  were 
selected  for  Implantation  from  the  first  four  rows  of  the  periodic  table, 
through  Br  (Z^  -  35),  except  for  In  and  Sb,  and  only  odd  elements  were 
studied  (columns  1,  3,  5  and  7).  All  odd  Z^  elements  through  Br  were  at¬ 
tempted,  that  Is,  were  Implanted  In  the  {111}  random  and  {110}  channeling 
orientations.  In  addition,  In  (Z^  ■  49)  was  Implanted  In  the  random  and 
channeling  orientations,  and  Sb  (Zj^  “SI)  was  implanted  In  the  random  orienta¬ 
tion,  because  they  are  used  as  Implantation  and  melt-grown  dopants  In  HgCdTe 
In  the  Industry. 

EXPERIMENTAL  TECHNIQUES  AND  ACCURACIES 

The  samples  of  HgCdTe  used  In  this  work  were  polished  solid  state  regrown 

{111}  and  {110}  oriented  material  annealed  for  three  or  four  weeks  at  350°C 

and  etched  to  remove  surface  damage.  All  Ions  were  implanted  at  room  tempera- 

•*2 

ture  with  current  densities  below  2  nAcm  ,  so  beam  heating  was  Insignifi¬ 
cant.  Ion  energies  are  accurate  to  within  2%,  and  Ion  fluences,  to  within 
about  10%. 

The  channeled  Implants  were  performed  by  first  aligning  the  HgCdTe  crys¬ 
tal  direction  with  the  Incident  Ion  beam  direction  to  within  0.1  degree  using 
Rutherford  backscatterlng  of  protons,  and  then  Implanting  the  selected  Ions  to 
the  desired  fluences.  The  angular  divergence  of  the  Ion  beam  during  scanning 
was  ±0,1  degree.  For  the  random  Implants,  {111}  wafers  were  tilted  7  or 
8  degrees  about  an  axis  from  which  a  flat  was  rotated  18  degrees  when  a  direc¬ 
tional  flat  was  apparent,  which  was  not  always.  In  those  cases,  the  azimuthal 
rotation  was  random. 

Secondary  Ion  mass  spectrometry  (SIMS)  measurements  were  made  using 
either  an  oxygen  primary  Ion  beam  and  positive  Ion  spectrometry,  or  a  cesium 
Ion  beam  and  negative  Ion  spectrometry.  Sputtering  rates  were  between  3  and 
10  nm/s.  The  depth  scales  of  SIMS  profiles  were  determined  from  surface 
profllometry  (Dektak)  measurements  of  the  crater  depths;  the  error  In  the 
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depth  scales  is  about  ^7%.  The  atom  densities  were  determined  by  setting  the 
integrals  of  these  unannealed  depth  profiles  equal  to  the  Ion  fluences. 

EXPERIMENTS 

Several  Implantation  studies  were  carried  out  under  this  program.  The 

first  was  an  exploratory  study  that  comprised  implanting  several  selected  Ions 

In  the  {111}  random  and  {110}  channeling  orientations  of  HgCdTe  at  150  keV  and 
1  “i  —2 

3  V  10  cm  .  The  ions  in  this  group  were  Be,  B,  F,  Na,  Al,  Cl,  and  R.  The 
{110}  channeling  depth  distributions  measured  using  SIMS  are  shown  in  Fig¬ 
ure  2-2.  Channeling  ranges  are  seen  to  vary  from  near  2  pm  for  B  to  about 
6  urn  for  F,  and  an  mdetermlned  depth  greater  than  7  um  for  Na.  The 
corresponding  {110}  random  profiles  are  shown  In  Figure  2-3.  These  {110} 
random  profiles  exhibit  significant  channeling  tails. 

Next,  H  was  studied  because  little  damage  Is  done  to  the  HgCdTe  by  this 
light  Ion  and  greater  random  ranges  can  be  achieved  than  for  heavier  Ions. 
Implanted  H  had  been  shown  to  produce  p-j  junctions  In  HgCdTe.  H  Is  a  diffi¬ 
cult  element  to  profile  using  SIMS  so  a  significant  SIMS  technology  study  was 
required.  We  were  successful  ^d  the  results  were  published  in  an  Applied 
Physics  Letter.  Subsequently,  H,  as  H  and  H,  was  channeled  In  the  {110} 

direction  of  HgCdTe  and  one  profile  is  Included  In  Figure  2-2  for  comparison 

2 

with  the  other  depth  distributions  shown  there.  The  channeled  H  profile  Is 
compared  with  a  {111}  random  profile  In  Figure  2-4.  H  Is  seen  to  have  a 
channeling  range  greater  only  than  B,  and  shallower  than  the  other  six  In  this 
group  of  elements. 

The  third  set  of  Implants  studied  comprised  the  18  elements  from  columns 

1,  3,  5,  and  7  of  the  partial  periodic  table  shown  In  Figure  2-1,  and  that  are 

potential  dopants  In  HgCdTe.  These  18  elements  were  implanted  at  200  keV  In 
the  {111}  random  and  the  {110}  channeling  orientations  of  preannealed  SSR 
HgCdTe.  {110}  channeling  profiles  are  shown  for  some  of  these  elements  In 

Figure  2-5,  to  a  depth  of  5  pm.  The  {110}  channeling  profiles  that  are  deeper 

than  5  pm  are  shown  separately  for  Individual  elements  In  Figure  2-6  for  LI, 
Figure  2-7  for  F,  Figure  2-8  for  Na  and  Figure  2-9  for  Al.  Other  channeled 
depth  distributions  are  shown  in  Figure  2-10  for  Cl,  In  Figure  2-11  for  K,  in 
Figure  2-12  for  Br,  and  in  Figure  2-13  for  In.  In  Is  of  Interest  because  of 
Its  recent  use  as  a  dopant  In  HgCdTe.  The  good  detection  limits  for  these 
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Figure  2-1.  Partial  Periodic  Table  of  the  Elements,  Shovdng  the  First  Four 
Rows  (Through  Kr)  Plus  In  and  Sb  and  the  Three  Substrate 
Components,  Hg,  Cd ,  and  Te.  The  potential  dopants  in  columns  1, 
3,  5  and  7,  which  were  also  the  ions  Implanted  and  profiled  in 
this  work,  are  outlined  in  boxes. 
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Figure  2-2.  Channeled  Depth  Distributions  for  Selected  Ions  Implanted  in  the 
{110}  direction  of  HgCdTe  at  150  keV 


RPT 60411 


25 


Figure 


Figure  2 


SAMT*  NiCCARCH  CENTIR 


-3.  Random  Depth  Distributions  for  Selected  Ions  Implanted  in 
{110}  HgCdTe  at  150  keV 


.  {Ill}  Random  and  {110}  Channeled  Depth  Distributions  for  Ions 
Implanted  in  HgCdTe  at  200  keV 
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Figure  2-5.  Channeled  Depth  Distributions  for  Selected  Ions  Implanted  in  the 
{110}  Direction  of  HgCdTe  at  200  keV  (to  5  um) 
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Figure  2-6.  (Ill)  Random  and  {110}  Channeled  Depth  Distributions  for 
Li  Implanted  in  HgCdTe  at  200  keV 


RPT60411 


HUGHES 


t  ,n17 


IONS 
200  ImV 
3x  10*5em-2 

SSR  HgCdT*  SUBSTRATE  (ANNEALED) 
RANDOM  AND 

<110>  CHANNELED  ORIENTATIONS 
ROOM  TEMPERATURE  IMPLANTS 
NO  POST  IMPLANT  ANNEAL 


< 

“  10’® 


J110)  CHANNELED 


\RANDOM 

\ 

\ 

- AMORPHOUS 

\XTAL 

I  \  I  i 


DEPTH. 


Figure  2-7*  {111}  Random  and  {110}  Channeled  Depth  Distributions  for 

F  Implanted  in  HgCdTe  at  200  keV 
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Figure  2-8.  {Ill}  Random  and  {110}  Channeled  Depth  Distributions  for 
Na  Implanted  in  HgCdTe  at  200  keV 
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Figure  2-9.  {Ill}  Random  and  {110}  Channeled  Depth  Distributions  for 
A1  Implanted  in  HgCdTe  at  200  keV 


Figure  2-10.  {Ill}  Random  and  {110}  Channeled  Depth  Distributions  for 
Cl  Implanted  in  HgCdTe  at  200  keV 
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-11.  {Ill}  Random  and  {110}  Channeled  Depth  Distributions  for 
K  Implanted  in  HgCdTe  at  200  keV 


-12.  {Ill}  Random  and  {110}  Channeled  Depth  Distributions  for 
Br  Implanted  in  HgCdTe  at  200  keV 
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Figure  2-13.  {Ill}  Random  and  {110}  Channeled  Depth  Distributions  for 
In  Implanted  In  HgCdTe  at  200  keV 

Ions  using  SIMS  allow  these  profiles  to  be  measured  even  though  the  maximum 
channeling  ranges  are  observed  only  for  atom  densities  less  than  10^^  cm“^,  or 
a  few  parts  per  billion. 

The  final  group  of  Implants  performed  under  this  program  was  a  set  of 
selected  Ions  Implanted  at  200  keV  but  at  1  x  10^^  cm”^  to  be  used  at  SBRC  for 
capacitance-voltage  profiling.  These  Ions  and  the  electrical  results  are 
reported  elsewhere  In  this  final  report.  Two  full  {110}  HgCdTe  wafers  were 
Implanted  In  the  channeling  orientation  with  each  Ion. 

The  maximum  channeling  ranges  for  Be  and  Cl  are  seen  In  Figure  2-2  to  be 
approximately  the  same  (-3.7  um)  at  150  keV,  but  a  much  larger  fraction  of  the 
implanted  Be  ions  exist  In  the  deep  channeled  portion  of  the  profile.  These 
observations  Illustrate  that  the  electronic  stopping  of  Be  and  Cl  at 
150  keV  are  nearly  identical,  accoimtlng  for  the  nearly  identical  channeling 
ranges,  but  the  small  size  of  the  Be  Ion  compared  with  the  Cl  Ion  accounts  for 
the  larger  fraction  of  more  deeply  channeled  Be  Ions.  Fewer  Ions  are  dechan- 
neled  as  they  proceed  through  the  HgCdTe  lattice.  Ms  Is  seen  to  be  an  example 
of  a  small  Ion  with  a  small  electronic  stopping,  so  Na  channels  deeply.  B  Is 

.// 


RPT 60411 


31 


SANTA  SANSAAA  Rf  SCAACM  CINTf  R 


a  small  Ion  with  a  relatively  high  value  of  electronic  stopping,  so  It  chan¬ 
nels  well  but  not  deeply.  Large  Ions  with  high  electronic  stopping  show  rela¬ 
tively  little  channeling  profile,  as  Is  the  case  for  some  of  the  higher  Z  Ions 
studied. 

We  measured  a  Rutherford  backscatterlng  7. 2Z  In  the  {110}  direc¬ 

tion  for  the  SSR  HgCdTe  material,  which  Is  a  figure  of  merit  for  the  crystal 
quality  of  the  substrate  material.  Values  of  X^in  Bahlr  and  Kallsh* 

for  the  {100}  of  RgCdTe  vary  from  20  to  30Z.  Gianneled  Ion  depth  distribu¬ 
tions  may  be  another  way  to  measure  the  crystal  quality.  In  this  case,  by 
examination  of  {110}  channeled  profiles.  The  presence  of  Hg  Interstitials, 
for  example,  would  block  the  channels  and  prevent  channeling  Ions  from  reach¬ 
ing  their  maximum  range,  or  contribute  to  an  Increased  dechanneled  portion  of 
the  channeled  depth  profile.  We  studied  this  possibility  briefly  by  measuring 
channeled  B  profiles  In  HgCdTe  substrates  prepared  under  various  annealing  and 
etching  conditions  as  Illustrated  In  Figure  2-14.  The  data  In  this  figure 
Indicate  that  the  crystal  quality  of  X  substrates  that  are  annealed  and  either 
standard  or  deep  (~20  um)  etched  are  the  same  within  experimental  accuracy  and 
are  both  better  than  unannealed  crystals.  The  long  anneal  Is  performed  to 
convert  the  n-type  as-grown  material  to  p  type. 

The  results  of  another  experiment  are  shown  In  Figure  2-14.  200-keV  B 

was  Implanted  Into  HgCdTe  In  the  {110}  channeling  orientation  of  HgCdTe  at 

doses  of  1  X  10^^,  3  x  10^^,  and  1  x  10^^  cm~^  to  determine  at  what  atom 

density  the  channeling  profile  saturates  for  B  In  HgCdTe.  This  saturation  Is 

caused  by  blocking  of  the  channels  by  Interstitials  that  result  from  Implanta- 

17  -3 

tlon  damage.  This  saturation  Is  seen  to  occur  at  about  10  cm  (at  about 

1 3  —2 

2  wm),  or  for  a  fluence  of  about  5  x  lO*-'  cm  .  This  saturation  occurs  at 

lower  fluences  for  heavier  Ions  and  the  profiles  for  all  Ions  heavier  than  B 

shown  In  this  work  Illustrate  the  saturation  density  for  {110}  channeling 

13  -2 

because  the  Implant  fluences  were  3  x  10  cm  . 

SIMS  profiles  of  Implants  Into  LPE  HgCdTe  were  less  consistent  and  had 
poorer  backgrounds  than  profiles  measured  In  the  SSR  material  used  In  this 
work.  Occasionally,  profiles  In  LPE  HgCdTe  were  ragged  as  If  there  were 
Internal  Interfaces,  Inclusions  or  complexes  within  the  LPE  layers. 
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Figure  2-lA.  Channeled  Depth  Distributions  for  200  keV  B  Implanted  Into 
Different  HgCdTe  Substrates 

sms  MEASUREMENT  TECHNOLOGY 

Some  problems  were  encountered  with  the  SmS  measurement  technology. 

High  backgrounds  of  some  elements  were  encountered  In  various  HgCdTe  materi¬ 
als,  especially  In  LPE  material.  Solid  state  regrown  (SSR)  material  and 
vertical  Bridgman  (VB)  material  exhibited  lower  background  of  most  elements. 
The  many  Isotopes  of  Cd,  Te,  and  Hg,  end  their  molecular  combinations,  and  the 
combination  of  all  of  those  with  oxygen  cause  Interferences  with  certain  Ions 
that  can  easily  be  determined.  This  problem  Is  much  less  significant  for 
negative  sms  using  Cs  because  Cd  and  Hg  do  not  yield  negative  Ions,  and  there 
are  no  doubly  Ionized  negative  Ions.  In  positive  sms,  all  of  these  Ions 
exist  In  substantial  quantities  and  their  doubly  Ionized  components  Interfere 
with  an  additional  set  of  elements.  To  learn  all  that  we  could  and  to  opti¬ 
mize  detection  sensitivity,  we  profiled  all  elements  with  both  positive  0  sms 
and  negative  Cs  SmS,  and  with  and  without  offset  voltages  In  many  cases  where 
Interferences  were  observed,  to  suppress  molecular  Interferences. 

The  systematic  measurement  of  profiles  of  Implanted  elements  using  0  and 
Cs  beams  In  a  Cameca  SIMS  Instrument  performed  under  this  program  and  a 
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parallel  IR&D  program  has  led  to  one  serious  observation.  The  substantially 
higher  energy  (14  keV)  of  the  more  massive  Cs  Ions  (133  amu)  compared  with  the 
0  Ions  (4  keV  and  16  amu)  causes  enough  greater  forward  momentum/energy  trans¬ 
fer  and  Impurity  substrate  mixing  to  cause  the  depth  distributions  to  be 
significantly  skewed  to  the  deep  side  compared  with  the  0  SIMS  profiles.  The 
widths  of  the  depth  distributions  obtained  using  a  Cs  beam  are  too  large  and 
the  deep  sides  of  the  profiles  extend  too  deeply  Into  the  HgCdTe  substrate. 
While  there  may  be  some  small  effect  for  4-keV  0,  It  Is  believed  to  be  negli¬ 
gible.  The  profiles  for  elements  beginning  at  about  the  transition  I  metals 
measured  using  Cs  SIMS  yield  values  of  the  are  slightly  too  large  (assuming 
0  profiles  are  to  be  correct),  values  of  ARp  that  are  too  large,  values  of  the 
third  moment  that  are  too  positive,  and  what  appear  to  be  channeling  tails 
when  none  or  much  smaller  ones  are  shown  using  0.  The  errors  in  depth  ap¬ 
proach  about  40Z  as  the  Z  of  the  Impurity  approaches  73,  Ta,  the  heaviest 
element  studied. 

The  profiles  reported  here  were  all  obtained  using  positive  0  SIMS  except 
In  cases  where  a  comparison  showed  Cs  profiles  not  to  be  very  distorted, 
especially  when  more  dynamic  range  was  achieved  using  Cs  SIMS,  which  was 
rarely.  This  limitation  Is  not  severe  because  all  elements  can  be  profiled 
using  0  SIMS;  only  some  elements  with  high  Ionization  potentials  yield  poor 
sensitivity  using  0  SIMS  and  these  are  primarily  the  noble  gases.  Most  of  the 
elements  In  columns  4,  5,  6,  and  7  profile  quite  well  using  0  SIMS,  in  some 
cases,  better  than  can  be  achieved  using  Cs  SIMS.  Cs  SIMS  Is  not  able  to 
measure  depth  profiles  for  the  elements  with  zero  or  very  small  electron 
affinity  (under  the  experimental  conditions).  The  sensitivity  factors  for 
SIMS  do  not  follow  the  often  assumed  rule  of  proportionality  between  Ioniza¬ 
tion  potential  or  electron  affinity  and  log  sensitivity  factor  for  Ionization 
potentials  above  9  eV  or  for  large  values  of  electron  affinity,  for  some 
substrates  (which  do  not  Include  HgCdTe).  As  part  of  this  program  and  a 
parallel  IR&D  program,  we  have  measured  sensitivity  factors  or  about  45  ele¬ 
ments  (from  H  to  Ta)  Implanted  In  HgCdTe  using  both  0  and  Cs  SIMS  and  plotted 
them  versus  Ionization  potential  and  electron  affinity.  It  Is  from  the  re¬ 
sults  of  these  studies  that  the  observations,  conclusions,  or  statements  made 
above  are  drawn. 
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In  a  few  cases,  molecular  ions  comprised  of  the  Ion  of  Interest  and  one 
from  the  substrate  or  the  primary  ion  beam  were  studied.  More  work  could  be 
done  In  this  area,  to  optimize  the  measurement  technology.  If  It  were  deemed 
significant  enough  to  expend  the  funds.  We  achieved  sufficient  detection 
sensitivity  to  obtain  channeling  profiles  for  some  but  not  all  of  the  18 
elements  of  Interest  (plus  Be),  high  backgrounds  were  often  encountered  for  P 
and  Cl. 

DISCUSSION 

Ion  energies  studied  were  150  and  200  keV,  energies  compatible  with 
commercial  Implanters  and  standard  In  the  Industry,  and  for  some  of  the  heav¬ 
ier  Ions,  250,  300,  or  400  keV.  Ion  fluences  for  channeled  profiles  were 
13  —2  13  —2 

generally  3  x  lo  cm  ,  but  1  x  10  cm  was  used  in  some  cases  for  samples 

that  were  subsequently  used  for  electrical  measurements,  for  example,  capacl- 

1 3 

tance-voltage.  Ion  fluences  studied  for  random  Implants  varied  from  1  x  10 
to  5  X  10^^  cm”^,  except  for  H,  for  which  3  x  10^^  and  3  x  10^^  cm-2  ggyg 

better  profiles  using  SIMS,  because  of  the  high  background  of  H  In  HgCdTe  and 
Interference  from  the  vacuum  ambient.  H  Is  so  light  that  the  Implant  damage 
Is  comparable  for  this  higher  fluence.  D  has  such  good  detection  sensitivity 
that  It  can  be  profiled  as  a  substitute  for  H  at  lower  fluences.  These  fluen¬ 
ces  produce  densities  toward  the  high  end  of  the  range  of  interest  for  device 
fabrication  in  the  random  portion  of  Implantation  depth  distributions,  but  the 

densities  produced  In  the  deeper  channeling  portion  of  the  profiles  (10  to 
17  —3 

10  cm  )  are  of  Interest  for  device  fabrication.  The  column  5  elements  were 
especially  difficult  to  SIMS  In  HgCdTe.  N  exists  in  the  HgCdTe  and  In  the 
vacuum  ambient.  P,  As,  and  Sb  exhibited  high  backgrounds  In  HgCdTe,  largely 
from  Interferences. 

Suitable  random  profiles  were  obtained  for  all  19  elements  studied. 
Channeling  profiles  were  not  attempted  for  Sb.  Channeling  profiles  were 
difficult  to  study  for  Qi  and  Ga  because  of  Interferences  from  doubly  ionized 
Te  and  the  oxide  of  Te ,  respectively. 

CONCLUSIONS 

We  achieved  all  of  our  goals  for  this  study,  with  the  exception  of  being 
able  to  obtain  adequate  SIMS  detection  sensitivity  to  measure  {110}  channeled 
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profiles  for  a  few  of  the  ions  selected  for  Implantation.  We  were  able  to 
profile  H  In  RgCdTe  and  obtained  both  random  and  {110}  channeled  depth  distri¬ 
butions.  We  were  able  to  Implant  and  to  SIMS  profile  the  18  elements  selected 
as  potential  dopants  In  HgCdTe  from  columns  1,  3,  5,  and  7  of  the  periodic 
table.  The  resulting  channeled  depth  distributions  and  ranges  exhibited  the 
significant  effects  of  the  Z^-dependent  electronic  stopping  seen  for  a  SI 
substrate.  We  determined  what  channeled  depths  and  associated  doping  densi¬ 
ties  can  be  achieved  for  {110}  channeling  and  the  ranges  and  other  range 
parameters  for  {111}  random  Implantation.  Maximum  {110}  channeling  ranges 

vary  from  2  to  50  um  for  200  keV  Implantation  energy,  and  saturation  atom 

18  —3  17  —3 

densities  In  the  deep  channeled  region  vary  from  10  cm  for  H  and  10  cm 
for  B  to  10^^  cm~^  for  the  heavier  Ions  that  channel  deeply.  We  determined 
that  preannealed  HgCdTe  material  has  a  better  crystal  quality  than  unannealed 
material  (via  the  degree  of  B  Ion  channeling),  and  that  material  subjected  to 
standard  and  deep  etching  are  of  equivalent  quality,  as  measured  by  this 
criterion.  We  measured  a  Rutherford  backseat terlng  of  7.2%  for  our 

HgCdTe  material  In  the  {110}  direction.  We  found  that  LPE  HgCdTe  Is  of  poorer 
quality  than  SSR  HgCdTe;  we  observed  higher  Impurity  backgrounds,  poorer 
consistency,  and  apparent  burled  Interfaces  and  Inclusion.  We  implanted  a 
number  of  full  wafers  of  HgCdTe  with  selected  Ions  In  the  {110}  channeling 
orientation  for  electrical  measurements  to  be  made  at  SBRC.  We  learned  much 
about  measuring  depth  distributions  In  HgCdTe  using  SIMS  and  established  a 
SIMS  technology  In  HgCdTe.  This  technology  requires  the  use  of  0  SIMS  rather 
than  Cs  SIMS  for  profiling  elements  heavier  than  about  Fe.  We  found  that  SIMS 
sensitivity  factors  do  not  follow  the  rule  of  proportionality  between  log 
sensitivity  factor  and  Ionization  potential  or  electron  affinity. 
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2.2  DOPING  PROFILES  OF  CHANNELED  ION  IMPLANTATION  INTO  Hgg  yCdg  3Te 
Introduction 

This  task  pursued  obtaining  higher  resolution  doping  profiles  of  HgCdTe 
wafers  implanted  In  the  channeling  {110}  direction.  A  comparison  is  made 
between  the  damage-induced  donor  levels  and  the  implanted  species  level.  This 
provides  some  Insight  into  the  damage  mechanism  and  separation  of  nuclear 
versus  electronic  stopping. 

Sample  Preparation  and  Experiment 

An  Ingot  of  X  »  0.3  HgCdTe  grown  by  Solid  State  Recrystalllzatlon,  was 

oriented  and  sliced  to  expose  the  {110}  surface.  These  wafers  were  then 

polished  to  remove  saw  damage  and  Implanted  with  various  Ions  at  200  keV  to  a 
1 3  —7 

dose  of  1  X  10  cm  .  Ions  studied  were  H,  Be,  0,  Cl,  P,  SI  and  In.  A  dip 

etch  was  used  to  obtain  a  tapered  surface.  Etching  while  piUllng  the  wafer 

out  at  a  constant  rate  results  In  the  far  edge,  which  sees  the  etch  longest, 

having  the  most  material  removed.  The  nearest  edge  transitions  smoothly  to 

the  original  surface.  The  final  amount  of  material  removed  from  the  far  edge 

was  2.75  ±  0.25  pm  (the  wafers  are  about  1  cm  In  diameter)  for  H,  Be  and  Cl, 

and  5.0  ±  0.5  ym  for  0,  P,  Si  and  In.  The  wafers  then  received  1500A  of 

PHOTOX"  SIO2  and  deposition  of  400A  Ti  followed  by  4000A  Au  to  form  capacitors 

—4  2 

with  and  area  of  4  x  10  cm  and  a  center-to-center  capacitor  separation  of 

145  pm.  This  gives  a  change  In  depth  probed  per  capacitor  of  4.5  x  10  pm 
_2 

and  9.1  X  10  pm  for  the  two  etch  depths.  Capacitance-voltage  data  was  then 

taken  at  1  MHz  to  yield  versus  depth.  The  background  donor  level  due 

15  —3 

to  residual  impurities  Is  ~3  x  10  cm  . 

Results 

Figures  2-15  through  2-21  show  depth  profiles  of  |Njj-N^|  from  C-V  com¬ 
pared  to  the  Implant  species  concentration  from  SIMS  for  H,  Be,  0,  P,  Cl,  Si 
and  In,  respectively.  The  first  thing  to  be  noted  as  the  transition  Is  made 
toward  heavier,  larger  atoms  Is  the  decrease  in  the  channeled  portion  of  the 
SIMS  curve.  All  samples  have  some  random  component  due  to  the  Impossibility 
of  all  Incoming  Ions  being  channeled.  For  larger  Ions  more  collide  with  the 
lattice  at  the  surface  rather  than  being  channeled.  There  are  significant 
exceptions  to  this,  however.  Both  Si  and  0  show  minimal  channeling  while  H 
and  Be  (which  are  lighter)  and  Cl  and  P  (which  are  heavier)  show  a  fair  amount 
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-15.  Atomic  and  Donor  Concentration  Depth  Profiles  for  Hydrogen 


Figure  2-16.  Atomic  and  Donor  Concentration  Depth  Profiles  for  Beryllium 
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2-17.  Atomic  and  Donor  Concentration  Depth  Profiles  for  Oxygen 


-18.  Atomic  and  Donor  Concentration  Depth  Profiles  for  Silicon 
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19.  Atomic  and  Donor  Concentration  Depth  Profiles  for  Phosphorus 


2-20.  Atomic  and  Donor  Concentration  Depth  Profiles  for  Chlorine 
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Figure  2-21.  Atomic  and  Donor  Concentration  Depth  Profiles  for  Indium 

of  channeling.  A  similar  effect  is  seen  in  silicon.^  The  "Zl”  channeling 

probability  is  dependent  on  the  radial  electron  density  which  varies  somewhat 

2  3 

cyclicly  with  mass.  '  It  appears  that  the  same  dependence  holds  in  HgCdTe. 
The  absolute  values  of  Zl  will  of  course  be  different  due  to  different  channel 
size,  bond  strength  and  atoms  present  in  the  lattice. 

The  next  point  of  interest  is  a  comparison  of  level  and  peak  position  in 
the  "random"  portion  of  the  profile.  In  Si  and  GaAs  the  ratio  of  the  damage 
level  to  ion  concentration  goes  up  with  increased  mass.  This  is  due  to  heav¬ 
ier  ions  carrying  less  of  the  recoil.  The  lattice  absorbs  more  of  this  energy 
propagating  a  higher  damage  level.  This  effect  is  not  seen  in  the  HgCdTe  sam¬ 
ples.  The  ratio  of  damage  to  ions  was  10  ±  1  for  all  ions  except  H  and  Cl. 

It  is  very  likely  that  the  reason  for  this  is  a  "self-annealing"  effect  in  the 
crystal.  This  effect  has  been  seen  on  a  microsopic  scale  by  T.  James^  in  a 
comparison  of  Implants  performed  at  300K  and  77K.  The  exception  of  H,  where 
the  peak  damage  level  is  a  factor  of  10  below  the  ion  concentration  is  to  be 
expected  due  to  its  small  mass  and  hence  very  little  energy  transferred  to  the 
lattice.  Cl,  on  the  other  hand,  shows  a  ratio  of  100.  This  is  substantially 
higher  than  the  other  ions  and  has  been  seen  on  previous  samples. 
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The  shape  of  the  damage  profile  follows  that  of  the  random  portion  of  the 
Ion  profile.  In  the  case  of  Be  the  damage  level  actually  drops  below  the  Ion 
level  In  the  channeled  region  of  the  sample.  The  Indication  then  Is  that 
nuclear  stopping  Is  responsible  for  Induced  damage  levels  and  that  the  weaker 
electronic  stopping  forces  do  not  distort  the  lattice  enough  to  produce  elec¬ 
trically  active  sites.  The  exceptions,  again,  are  H  and  Cl.  For  H,  the 
damage  level  Is  fairly  constant  and  below  the  peak  Ion  concentration.  This 
can  be  explained  by  Its  small  mass  as  mentioned  earlier.  In  the  case  of  Cl, 
however,  there  Is  a  "shoulder”  In  the  damage  profile,  beyond  the  random  por¬ 
tion  of  the  SIMS  curve,  which  appears  to  follow  the  channeled  Ion  distribu¬ 
tion. 

Conclusion 

Several  trends  are  apparent  In  this  data.  Damage- Induced  donor  levels 
are  Indicated  to  be  caused  by  nuclear  rather  than  electronic  stopping.  The 
level  of  this  damage  may  initially  scale  with  Increased  mass,  however  a  self¬ 
annealing  effect  is  apparent,  which  results  In  all  Ions  resulting  In  approxi¬ 
mately  the  same  damage  level.  Hydrogen  is  an  exception  due  to  It's  small 
mass.  The  behavior  of  Cl,  which  does  not  follow  these  trends,  will  require 
further  Investigation.  Beyond  the  region  of  electrical  damage  and  nuclear 
stopping  the  depth  to  which  Ions  are  channeled  depends  roughly  on  the  radial 
electron  distribution. 
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Our  investigation  into  the  metal-HgCdTe  Interface  continues  with  new  data 
presented  on  the  Cu-HgCdTe  Interface  formation  and  a  parallel  study  of  the  Cu- 
CdTe  Interface.  These  experiments  were  performed  in  UHV  on  clean  cleaved 
samples  and  analyzed  using  PES.  The  role  of  the  weak  HgTe  bonding  in  the 
lattice  is  emphasized  in  explaining  the  interface  morphology  for  the  Cu-HgCdTe 
and  Cu-CdTe  systems.  The  differences  and  similarities  between  the  resultant 
morphologies  of  these  metal-semiconductor  interfaces  and  the  Ag  and  A1  metal- 
HgCdTe  systems  are  also  discussed,  and  the  role  of  the  Hg  bonding  in  each  is 
elucidated.  Section  3.1  of  this  report  contains  a  full  discussion  of  these 
results. 


The  calculation  of  the  band  structure  of  CdTe  along  the  [110]  direction 
is  reported  in  Section  3.2.  The  calculation  was  performed  using  the  nonlocal 
pseudopotential  of  Chellkowsky  and  Cohen  (see  references,  Section  2.2)  and 
gives  insight  into  the  role  of  the  actual  final  states  Involved  in  the  normal 
emission  process  in  Influencing  some  of  the  features  of  the  observed  spec¬ 
tra.  Much  of  the  observed  spectra  agree  with  the  approximation  that  the  final 
states  are  free  electron-like  and  therefore  are  parabolic,  but  certain  peaks 
in  the  data  appear  to  be  influenced  by  a  mixing  of  the  parabolic  bands  that 
takes  place  in  the  crystal  potential.  The  optical  transition  matrix  elements 
for  Interband  transitions  were  determined  using  the  calculated  wave  functions 
of  each  band.  These  matrix  elements  are  a  major  contribution  in  determining 
the  intensity  of  spectral  peaks  and  their  polarization  dependence  agreed  with 
the  measured  in  experiment. 


Section  3.3  contains  new  information  extracted  from  angle- re solved  photo- 
emission  (ARPES)  experiments  previously  conducted.  Attempts  at  the  mapping  of 
the  electronic  valence  bands  along  the  [100]  and  [111]  directions  from  the 
data  obtained  in  these  experiments  have  been  undertaken.  Due  to  the  experi¬ 
mental  complexity  resulting  from  off-normal  emission  in  obtaining  this  data, 
the  mapping  of  the  bands  in  these  cases  is  very  difficult  and  nontrivial. 

Preliminary  results  indicate  that  the  expected  theoretical  dispersion  of  these 

bands  is  found  in  these  angle-resolved  spectra.  The  application  of  these 

ARPES  results  into  the  investigation  of  the  alloying  effect  in  HgCdTe  is  also  < 

discussed. 
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At  the  beginning  of  this  study,  bonding  of  the  cation  In  the  lattice  was  j 

I 

Imperfectly  understood  and,  thus,  the  consequences  of  this  were  not  fully  i 

appreciated.  As  we  all  know,  HgCdTe  crystallizes  In  the  zlncblende  crystal 
structure  which  can  be  thought  of  as  two  Interpenetrating  face-centered  cubic 
sublattices,  one  displaced  from  the  other  by  1/4  [111].  Te  occupies  the  sites 
of  one  of  the  sublattices,  and  Hg  and  Cd  occupy  sites  on  the  other.  The 
resulting  bond  configuration  Is  one  of  tetrahedral  symmetry  where  the  outer 
valence  electrons  on  the  anion  and  cation  hybridize  to  form  sp  -hybrids  to 
take  place  In  covalent  bonding.  When  one  does  theoretical  calculations  to 
determine  the  electronic  valence  bands,  and  hence  derive  many  important  param¬ 
eters  such  as  electronic  transport  and  ease  of  dislocation  formation,  one  must 
employ  the  crystal  potential  of  the  lattice.  To  do  this  calculation  for  most 
ternary  alloy  semiconductors,  for  example  an  alloy  system  in  which  there 
exists  two  different  cations  located  on  the  cation  sublattice  a  composition- 
ally  weighted  average  of  the  cation  potentials  Is  employed  on  the  cation 
sublattice.  This  virtual  crystal  approximation  (VGA)  typically  yields  theo¬ 
retical  parameters  that  are  In  close  agreement  with  experimentally  determined 
values.  As  one  can  see,  there  are  no  provisions  In  this  theoretical  framework 
to  account  for  any  differences  In  the  bonding  between  the  anion  and  the  two 
different  nearest  neighbor  cations. 

In  HgCdTe,  some  of  the  behavior  of  the  crystal  suggests  a  very  periodic 
lattice,  whereas  other  observed  phenomena  was  reflective  of  a  more  random 
lattice  where  the  bonding  between  Hg  and  Te  Is  distinguishable  from  that  of 
CdTe.  The  extremely  high  electron  mobilities  in  the  conduction  band  minimum 
(CBM)  Is  Indicative  of  highly  nonlocallzed  electronic  wave  functions  charac¬ 
teristic  of  a  periodic  rather  than  random  lattice.  The  ease  of  breaking  the 
HgTe  bond  In  the  alloy  resulting  In  the  migration  of  Hg  to  the  surface  and  Its 
subsequent  evaporation  is  suggestive  of  a  different  bonding  for  Hg  and  Cd  in 
the  alloy,  and  this  difference  should  result  In  a  different  potential  for  Hg 
and  Cd  and  hence  be  reflected  In  the  electronic  structure. 

The  techniques  of  FES  and  ARPES  available  to  our  research  group  were  able 
to  probe  the  valence  electronic  states  that  are  both  highly  localized  and 
nonlocallzed,  and  therefore  revealed  several  very  Interesting  features 
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concerning  the  bonding  In  the  lattice.  When  one  observes  the  nonlocallzed 
states  at  the  VBM  In  the  region  of  r  In  the  Brlllloun  zone,  which  are  com¬ 
prised  mostly  of  Te-p  states,  across  the  composition  range  from  HgTe  to  CdTe, 
these  states  remain  fixed  In  energy  and  hence  are  fairly  Insensitive  to  compo¬ 
sitional  changes.  The  CBM  Is  predominately  cation  s-derlved  and  since  It  also 

occurs  at  T ,  these  very  nonlocallzed  states  sample  an  average  potential  of  the 

2 

entire  lattice.  Since  the  Hg6s  electrons  are  bound  1.4  eV  deeper  than  the 
CdSs  ,  substituting  Ilg  for  Cd  In  CdTe  will  lower  the  conduction  band  and  hence 
this  CBM  moves  linearly  with  composition  as  one  would  expect  In  the  approxima¬ 
tion  of  a  virtual  crystal.  This  advantageous  movement  of  the  CBM  while  the 
VBM  remains  fixed  Is  idiat  gives  HgCdTe  Its  very  attractive  characteristic  of 
bandgap  tunablllty. 

Insight  Into  a  very  Interesting  feature  of  the  bonding  In  the  alloy  Is 
gained  of  one  looks  at  the  electronic  valence  band  states  located  between  4 
and  6  eV  below  the  VBM.  These  states  are  located  at  the  zone  boundary  and 
therefore  are  very  localized  and  hence  sensitive  to  Hg  and  Cd  percentage.  Two 
distinct  peaks  are  seen  that  are  fixed  in  energy  but  change  In  Intensity  as 
the  composition  changes.  The  deeper  bound  peak  can  be  attributed  to  Hg 
s-states  at  5.4  eV  and  the  lesser  bound  peak  at  4.6  eV  Is  attributed  to  Cd 
8-states.  This  represents  a  deviation  from  virtual  crystal  behavior  because 
one  would  expect  these  states  to  change  linearly  with  composition  In  the 
framework  of  the  VCA.  This  selective  breakdown  of  the  VCA  reveals  that  the 
HgTe  bonds  are  different  from  the  CdTe  In  the  lattice,  and  this  unusual  behav¬ 
ior  Indicates  that  a  more  sophisticated  picture  of  the  bonding  In  this  alloy 
must  be  found.  One  must  take  Into  account  the  difference  In  the  cation  poten¬ 
tials  In  evaluating  the  alloy  electronic  structure.  The  coherent  potential 

a 

approximation  (CPA)  takes  this  Into  account,  and  calculated  density  of  states 
using  the  CPA  give  much  better  agreement  to  experimental  data  than  using  the 
VCA. 

The  selective  breakdown  of  the  VCA  also  gives  Insight  Into  many  of  the 
observed  mechanical  deficiencies  of  HgCdTe,  as  well  as  the  point  defect  equi¬ 
libria  which  gives  rise  to  acceptor-type  behavior.  The  deeper  bond  Hg 
s-valence  bonding  states  are  retained  nearer  to  the  Hg  nucleus  than  chat  of 
the  lighter  Cd ,  and  hence  the  %  does  not  produce  as  strong  of  a  covalent  bond 
to  Te  than  does  Cd.  This  results  In  a  HgTe  bond  chat  Is  as  much  as  30%  weaker 
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than  the  CdTe  in  the  HgCdTe  alloy.  There  is  no  other  semiconductor  with  such 
a  large  difference  In  the  pseudobinary  constituent  bonds.  Further  Information 
Is  gathered  from  the  UPS  spectra  of  different  alloy  compositions  by  looking  at 
the  behavior  of  the  Cd4d  and  HgSd  core  levels.  A  key  finding  In  our  work  was 
that  upon  adding  Hg  to  CdTe,  the  Cd4d  becomes  more  deeply  bound,  whereas  the 
addition  of  Cd  to  %Te  decreases  the  HgSd  binding  energy.  Therefore  the  UgTe 
bond  weakens  further  and  the  CdTe  strengthens  In  the  alloy  as  compared  to 
their  energy  In  the  respective  binaries.  Consequently,  adding  Gd  to  HgTe 
exacerbates  the  Instability  of  the  alloy,  and  can  give  rise  to  greater  densi¬ 
ties  of  point  defect  such  as  Hg  vacancies,  and  the  breaking  of  the  HgTe  can 
result  In  Hg  out-dlf fusion  from  the  lattice.  This  has  severe  ramifications 
Into  the  structural  Integrity  of  the  alloy,  and  the  history  of  the  sample 
during  growth  and  annealing  Is  very  critical.  In  all  of  our  studies  of  HgCdTe 
as  listed  In  Section  4.1,  this  bonding  nature  of  the  alloy  Is  manifested  In 
studies  of  the  Interfaces,  surfaces,  and  bulk  properties  of  HgCdTe. 

Experimental  determination  of  the  electronic  valence  band  structure  along 
the  major  symmetry  directions  [110],  [100],  and  [111]  using  ARPES  has  also 
been  performed  under  this  contract.  Accurate  knowledge  of  the  band  structure 
of  a  material  enables  many  materials  properties  to  be  estimated.  Properties 
ranging  from  effective  mass  In  the  conduction  band  to  vacancy  formation  ener¬ 
gies  may  be  estimated  from  a  detailed  picture  of  the  bands.  Accurate  knowl¬ 
edge  of  the  electronic  bands  Is  also  essential  when  deriving  the  transport 
properties  of  superlattices,  which  are  currently  of  great  Interest  in  HgCdTe 
technology.  The  theoretical  properties  calculated  for  these  materials  are  not 
as  precise  without  this  kind  of  knowledge  derived  from  our  study. 

In  the  ARPES  studies  performed  on  this  contract,  the  acquired  spectral 
data  points  were  mapped  along  the  primary  directions  In  reciprocal  space,  and 
the  experimental  points  were  compared  to  theoretical  calculations.^  The 
Initial  state  bands  are  determined  by  assuming  a  free  electron  final  state  in 
the  conduction  band.  This  first  order  approximation  holds  for  some  of  the 
bands,  or  rather,  yields  Initial  state  bands  that  agree  well  with  the  theoret¬ 
ically  predicted  resiilts,  but  for  others  the  assumption  of  a  free  electron 
final  state  and  folding  back  Into  the  first  zone  of  the  primary  reciprocal 
lattice  vector  (e.g.  G  -  4ifa[110]  In  the  case  of  [110]  normal  emission)  yield 
Inaccurate  Initial  state  values.  In  Section  2.2  of  this  report,  a  more 
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sophisticated  method  of  determining  the  initial  states  is  presented.  Nonlocal 
pseudopotential  calculations  of  the  conduction  bands  up  to  30  eV  above  VBM 
indicate  a  mixing  of  the  parabolic  bands  in  the  crystal  potential  and  the 
incorporation  of  other  reciprocal  lattice  vectors  in  the  electron  scattering 
provide  for  a  more  accurate  initial  state-band  determination.  Another  utility 
of  our  ARPES  work  is  that  it  allows  one  to  determine  more  subtle  features  in 
the  valence  band  such  as  alloy  broadening.  Closely  related  is  the  ability  to 
determine  the  effect  of  Hg  loss  by  experimentally  examining  ARPES  as  a  func¬ 
tion  of  Hg  loss.  The  details  of  this  experiment  are  given  in  Section  2.3  of 
this  report. 

Native  Oxidation  of  HgCdTe 

One  of  the  primary  problems  one  must  solve  to  Increase  the  reliability, 
optimize  the  processing  technology,  and  improve  device  yield  is  to  Improve  the 
passivation  of  the  HgCdTe  surface.  The  oxidation  of  this  ternary  alloy  system 
is  very  complex  and  the  resultant  equilibrium  oxides  are  determined  by  the 
intensive  thermodynamic  state  variables  that  the  system  is  subjected  to  during 
growth  of  the  oxide.  There  are  many  equilibrium  or  metastable  oxides  that  may 
be  theoretically  produced,  but  most  are  subject  to  decomposition  even  at 
relatively  low  temperatures,  or  upon  the  deposition  of  other  passivation 
layers.  This  unreliability  in  the  structural  and  hence  dielectric  properties 
cast  doubt  as  to  the  benefit  of  producing  native  oxides  at  all  on  HgCdTe 
surfaces.  The  very  low  sticking  coefficient  for  oxygen  on  HgCdTe  also  renders 
this  a  very  difficult  process  to  achieve.  We  suggested  very  early  on  in  this 
contract  that  it  might  be  beneficial  to  minimize  the  native  oxides  of  HgCdTe 
and  deposit  a  layer  of  an  external  oxide,  e.g.  SIO2,  as  gently  as  possible  to 
passlfy  the  surface. 

To  understand  more  fundamentally  the  kinetics  of  oxide  formation  and  the 
initial  stages  of  uptake  of  oxygen  on  the  surface  of  HgCdTe,  PES  was  employed 
as  a  sensitive  surface  probe  to  study  this  interaction.  Because  of  the  dif¬ 
ference  in  the  cation  bonding  in  the  lattice,  a  sensitive  study  of  the  oxygen 
chemistry  with  HgCdTe  should  also  reveal  a  difference  in  the  kinetics  and 
activity  of  the  system  as  a  function  of  composition.  There  were  two  primary 
areas  of  research  in  this  study,  the  identity  and  character  of  the  oxygen 
species  responsible  (e.g.  0,  O2,  etc),  and  the  rate  of  oxygen  uptake  on  the 
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surface.  Initial  experiments  Indicated  that  molecular  oxygen  that  was  unex¬ 
cited  did  not  stick  appreciably  to  the  HgCdTe  surface  over  a  period  of  many 
hours.  When  the  Ion  gauge  filament  (electron-filament)  was  turned  on  In  line 
of  sight  of  the  sample  oxygen  uptake  was  stimulated.  The  exact  nature  of  the 
excited  oxygen  species  Is  not  known »  but  It  most  likely  appears  to  be  molecu¬ 
lar  oxygen  which  may  be  vlbratlonally  excited. 

The  complex  oxygen  chemistry  of  HgCdTe  tends  to  support  the  bonding 
theories  discussed  In  the  previous  section.  It  was  found  that  the  rate  of 
oxide  production  of  HgQ. 69^‘^0. 31^®  HgTe ,  therefore  the 

rate  of  growth  scales  with  the  nunber  of  HgTe  bonds  in  the  alloy.  This  oxide 
growth  Is  accompanied  by  a  release  of  Hg  from  the  lattice  and  an  Incorporation 
of  the  Cd  In  the  oxide.  The  free  Te  to  oxidizes  once  the  bond  to  the  Hg  has 
been  broken.  The  oxidation  rate  of  CdTe  was  faster  than  that  of  the  alloy, 
indicating  a  lowering  of  the  activity  of  the  CdTe  sites  In  the  alloy.  One  of 
the  major  conclusions  than  can  be  drawn  from  our  oxidation  studies  Is  that 
even  the  first  one  or  two  atomic  layers  alter  the  HgCdTe  composition  at  the 
surface,  which  could  have  deleterious  consequences  for  any  devices  passivated 
with  native  oxide. 

Mechanical  Defects 

Observations  Into  the  change  of  electronic  behavior  of  HgCdTe  upon  me¬ 
chanical  stress  and  the  question  of  mechanical  stability  of  the  lattice 
prompted  our  TEM  studies  Into  the  mechanical  defect  structure  of  HgCdTe.  Also 
desirable  in  our  studies  was  Information  into  the  heterojunction  widths  be¬ 
tween  the  substrate  CdTe  and  HgCdTe  LPE  layers.  Three  primary  studies  were 
conducted  on  HgCdTe  under  this  contract.  The  determination  of  the  defects  In 
the  lattice  of  SSR  p-type  HgCdTe  giving  rise  to  etch  pits  produced  by  the 
Pollsar^  etch-2  was  the  first  study.  The  defect  structure  at  a  cleaved  HgCdTe 
surface  was  also  Investigated.  The  desire  to  observe  at  LPE  layers  In  cross 
section  prompted  another  study  Into  the  difference  of  Ion  milling  Induced 
defects  In  l^CdTe  and  CdTe. 

A  traditional  method  of  quickly  determining  the  quality  of  a  single 
crystal  of  known  orientation  Is  to  etch  the  surface  with  a  suitable  defect 
etchant  to  produce  etch  pits  on  the  surface.  The  concentration  and  distribu¬ 
tion  of  these  etch  pits  yields  valuable  information  into  the  defect  structure 
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of  the  lattice.  Because  of  the  ambiguity  that  may  exist  as  to  the  origin  of 
the  defects  giving  rise  to  the  etch  pits,  we  undertook  a  study  using  TEM  to 
identify  the  defects.  In  a  predominate  number  of  etch  pits,  Te  precipitates 
were  found  to  most  likely  have  given  rise  to  the  etch  pit.  It  is  important  to 
emphasize  that  the  material  study  was  SSR  p-type  HgCdTe,  and  the  [111]  A-face 
(HgCd)  give  rise  to  the  etch  pits  upon  etching  with  the  Polisar^  etch-2.  Near 
the  edge  of  the  ingot  the  etch  pits  lie  in  an  array  strongly  reminiscent  of 
the  dendritic  arms  in  the  original  cast  ingot.  These  observations  suggest 
that  the  SSR  ingot  is  strongly  influenced  by  the  casting  stage. 

Another  study  of  great  importance  was  to  understand  the  defect  structure 
induced  on  the  [110]  face  upon  cleaving  because  this  is  the  surface  studied  in 
PES  and  hence  is  of  great  theoretical  importance.  The  face  was  found  to  be 
rough  on  an  optical  scale,  and  the  fracture  surface  exhibited  numerous  steps 
and  the  fracture  appeared  to  be  brittle  in  some  regions.  Upon  examining  the 
surface  with  TEM,  it  was  found  that  at  the  portion  of  the  surface  that  a  step 
appeared,  a  high  density  of  dislocations  had  been  nucleated  at  the  step.  The 
crack  appeared  to  proceed  slightly  off  the  ideal  cleave  face,  arrest  with  a 
corresponding  nucleation  of  many  dislocations,  and  then  proceed.  Thus,  the 
cleavage  face  was  found  to  be  much  more  defective  than  the  bulk.  This  might 
explain  why  the  surface  often  converts  to  n-type  upon  cleaving,  and  future 
experiments  looking  into  this  are  planned  at  Stanford. 

Another  important  mechanical  feature  of  HgCdTe  is  the  abruptness  of  the 
transition  region  between  the  CdTe  substrate  and  the  LPE  HgCdTe.  These  inter¬ 
faces  are  difficult  to  image  with  the  TEM  because  it  is  difficult  to  make  a 
TEM  sample  without  introducing  artifacts.  We  employed  ion  milling  to  thin  the 
material,  and  although  we  introduced  many  artifacts,  the  distribution  of  the 
ion  milling  induced  defects  in  the  CdTe/HgCdTe  heterojunction  was  very  sur¬ 
prising.  Detailed  analysis  of  the  defects  in  the  HgCdTe  and  CdTe  revealed  the 
probable  reason  for  a  2  to  3  order  of  magnitude  higher  defect  density  in  the 
CdTe;  the  stacking  fault  widths  in  the  CdTe  were  much  wider  than  in  HgCdTe. 
This  lends  support  to  the  theoretical  predictions^  of  a  weak  HgTe  bond  and  an 
increase  in  the  metalllclty  of  this  bond  because  the  stacking  fault  should  be 
narrower  for  less  directional  or  rigid  bonds  and  the  number  of  HK  vacancies 
greater  than  a  weaker  HgTe  bond.  Both  of  these  characteristics  should 
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increase  the  rate  of  recovery  of  the  hgCdTe  as  compared  to  the  CdTe,  which  was 
observed  experimentally. 

Metal  Overlayers 

A  relatively  neglected  Interface  up  to  this  point  In  HgCdTe  Is  that  of 
the  metal/HgCdTe  system.  Most  processors  of  HgCdTe  devices  tend  not  to  worry 
about  metal  contacts  because  they  are  for  the  most  part  ohmic  contacts  on 
n-type  material.  However,  ohmic  contacts  on  n  as  well  as  p-type  materials  are 
Important.  Similarly,  the  ability  to  produce  reliable  Schottky  barriers  is 
Important  for  optimizing  devices.  But  In  light  of  the  chemistry  that  occurs 
between  HgCdTe  and  even  two  monolayers  of  oxygen,  the  stoichiometry  of  the 
HgCdTe  In  the  vicinity  of  the  surface  may  be  drastically  altered  upon  deposi¬ 
tion  of  a  reactive  or  even  unreactlve  metallic  species,  which  could  produce 
deleterious  effects  on  the  electronic  bands  In  the  vicinity  of  the  surface. 

We  have  studied  the  chemistry  of  the  HgCdTe  surface  upon  deposition  of  A1 ,  Ag , 
and  Cu.  The  A1  Is  predicted  thermodynamically  to  react  strongly  with  the  Te, 
and  the  Ag  should  not  react,  and  Cu  should  fall  between  them  In  reactivity 
according  to  heats  of  formation  data.  In  addition  to  the  chemistry  of  these 
surfaces,  pinning  of  the  ferml  level  was  also  monitored  as  a  function  of 
deposition. 

'Ihere  are  several  Interesting  observations  that  we  have  made  Into  the 
characteristics  of  each  Interface.  Wie  have  found  that  the  reactive  A1  case 
drastically  depletes  the  surface  of  Hg  as  the  first  few  monolayers  are  depos¬ 
ited.  This  Is  not  the  case  for  Ag  or  Cu.  In  the  Ag  case,  the  Ag  Is  seen  to 
replace  the  in  the  top  approximate  lOOOA,  it  therefore  diffuses  very  deeply 
into  the  bulk.  For  the  Al/HgCdTe  Interface,  the  first  3  to  5  monolayers  are 
Incorporated  Into  an  AlTe  reacted  layer,  with  further  A1  deposition  being 
Incorporated  as  metallic  Al.  The  deposition  of  Cu  onto  HgCdTe  and  CdTe  re¬ 
veals  an  enhanced  Te  and  Cd  Intermixing  into  the  Cu  for  HgCdTe  compared  to 
CdTe.  These  systems  were  discussed  In  Section  2.1  of  this  report.  The  ferml 
level  Is  also  pinned  for  each  case.  For  Ag  on  HgCdTe  the  bands  bend  up  0.05 
to  0.10  eV  which  indicates  a  lower  binding  energy  for  this  case.  In  the  Al  on 
HgCdTe  case,  all  core  level  and  VB  features  shift  0. 3  ±  0.05  eV  to  higher 
binding  energy,  l.e.  the  ferml  level  moves  upward  from  the  VBM  energy.  The  Cu 
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case  Is  still  beln;;  investigated.  We  have  shown  unique  behavior  in  HgCdTe 
metal  contacts,  much  more  remains  to  be  done  to  elucidate  these  findings. 


Cd  Loss  in  BrM  Etched  Surfaces 

In  extending  our  studies  to  HgCdTe  surfaces  of  practical  Interest,  we 
have  observed  with  PES  profound  Cd  loss  upon  etching  these  surfaces  with 
BrM.  Losses  between  5%  and  20%  have  routinely  been  observed.  These  techno¬ 
logically  Important  surfaces  (e.g.  the  [111]  face  for  LPE  and  MBE)  must  be 
cleaned  before  the  initial  processing  steps  can  be  performed  to  build  up  the 
device  structure.  That  we  see  Cd  loss  upon  etching  HgCdTe  with  BrM  is  not  too 
surprising  in  light  of  the  etching  rate  experiments  performed  by  Talasek,^  but 
the  magnitude  of  this  loss  suggests  that  the  material's  crystalline  structure 
and  electrical  properties  could  be  greatly  affected,  the  significance  of  our 
studies  outlined  in  the  next  section  In  studying  cation  loss,  and  Its  effect 
on  the  electronic  structure,  Is  underlined  by  this  finding  on  practical  sur¬ 
faces. 

Electron  Beam  Induced  Desorption 


Because  HgCdTe  is  seen  to  lose  the  cation  species  when  subjected  to 
conditions  identical  to  those  seen  during  processing,  it  is  very  important  to 
study  the  effect  of  this  loss  of  stoichiometry  on  the  electronic  structure  of 
the  alloy.  Such  studies  are  also  of  fundamental  Importance.  We  have  men¬ 
tioned  many  experiments  in  this  report  that  result  in  a  loss  of  Hg  from  the 
lattice.  Depositing  metals  on  the  surface  or  reacting  the  surface  with  ex¬ 
cited  oxygen  cause  the  surface  region  of  HgCdTe  to  be  depleted  of  Hg  to  vary¬ 
ing  depths.  Etching  free  HgCdTe  surfaces  with  bromine  methanol  (BrM)  results 
in  profound  Cd  loss  from  the  surface  region.  A  very  powerful  study  may  be 
conducted  if  one  of  the  cations  can  be  removed  in  a  well-controlled  manner  and 
the  resulting  surface  probed  using  ARPES  to  reveal  the  subtle  effects  of 
cation  loss  on  the  electronic  structure  of  the  altered  HgCdTe. 

Because  we  can  remove  Hg  from  the  lattice  by  impinging  the  surface  of 
HgCdTe  with  an  electron  beam,  and  the  electron  beam  parameters  can  be  con¬ 
trolled  precisely,  we  are  able  to  surgically  remove  Hg  and  study  the  resultant 
electronic  structure.  We  can  also  monitor  the  crystallinity  of  the  lattice 
with  low  energy  electron  diffraction  (LEED)  to  give  us  critical  information  as 
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to  the  structural  nature  of  the  lattice  after  profound  cation  loss.  At  the 
Initial  stage  of  Hg  loss  (approximately  3%) ,  the  surface  structural  crystal¬ 
linity  and  electronic  structure  are  Intact,  and  the  surface  becomes  degenerate 
p-type.  These  parameters  remain  Intact  up  to  about  9%  loss,  but  at  greater  Hg 
loss  the  surface  crystal  and  electronic  structure  gradually  lose  the  crystal¬ 
line  character  of  a  covalently  bonded  lattice.  The  selection  rules  dictated 
by  a  crystalline  solid  are  lost  and  the  valence  band  loses  Its  dispersion. 
These  observations  reveal  that  the  attractive  electronic  properties  of  HgCdTe 
can  be  drastically  altered  upon  processing  the  alloy  when  making  devices. 

Imnllcatlons  Into  Future  Research 


HgCdTe  Is  of  extreme  Importance  because  Its  bandgap  Is  tuned  by  compo¬ 
sition  to  any  value  from  0  to  1.5  eV.  This  has  given  It  great  potential  as  an 
Infrared  (IR)  detector.  However,  great  difficulties  have  arisen  In  achieving 
the  potential  of  HgCdTe  because  of  Its  "structural"  weakness  which  make  It 
difficult  to  process  and  passivate  devices  or  device  arrays.  A  key  contribu¬ 
tion  from  this  contract  was  establishing  of  a  one  to  one  relationship  between 
the  bandgap  tunablllty  and  the  structural  weakness  In  HgCdTe.  As  was  estab¬ 
lished  under  the  contract  both  properties  are  due  to  the  peculiar  nature 
states  of  Hg.  In  order  to  dully  exploit  the  results  of  this  contract.  It  Is 
critical  that  the  key  role  of  Hg  be  understood.  Doing  so,  several  avenues  of 
work  are  clearly  called  for.  This  will  be  discussed  below. 

One  avenue  Is  to  explore  materials  which  have  the  bandgap  "tunablllty"  of 
HgCdTe  but  which  may  have  strong  lattices.  Since  the  band  tunablllty  Is 
uniquely  due  to  Hg,  this  must  be  retained  In  the  lattice.  An  Important  result 
from  this  contract  was  the  finding  that  should  have  more  mechanical 

strength  than  HgCdTe  while  retaining  the  Eg  tunablllty.  It  Is  very  Important 
that  this  be  followed  up  by  work  on  ^  contract  Is  being  planned 

by  NASA  and  the  Army  to  get  such  work  underway.  However,  much  more  support  Is 
necessary  to  fully  exploit  this  opportunity  which  was  revealed  by  the  present 
contract.  One  area  of  Importance  not  covered  by  the  NASA-Army  contract  Is 
work  on  surface  passivation  and  metallization. 

There  have  been  Informal  suggestions  that  alloying  HgCdTe  with  up  to 
about  10%  MnTe  might  give  material  with  the  same  IR  properties  but  more 
strength  than  HgCdTe.  For  completeness.  It  would  be  Important  to  pursue  such 
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possibilities.  However,  this  seems  less  promising  than  Hgj_jjZn^Te.  It  ap¬ 
pears  increasingly  clear  that  exploration  and  exploitation  of  CdTe-HgTe  super¬ 
lattices  will  be  a  very  important  future  IR  project;  however,  ZnTe-HgTe  super¬ 
lattices  should  not  be  overlooked  from  the  reasons  given  above.  Here  ques¬ 
tions  of  the  stability  of  the  HgTerCdTe  (HgTe:ZnTe)  interfaces  are  very  impor¬ 
tant.  The  present  contract  gives  a  good  background  into  these.  The  experi¬ 
mental  and  theoretical  approaches  established  here  should  be  exploited  for 
future  work  on  these  critical  Interfaces.  This  contract  has  established  that 
because  of  the  Hg  "Instability"  in  HgCdTe,  Interfaces  on  this  or  other  Hg 
containing  materials  (such  as  Hgj_j^Zn^Te)  provide  very  unique  problems.  Three 
types  of  Interfaces  are  Important. 

•  semiconductor : Insula tor 

•  semiconductor :metal 

•  semiconductor : semiconductor 

The  superlattices  referred  to  in  the  preceding  paragraph  are  one  example 
of  the  semiconductor : semiconductor  interface.  Much  effort  has  been  spent  to 
date  in  studying  a  very  limited  number  of  these  interfaces.  However,  this 
past  work  has  principally  been  empirical  where  a  very  specific  problem  has 
been  attached  in  the  context  of  a  given  application.  It  would  be  much  more 
efficient  to  attack  these  problems  as  well  in  a  very  generic  way.  To  do  this, 
a  number  of  avenues  opened  by  the  present  contract  should  be  exploited. 

Within  this  context  we  will  now  make  more  specific  suggestions. 

The  Hg  Instability  is  a  key  problem  with  HgCdTe  (or  other  long  wave¬ 
lengths,  Hg  containing  IR  detectors).  We  have  shown  in  this  contract  that 
this  is  particularly  critical  at  the  Interfaces  listed  above.  We  have  also 
established  that  Hg  loss  can  be  produced  at  surfaces  in  a  very  controlled 
manner  by  an  electron  beam.  Studies  of  the  effects  of  such  loss  should  be 
done  on  a  wide  range  of  parameters  Important  to  device  operation. 

This  contract  has  studied  in  some  detail  the  Cd  loss  produced  on  standard 
(bromine-methanol)  preparation  of  HgCdTe  surfaces.  The  Hg  loss  would  be 
expected  to  effect  properties  in  a  similar  way  to  this  loss;  however,  detailed 
studies  of  this  loss,  its  cause  and  effects,  are  also  of  importance. 

The  disturbances  and  loss  of  the  Hg  In  HgCdTe  lattice  caused  by  vapor 
deposition  of  metals  has  been  established  and  studies  in  this  report.  The  IR 
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detector  technology  has  advanced  to  the  point  where  such  phenomena  may  effect 
devise  operation.  Thus,  these  phenomena  should  be  studied  In  a  systematic 
fashion. 

It  has  also  been  established  that  oxidation  disturbs  the  HgCdTe  sur¬ 
face.  This  has  strong  Implications  Into  passivation;  as  does  the  realization 
that  deposition  of  a  foreign  oxide  can  disturb  the  surface.  It  Is  critical 
that  this  be  studied  In  a  systematic  way  In  order  to  reduce  the  ultimate  costs 
In  developing  IR  detector  arrays. 

These  are  but  a  few  of  the  major  technological  problems  presently  facing 
II-VI  semiconductor  IR  technology.  The  successful  exploitation  of  these 
exotic  II-VI  semiconductor  alloys  Is  dependent  on  many  of  these  problems  being 
solved.  The  fundamental  understanding  Into  the  bonding  of  HgCdTe  In  this 
present  contract  has  provided  an  Impetus  to  study  these  new  ternary  alloy 
systems,  and  the  future  successful  manufacture  of  devices  made  from  these 
semiconductors  Is  dependent  on  the  Increased  understanding  Into  the  peculiar 
nature  of  the  bonding  in  these  lattices. 
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3.1  ROLE  OF  Hg  BONDING  IN  METAL /Hg, _^Cd^Te  INTERFACE  FORMATION 

XX  X 

Abstract 

We  discuss  the  role  of  the  weak  bonding  of  Rg  In  the  lattice 

upon  Interface  morphology  when  metal  Is  deposited  on  the  surface.  The  Ag/ 
Hgj_jjCdjjTe ,  Cu/Hgj_^Cd^Te ,  Al/Hg^.j^Cd^Te ,  and  Cu/CdTe  Interfaces  are  studied 
In  ultra-high  vacuum  with  X-ray  and  ultraviolet  photoemlsslon  and  with  LEED. 
The  three  metals  are  found  to  yield  distinctly  different  Interface  morpholo¬ 
gies.  For  Ag/Hgj_^Cd^Te ,  the  Hg  instability  results  in  vacancy  diffusion  of 
2  3 

Ag  10  -10  A  into  the  semiconductor.  For  Ca  overlayers,  the  instability 
greatly  enhances  Te  and  Cd  Intermixing  into  the  Cu  for  ^Rj-x^x^®  compared  to 
CdTe.  The  more  reactive  A1  overlayers  deplete  the  weakly  bound  Hg  from  the 
Hg2_^CdjjTe  surface.  We  compare  the  metal/Hgj_jjCd^Te  Interfaces  with  the 
corresponding  metal/CdTe  and  metal/GaAs  Interfaces.  For  the  three  overlayer 
metals  on  Hg2_^Cd^Te  the  instability  of  the  Hg  bond  results  in  Interfaces  with 
a  much  greater  degree  of  intermixing  than  for  the  corresponding  metal/binary 
semiconductor  Interfaces. 

Introduction 

Hgi-j^CdxTe  is  a  semiconductor  alloy  with  a  bandgap  spanning  the  range  0 

to  1.5  eV  as  the  composition  is  varied  from  x  ■  0  to  x  ■  1.  As  such  it  is  an 

Important  candidate  for  IR  photodetection.  Intimately  related  to  this  bandgap 

tunabillty  is  a  marked  instability  in  the  HgTe  bond.^  The  weakness  of  this 

bond  is  reflected  in  the  much  lower  heat  of  formation  AH^  for  HgTe  than  for 

CdTe:  AH^(HgTe)  -  -7.6  kcal/mole,  while  AH^(CdTe)  -  -24.1  kcal/mole.^  Fur- 

3  1 

ther,  there  is  theoretical  as  well  as  experimental  evidence  that  the  HgTe 
bond,  which  is  already  weak  in  HgTe,  is  further  destabilized  in 
the  presence  of  the  Cd.  This  instability  plays  a  critical  role  in  electronic 
and  mechanical  properties  of  the  material,  both  for  the  bulk  and  for  surfaces 
and  interfaces.^ 

Both  for  applications  to  IR  photodetection  and  to  further  the  understand¬ 
ing  of  Schottky  barrier  formation  in  general,  it  is  of  interest  to  study  the 
Interface  properties  of  the  raetal/Hg2_jjCdjjTe  junction.  As  might  be  expected, 
the  weakness  of  the  HgTe  bond  plays  a  crucial  role  in  determining  the  nature 
(morphology  and  band  bending)  of  the  interface.  In  this  paper  we  explore  the 
behavior  of  the  metal/Hgj_jjCd^Te  Interface  morphology  through  photoemission 
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and  LEED  studies  of  Interfaces  of  Hg]^.](Cd^Te  with  three  metals  of  Increasing 
reactivity:  Ag,  Cu  and  Al;  the  Cu/CdTe  Interface  was  also  studied.  For  all 
three  metals,  the  HgTe  bond  Instability  Is  found  to  result  In  highly  nonabrupt 
Interfaces,  with  marked  depletion  of  Hg  from  the  semiconductor  near-surface 
for  the  more-reactive  overlayers,  and  Intermixing  of  the  semiconductor  and 
overlayer  for  the  less-reactive  cases. 

Table  3-1  shows  the  heats  of  formation  of  Te  of  the  noble  metals  and 

of  Al.  Also  shown  are  the  heats  of  alloying  AHg^2(C;M)  of  the  cations  (C)  Cd 

and  Eg  at  Infinite  dilution  In  the  overlayer  metals  (M)  ,  calculated  using  the 

4  5 

semlemplrlcal  model  of  Mledema.  *  Based  on  their  small  heats  of  formation 
with  Te,  the  noble  metals  would  be  expected  to  form  nonreactlve  Interfaces 
with  Hgj_^Cd^Te,  while  Al,  with  Its  large  AH£(Al2Te2),  should  In  contrast  be 
strongly  reactive  with  Hgj_^Cd^Te.  In  accordance  with  these  considerations, 
there  Is  clear  evidence  of  strong  overlayer-Te  formation  for  Al  overlayers  but 
not  for  Au,^  Ag,  or  Cu  overlayers.  However,  although  the  noble  metals  have 
comparable  heats  of  Te  formation,  Ag/Hgj_jjCd^Te  and  Cu/Hg^_^Cd^Te  prove  to 
have  dissimilar  morphologies,  while  Au/Hgj_jjCd^Te ,  for  which  we  refer  to  the 
results  of  Davis  et  al.,^  presents  a  porphology  distinct  from  both  the  Ag  and 
Cu  Interfaces.  An  understanding  of  these  differences  requires  consideration 
of  the  heats  of  solution  as  well  as  AH^.^  Finally,  the  large  heat  of 

Te  formation  for  Al  gives  the  Interface  of  Al  with  ®  morphology 

which  differs  from  Che  Interfaces  with  the  noble  metals. 


Table  3-1.  Heats  of  Te  Formation  AH^,  and  Heats  of  Solution  AHgQj^(C;M) 
of  Cation  C  in  Metal  M  in  the  Solid  Phase,  for  Various  Metals 


Metal  (M) 

Tellur Id e  (Te) 

AHj* 

(kcal/mole) 

AH  o^(Cd;M)** 
(kcal/mole) 

AH  i(Hg;M)** 
(kcal/mole) 

Cd 

CdTe 

-24.1 

Hg 

HgTe 

-7.6 

Ag 

Ag2Te 

-8.6 

-3.4 

-2.1 

Oi 

Ou^Te 

-10.0 

-1.7 

+0.3 

Au 

AuTe2 

-4.5 

-12.2 

-6.4 

Al 

-76.2 

+3.4 

+4.1 

*From  Ref.  2. 

**Calculated  from  Mledema' s  semlemplrlcal  model  (Refs.  4  and  5). 
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Experimental 

Single-crystal  bars  of  p-type  l^j_jjCd^Te  (x  ■  0.23  for  Ag  and  A1  over¬ 
layers,  X  ■  0.25  for  Cu)  and  n-CdTe  (In-doped  to  1  x  10^^/cm^)  with  cross- 
sectional  areas  of  5-5  mn  were  transferred  into  a  previously-baked  vacuum 
chamber  and  then  cleaved  in  ultra-high  vacuum  (base  pressure  6  x  10”^^  torr) 
to  reveal  a  [110]  face.  For  the  cleaved  surfaces  the  Fermi  level 

Ep  was  located  near  the  conduction  band  minimum  (CBM)  making  the  surfaces 

p 

n-type,  as  is  often  observed.  For  CdTe,  Ep  for  the  clean  cleaved  surface  was 

located  0.3  eV  below  the  CBM.  (The  bulk  bandgap  is  0.20  eV  for  x  «  0.23, 

0.22  eV  for  x  -  0.25,  and  1.49  eV  for  CdTe.)  Sequential  depositions  of  metal 
were  evaporated  from  a  tungsten  filament  onto  the  room  temperature  surface, 
which  was  studied  at  each  stage  by  X-ray  and  ultraviolet  photoelectron  spec¬ 
troscopy  (XPS  and  UPS)  using  respectively  a  Mg  Ka  X-ray  source  (hv  ■ 

1253.6  eV)  and  a  helium  lamp  (He  I  and  He  II,  hv  -  21.2  and  40.8  eV).  The 
photoelectrons  were  analyzed  with  a  double-pass  cylindrical-mirror  analyzer 
(CMA).  All  Interfaces  except  for  Al/Hgj_^Cdj^Te  were  also  characterized  by 
LEED  using  electron  energies  less  than  65  eV  and  beam  currents  less  than 
1  uA.  Vfe  observed  no  electron- beam- Induced  surface  disruption  for  Hg,  Cd  Te 
or  CdTe,  in  contrast  to  the  disruption  reported  for  CdTe  in  Ref.  9.  The 
amount  of  metal  deposited  is  given  in  units  of  a  monolayer,  ML,  which  we 
define  to  be  the  surface  density  of  atoms  on  the  Hg2_jjCdj^Te  [110]  face:  for  x 
■  0.23,  1  ML  ■  6.76  X  10^^  atoms/cm^  which  corresponds  to  1.15A  of  metallic 
Ag,  0. 80A  of  Cu,  and  1.12A  of  Al.  The  variation  of  surface  atomic  density 
with  x-value  over  the  entire  composition  range  is  IZ,  which  is  negligible  for 
our  purposes . 


Results 

In  subsections  A  to  C  below,  the  three  metal/Hg^_^Cd^Te  interfaces  are 
compared  with  the  respective  metal/CdTe  and  metal/GaAs  Interfaces  with  em¬ 
phasis  on  the  role  that  the  Hg  bonding  plays  for  Hgj_^Cd^Te.  Full  discussions 
of  the  Ag/Hgj_^CdjjTe^®’ and  Al/Hg^_j^CdjjTe^  ^  interfaces  have  been  pub¬ 
lished.  The  latter  interface  has  been  investigated  previously  by  Daniels  et 
12  13 

al.  and  Davis  et  al.  We  present  here  a  brief  preliminary  discussion  of 

Cu/Hgj_jjCd  Te  and  Cu/CdTe;  a  more  complete  analysis  will  be  presented  else- 
14 

where . 
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Before  presenting  the  results,  we  give  a  brief  overview.  The  differences 
In  Interface  morphology  mentioned  In  the  Introduction  are  reflected  In  the 
He  I  valence  band  spectra  in  Figure  3-1,  which  shows  spectra  for  the  clean 
Hg^_jjCdjjTe  surface  as  well  as  for  the  surface  covered  by  5  ML  of  each  over¬ 
layer  metal.  The  Ag4d  intensity  for  5  ML  Ag/Hg^.j^Cd^Te  is  much  less  than  the 
Cu3d  intensity  for  an  equal  coverage  of  Cu/Hg2_jjCd^Te.  The  Al/Hg2_j^Cd^Te 
spectrum  Is  striking  for  the  almost  complete  attenuation  of  the  Hg5d  level  at 
the  relatively  low  5  ML  coverage. 

Ag  overlayers 

Because  at  hv  *  21.2  eV  the  Ag4d,  Cd4d  and  Hg5d  states  have  similar  cross 
sections, the  comparatively  low  Intensity  of  the  Ag4d  In  Figure  3-1  Implies 
that  either  the  Ag  Is  not  sticking  to  the  surface  or  It  diffuses  well  Into  the 
semiconductor.  A  detailed  analysls^*^  of  the  XPS  core  levels  shows  that  the  Ag 
diffuses  10  -1-  A  Into  the  bulk  of  the  semiconductor.  Only  after  on  the  order 
of  100  ML  Ag  has  been  deposited  does  the  surface  attain  a  metallic  character, 
with  the  valence  band  spectrum  showing  an  occupied  density  of  states  up  to  the 
Fermi  level.  While  the  Te  and  Cd  concentrations  in  the  near-surface  remain 
essentially  constant  up  to  100  ML  Ag  deposition,  the  Hg  concentration  de¬ 
creases  by  3Z  with  each  ML  of  Ag  evaporated  onto  the  surface.^®  It  Is  thus 
reasonable  to  Infer  that  the  Ag  moves  into  the  lattice  by  displacing  and 
then  moving  Into  the  resulting  vacancies,  permitting  the  high  solubility 
associated  with  vacancy  diffusion.  After  Refs.  10  and  11  were  published  we 
repeated  the  study  of  Ag/Hg^_^Cd^Te ,  this  time  including  characterization  of 
the  surface  order  with  I£ED.  In  support  of  the  Hg-replacement  model  of  Ag 
diffusion,  a  LEED  pattern  showing  the  original  1x1  surface  order  Is  observed 
at  coverages  as  high  as  40  ML.  At  100  ML  the  LEED  pattern  is  gone.  Indicating 
that  the  surface  is  now  covered  with  polycrystalline  metallic  Ag  or  that  the 
concentration  of  Ag  in  the  Hgj_^Cd^Te  surface  at  this  coverage  is  high  enough 
to  disrupt  the  lattice  order  of  the  surface.  No  evidence  is  seen  for  signifi¬ 
cant  chemical  reaction  between  the  overlayer  and  the  semiconductor  (although 
the  data  do  not  completely  rule  out  the  occurrence  of  such  reaction). 

The  Importance  of  the  presence  of  easily  removed  Hg  in  formation  of  the 
Ag/Hgj_^Cd^Te  interface  Is  verified  by  a  comparison  with  Ag/CdTe  from 
Ref.  9.  With  no  weakly  bound  cations  which  Ag  can  replace  In  the  lattice,  Ag 
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Figure  3-1.  Valence  Band  Spectra  at  hv  *  21.2  eV  of  the  Clean  Hg,  Cd  Te 
Surface  as  Well  as  for  the  Surface  Covered  by  5  ML  of  Each 
Overlayer  Metal.  The  Cu/Hgj_j^Cd^Te  spectrum  is  for  the 
composition  x  ■  0.25;  the  other  three  spectra  are  for  x  ■ 

0.23.  A  spectrum  of  clean  x  ■  0.25  Hg^  Cd^Te  is  not  shown  as 
it  is  very  similar  to  the  clean  x  »  0.23  spectrum. 

should  not  diffuse  into  CdTe.  In  agreement  with  this  reasoning,  no  movement 
of  Ag  into  the  semiconductor  is  observed  for  Ag/CdTe;  this  interface  is  es- 
sentlally  abrupt,  in  dramatic  contrast  to  ^g/Hgj _j^Cd^Te  (a  small  Cd  movement 
into  the  overlayer  is  observed  in  Ref.  9,  but  this  may  well  be  due  to  the 
anomalous  instability  of  the  surface  observed  during  their  LEED  work).  The 
behavior  of  Ag  on  CdTe  closely  parallels  its  behavior  on  the  III-V's.  Ag/ 
InP^^  and  Ag/GaAs,^^  for  Instance,  have  been  shown  to  be  nonintermixing,  with 
the  only  Important  deviation  from  the  "ideal"  abrupt  interface  being  Ag  island 
I  formation  Instead  of  pure  laminar  growth. 

Cu  overlayers 

I  Because  of  the  similarities  between  Cu  and  Ag  in  m’tal-Te  heats  of  forma- 

I  tion,  the  two  metals  might  naively  be  expected  to  form  similar  interfaces  on 

Hg,  Cd  Te.  However,  as  was  mentioned  in  the  introduction,  heats  of  alloying 

1  X 
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will  also  play  a  role  In  determining  Interface  reactivity  (this  point  Is 
discussed  In  greater  detail  In  Section  4).  Thus  the  valence  band  spectra  for 
5  ML  of  the  two  metals  on  Hgi_„Cd^Te  In  Figure  3-1  differ  In  two  ways:  The 
Hg5d  Intensity  decreases  from  the  clean  spectrum  to  the  5  ML  coverage  much 
more  for  Cu  overlayers  than  for  and  while  the  Ag4d  Is  barely  visible  due 
to  diffusion  of  Ag  Into  the  bulk,  the  Cu3d  Intensity  Is  comparable  to  that  of 
the  Cd4d  and  I^5d  levels.  Since  at  hv  >  21.2  eV  all  three  subshells  have 
similar  photolonl ration  cross  sections, there  can  be  no  deep  diffusion  of  Cu 
Into  the  semnlconductor  as  there  Is  for  Ag. 

More  details  of  the  Interface  morphology  can  be  obtained  from  the  XPS 
core  levels.  Figure  3-2(a)  shows  as  a  function  of  Cu  coverage  the  peak  areas 
I  of  the  semiconductor,  normalized  to  the  zero-coverage  values  Iq.  If  the 
Interface  for  Cu  were  abrupt  and  the  metal  film  laminar,  1/Iq  would  fall 
exponentially  with  coverage,  with  a  1/e  length  of  order  20  ML.  Only  for  the 
Cd3d  peak  below  10  ML  is  such  attenuation  observed.  Above  10  ML  the  attenua¬ 
tion  of  the  Cd3d  peak  area  slows  so  that  at  100  ML  coverage  I/Iq  40%.  The 
deviation  from  exponential  attenuation  Is  even  more  marked  for  the  Te3d  peak 
area,  which  actually  rises  5%  by  the  10  ML  coverage,  then  levels  off  and  falls 
slowly  to  I/Iq  ■  85%  at  100  ML  coverage;  there  is  virtually  no  further  attenu¬ 
ation  during  the  next  200  ML  deposition.  The  attenuation  of  the  core  intensi¬ 
ties  Is  much  too  slow  to  be  explainable  by  formation  of  Cu  Islands  at  the 
Interface.  The  Interface  Is  thus  highly  nonabrupt ,  with  significant  intermix¬ 
ing  of  the  surface  up  to  coverages  of  several  hundred  ML.  The  Hg  behaves  very 
differently  from  the  Cd  or  Te ,  with  tte  Hg4f  peak  area  decreasing  to  I/Iq  = 

80%  by  0.6  ML,  and  then  attenuating  exponentially  until  10  ML  with  a  1/e 
length  of  9  ML  (llA).  The  escape  length  In  Cu  (or  a  Cu-Cd-Te  alloy)  for  elec¬ 
trons  photoemltted  with  1150  eV  kinetic  energy  from  the  Hg4f  subshell  should 
be  about  20A.  Thus  the  Hg4f  attenuation,  especially  below  0.6  ML,  is  too 
fast  to  be  an  escape  depth  effect  and  must  be  due  to  depletion  of  Hg  from  the 
surface. 

Some  intermixing  Is  observed  for  the  Cu/CdTe  Interface  as  well,  but  to  a 
much  smaller  degree  than  for  Cu/Hgj_jjCdjjTe.  Figure  3-2(b)  shows  the  peak 
areas  of  the  CdTe  core  levels  as  a  function  of  Cu  coverage,  normalized  as  In 
Figure  3-2(a).  Both  Te  and  Cd  core  levels  attenuate  faster  than  for  Cu/ 
Hgj_jjCd^Te.  For  Te3d,  I/Iq  ■  20%  at  100  ML,  with  little  attenuation 
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Figure  3-2(a).  Attenuation  of  the  XPS  Peak  Areas  I  of  the  Hgj_^Cd^Te  Core 
Levels  With  Increasing  Cu  Coverage.  The  peak  areas  are 
normalized  to  their  zero-coverage  values  Iq. 

Figure  3-2(b).  Attenuation  of  the  Semiconductor  Core  Level  Peak  Areas  as  in 
(a),  But  for  Cu/CdTe. 

thereafter.  The  Cd3d  attenuates  exponentially  with  a  15  to  20  A  (19  to  25  ML) 
escape  length  at  all  coverages,  with  no  signal  observable  beyond  100  ML. 

Thus,  little  or  no  Cd  Is  Intermixed  from  CdTe  Into  the  overlayer,  while  the 
concentration  of  Te  Intermixed  Into  the  overlayer  for  CdTe  Is  about  1 /A  of  the 
Intermixed  Te  concentration  for  Hgj.j^Cd^Te.  Formation  of  islands  of  Cu  on  the 
surface  would  also  slow  the  attenuation  of  the  Te3d  level  below  that  expected 
for  a  laminar  interface,  but  then  the  Cd  and  Te  intensities  would  be  expected 
to  attenuate  at  similar  rates.  Thus,  while  Is  landing  of  the  Cu  Is  not  un¬ 
likely,  the  dependency  of  I/Iq  on  the  atom  being  examined  unambiguously  Indi¬ 
cates  Intermixing  of  Te  Into  the  overlayer.  The  depletion  of  the  weakly 
bonded  in  l^j_j^Cdj^Te  provides  an  excess  of  Te  as  well  as  presumably  col¬ 
lapsing  the  lattice,  thus  causing  substantially  greater  intermixing  for 
Hgj_^Cd^Te  than  for  CdTe. 
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The  similarities  between  metal/CdTe  and  metal/GaAs  are  even  greater  for 
Cu  than  for  Ag.  For  Cu/GaAs,  both  Ga  and  As  Intermix  with  the  overlayer, 
with  the  anion  being  the  predominant  intermixed  species  seen  in  the  overlayer 
at  high  coverages.  Indeed,  the  As3d  peak  area  at  hv  >  85  eV  attenuates  suf¬ 
ficiently  slowly  that  at  100  ML  I/Iq  *  30%,  showing  that  the  anion  concentra¬ 
tion  in  the  Cu  overlayer  is  comparable  for  GaAs  and  CdTe.  The  concentration 
of  Ga  in  the  Cu  overlayer  is  much  smaller  than  the  concentration  of  As,  with 
1/Iq  -  2%  at  100  ML. 

A1  overlayers 

The  He  I  spectrum  for  Al/Hgj_jjCd^Te  in  Figure  3-1  shows  a  very  rapid  Hg5d 
signal  attenuation,  with  I/Iq  5%  at  5  ML,  implying  a  significant  depletion 
of  Hg  from  the  surface.  The  attenuation  of  the  semiconductor  core  level 
intensities  with  A1  coverage  is  shown  in  Figure  3-3.  The  Cd3d  attenuates  at 
a  rate  consistent  with  attenuation  through  an  abrupt  A1  overlayer,  while  the 
Te3d  signal  attenuates  much  more  slowly,  with  I/Iq  =  60%  at  40  ML.  This  slow 
attenuation  of  the  Te3d  intensity  is  consistent  with  the  Te-rich  surface- 
segregated  layer  reported  by  Davis  et  al.  Verifying  the  expected  reactivity 
of  this  interface,  at  coverages  below  2  ML  there  is  one  A12p  peak  correspond¬ 
ing  to  reacted  Al  (i.e.  shifted  to  higher  binding  energy  than  for  the  metallic 
A12p)  ,  while  at  higher  coverages  a  metallic  Al  contribution  arises  at  0. 9  ± 

0.2  eV  lower  binding  energy  and  becomes  dominant  as  the  interfacial  reacted 
layer  is  covered  up  by  unreacted  Al .  A  detailed  analysis  of  the  interface  is 
given  in  Ref.  11.  It  shows  that  as  the  first  2  to  3  ML  of  Al  are  deposited, 

Hg  is  depleted  to  50%  or  less  of  its  original  concentration  to  a  depth  of  10 
to  30A  into  the  semiconductor;  there  is  little  or  no  further  Hg  depletion  at 
higher  coverages.  The  Hg-depleted  semiconductor  surface  is  covered  by  a  3  to 
5  ML  thick  AlTe  reacted  layer,  which  in  turn  is  covered  by  metallic  Al.  A 
several  ML  thick  Te-rich  reacted  layer  is  segregated  to  the  surface  of  the 
metallic  Al  layer.  A  similar  "'-'terface  morphology  is  seen  for  the  reactive 
In^^  and  overlayers. 

For  Al/Hg^  _jjCd^Te ,  a  relatively  sharp  peak  originating  from  the  Te  p-like 

states  of  the  Te-rich  layer  floating  on  the  metallic  Al  evolves  in  the  UPS 

1  1  22 

valence  band  at  high  coverages  3.2  eV  below  Ep.  Patterson  and  Williams 
show  a  very  similar  peak  evolving  at  high  coverages  for  Al/CdTe,  which 
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Figure  3-3.  Attenuation  of  the  XPS  Peak  Areas  I  of  the  Hgj_^Cd  Te  Core 
Levels  With  Increasing  A1  Coverage.  The  peak  areas  are 
normalized  to  their  zero-coverage  values  Iq. 

suggests  the  presence  of  outdiffused  Te  in  or  on  the  A1  for  this  interface  as 
well.  They  also  show  evidence  for  the  presence  of  metallic  Cd  dissociated 
from  the  CdTe,  as  manifested  by  the  growth  of  a  Cd4d  peak  at  0.6  eV  lower 
binding  energy  than  for  Cd  bound  to  Te.  (Note  that  the  formation  of  metallic 
Cd  does  not  imply  Cd  intermixing  with  the  overlayer.  For  Al/CdS,  for  ex¬ 
ample,  the  dissociated  metallic  Cd  is  localized  at  the  semiconductor  sur¬ 
face.)  No  metallic  Cd  is  seen  for  Al/Hgj_^Cd^Te ,  suggesting  that  the  presence 
of  is  inhibiting  formation  of  dissociated  metallic  Cd  for  the  Hgj^^Cd^^Te 
Interface.  One  contributing  factor  may  be  the  enhancement  of  the  CdTe  bond 
strength  in  Hgj.j^Cd^^Te  compared  to  CdTe.^  Probably  mor<'  important  is  the 
nonstolchlometry  of  the  first  10  to  30A  of  the  semiconductor  due  to  Hg  deple¬ 
tion  at  low  A1  coverages.  For  x  -  0.23  a  depletion  of  half  the  Hg  at  the 
surface  creates  a  large  excess  of  Te,  which  is  then  available  to  help  trap  any 
Cd  which  might  otherwise  (as  in  Al/CdTe)  have  moved  into  the  overlayer.  This 
Te  excess  also  provides  Te  for  the  AlTe  interfacial  reacted  layer  and  for  the 
Te-rlch  segregated  surface  layer. 
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For  Al/GaAs,^^  a  limited  metal-cation  exchange  reaction  occurs,  with  a 

small  amount  of  released  Ga  present  in  the  overlayer  near  the  semiconductor 

surface.  The  reaction  appears  to  be  confined  to  within  at  most  a  few  crystal 

25 

layers,  reflecting  the  large  reactivity  of  the  interface.  This  decrease  in 

interface  width  with  increasing  reactivity  for  GaAs  appeal’s  to  hold  for 

Hgi_xCd^Te  as  well  if  the  interface  width  is  considered  to  be  the  width  over 

which  the  composition  changes  from  semiconductor  to  bulk  metal.  However, 

(and  probably  CdTe)  have  a  surface  segregated  anion  layer  which 

2  6 

does  not  appear  in  GaAs.  It  has  been  observed  more  generally  that  the  ratio 
of  outdlf fused  anion:  outdlf fused  cation  decreases  with  increasing  overlayer 
reactivity  for  III-V's  but  not  for  II-VI’s. 

Discussion 

For  the  three  overlayer  metals  discussed  in  this  paper,  the  trends  in 
interface  morphology  for  CdTe  parallel  those  for  GaAs.  Ag  forms  a  noninter- 
mlxed  interface  on  both  semiconductors,  while  for  Cu  there  is  a  significant 
concentration  of  the  anion  in  the  overlayer  up  to  several  hundred  ML  cover¬ 
age.  This  similarity  between  the  two  semiconductors  for  the  noble  metals  also 
holds  for  Au.  Au/GaAs  is  at  least  as  intermixing  as  Cu/GaAs:  probed  with  hv 

•  100  eV  (Giving  high  surface  sensitivity),  I/Iq  ■  ^0  to  50%  for  As3d  and  ■ 

28 

10%  for  Ga3d  after  80  ML  coverage.  In  comparison,  for  Au/CdTe  at  17  ML  cov¬ 
erage,  I/Iq  *  60%  for  the  Te4d  and  •  45%  for  the  Cd4d  (probed  with  high  sur¬ 
face  sensitivity),  making  Au/CdTe  somewhat  more  intermixed  than  Cu/CdTe.  For 
both  semiconductors  the  reactive  overlayer  A1  forms  more  localized  interfaces 
than  for  the  noble  metals,  but  the  limiting  of  the  interface  width  is  much 
more  pronounced  for  GaAs. 

With  the  introduction  of  Ife  into  the  lattice,  the  overlayers  formed  on 
Hgj.xCdx^®  noble  metals  and  A1  differ  significantly  from  the  respec¬ 

tive  interfaces  for  CdTe  and  GaAs.  In  all  cases  the  instability  of  Hg  in  the 
lattice  serves  to  delocalize  the  interface  to  widths  of  a  few  tens  of  ang¬ 
stroms  for  A1  (plus  the  surface  segregated  Te-rich  layer),  and  hundreds  of 
angstroms  for  the  noble  metals.  The  mechanism  of  the  intermixing  enhancement 
is  different  for  the  three  overlayers  studied  here,  as  well  as  for  Au.^  Fur¬ 
thermore,  the  Cd  signal  at  high  coverages  is  greater  in  Au/CdTe  than  in 
Cu/CdTe,  while  for  Hgj_jjCd^Te  the  situation  is  reversed,  with  a  greater  Cd 
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signal  in  Cu/Hg,  Cd  Te  than  in  Au/Hg,  Cd  Te./  These  differences  can  be 

A  X  X  XXX  ‘ 

understood  qualitatively  through  consideration  of  the  heats  of  alloying  of  the 
semiconductor  cations  in  the  overlayer  metal.  The  heats  of  Te  formation  and 
cation-overlayer  are  shown  in  Table  3-1. 

The  relative  abruptness  of  the  Ag/CdTe  interface  is  due  to  a  lack  of 
chemical  reaction  because  of  silver's  low  heat  of  Te  formation  AH£(Ag2Te)  and 
allo)d.ng  AHgQj^(Cd;Ag)  .  For  Hg^_^Cd^Te,  the  weakly  bound  Hg  is  driven  out  by 
the  heat  of  condensation  of  the  Ag  on  the  surface.  This  leads  to 

formation  of  Hg  vacancies  through  which  the  Ag  can  diffuse,  driven  by  the  con¬ 
centration  gradient  in  the  chemical  potential  across  the  Interface.  Te  forma¬ 
tion  and  alloying  can  provide  a  barrier  to  this  diffusion,  but  because 

AH^(Ag2Te)  and  AH  .(Cd;Ag)  are  small  the  barrier  is  low,  permitting  the  ob- 

7  ^ 

served  motion  of  Ag  10  -10  A  into  the  Hgi_„Cd  Te. 

X  X 

Au  also  has  a  small  heat  of  Te  formation,  but  its  heat  of  alloying  is 
much  larger:  AH^^j^C Cd ; Au)  ■  -14  kcal/raole,  which  is  four  times  the  value  of 
AHsoi( Cd ;Ag) .  The  large  AH^^j^C Cd ;Au)  is  sufficient  to  disrupt  the  CdTe  sur¬ 
face,  pulling  Cd  into  the  overlayer.  The  presence  of  Te  in  the  over layer  is 
probably  also  favored  by  alloying  considerations,  but  this  cannot  be  verified 
without  information  on  AH^^j^CTe  ;Au) .  In  contrast,  for  Cu/CdTe  little  or  no  Cd 
is  observed  in  the  overlayer  due  to  the  small  heat  of  alloying  for  Cu: 
AHsoi(Cd;Cu)  ■  2  kcal/mole.  The  slightly  greater  heat  of  Te  formation  for  Cu 
than  for  Ag  apparently  is  sufficient  to  pxill  Te  into  the  overlayer  for  Cu/CdTe 
while  leaving  Ag/CdTe  abrupt. 

In  the  case  of  Hgj_jjCd^Te  interfaces  with  Cu  and  Au ,  the  Hg  instability 
alters  the  intermixing  compared  to  CdTe.  For  Cu/Hgj_^Cd^Te ,  formation  of  Hg 
vacancies  occurs  to  a  greater  degree  than  for  the  less  reactive  Ag/ 

Hgj_^Cd^Te.  The  resulting  surface  lattice  collapse  probably  enables  both  Cd 
and  Te  to  intermix  with  the  Cu  overlayer,  and  may  also  serve  as  a  barrier  to 
Cu  indlffuslon.  The  small  heat  of  CuCd  alloying  now  would  presumably  make  it 
energetically  favorable  for  the  Cd  to  move  to  the  surface  of  the  Cu  rather 
than  mixing  uniformly  with  the  Cu.  In  this  case  the  Cd  core  level  intensity 
would  be  greater  than  for  the  same  amount  of  Cd  mixed  uniformly  into  the  over¬ 
layer.  The  observed  strong  Cd  signal  intensity  for  Cu/Hgj_j^Cd^Te  is  consis¬ 
tent  with  such  Cd  segregation  to  the  surface  (although  this  evidence  is  not 
conclusive).  For  Au/Hgj_^Cd^Te ,  on  the  other  hand,  the  large  heat  of  AuCd 
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alloying  will  tend  to  mix  the  Cd  evenly  into  the  overlayer,  or  even  keep  it 
localized  near  the  Au/Hg^_j^Cd^Te  interface.  As  can  be  seen  from  Table  3-1, 
the  heat  of  AuHg  alloying  is  nonnegliglble,  so  that  some  Hg  also  might  be  seen 
in  the  overlayer,  despite  the  large  Hg  vapor  pressure.  Davis  et  al.^  show 
that  for  40  ML  Au/Hgj_^Cd^Te ,  I/Iq  “  20  to  30%  for  Cd4d  and  I/Iq  •  15%  for 
Hg5d  (both  probed  with  hv  ■  70  eV  giving  surface  sensitivity  on  the  order  of 
5  ML).  We  attribute  the  slow  attenuation  to  mixing  of  the  semiconductor 
cations  into  the  Au. 

Finally,  for  Al/Hg,  Cd  Te  a  strong  AlTe  reaction  is  observed.  Conse- 
quently  there  is  a  Hg  depletion  of  50%  or  more  in  the  first  10  to  30A  of  the 
semiconductor.  This  strong  alteration  of  lattice  stoichiometry,  which  does 
not  have  a  direct  parallel  for  the  binary  semiconductors  CdTe  and  GaAs,  must 
disrupt  the  lattice  and  provide  an  excess  of  Te.  A  3  to  5  ML  AlTe  reacted 
layer  on  top  of  the  Hg-depleted  region  is  covered  by  metallic  Al,  with  a  seg¬ 
regated  Te-rlch  layer  floating  on  the  surface.  Although  the  Hg  loss  and 
resulting  Te  excess  is  much  greater  for  Al  overlayers  than  for  Cu  overlayers, 
the  reactivity  of  Al  with  Te  tends  to  “chemically  trap"  the  Te ,  making  the 
presence  of  Te  in  or  on  the  overlayer  less  than  for  Cu.  The  positive  heat  of 
alloying  AHggj(Cd;Al)  of  Cd  in  Al  (Table  3-1),  and  the  presence  of  the  Te 
excess  at  the  interface,  keep  Cd  from  entering  the  overlayer  to  any 
significant  extent. 

Conclusions 

The  instability  of  the  Hg  bond  in  Hg,  Cd  Te  compared  to  the  Cd  bond 
plays  a  pivotal  role  in  determining  the  morphology  of  metal/ Hgj_j^Cd^Te  inter¬ 
faces.  Band  bending,  which  was  not  emphasized  in  this  paper,  also  is  influ¬ 
enced  by  the  Hg  bonding.  For  Ag/Hgj_jjCdj^Te ,  fo'  instance,  the  indif fused  Ag 
or  possibly  vacancies  themselves  dope  the  semiconductor  p-type  at  low 
coverages,  bringing  Ep  closer  to  the  valence  band  maximum. The  weak  Hg 
bonding  causes  metal/ Hgj_jjCd^Te  interfaces  to  show  behavior  qualitatively  dif¬ 
ferent  from  that  of  the  binary  II-VI  and  III-V  semiconductors  CdTe  and  GaAs. 
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3.2  CALCULATION  OF  CdTe  BAND  STRUCTURE 

As  an  aid  in  completing  the  analysis  of  the  angle-resolved  data,  calcula¬ 
tion  of  the  band  structure  of  CdTe  using  the  nonlocal  pseudopotential  of 
Chelikowsky  and  Cohen^  was  performed.  The  spectrum  of  computed  energies  and 
wave  function  spanne’  the  valence  bands  and  the  high  lying  conduction  bands  up 
to  30  eV  above  the  valence  band  maximum. 

The  results  of  this  calculation  shed  insight  into  the  role  of  the  actual 
final  states  Involved  in  the  normal  emission  from  CdTe  (and  by  analogy,  from 
Hg2_jfCdjjTe)  in  shaping  the  features  observed  in  the  spectra.  While  a  simple 
model  of  the  final  states  that  approximates  them  by  free  electron  parabolas 
successfully  accounts  for  much  of  the  observed  structure,  certain  peaks  in  the 
data  reject  Instead  the  mixing  of  the  parabolic  bands  that  takes  place  in  the 
crystal  potential.  The  calculation  sheds  light  on  this  aspect  in  a  semiquan- 
titative  way  (trends  are  reproduced  but  quantitative  comparison  is  not  possi¬ 
ble).  The  calculated  valence  and  high- lying  conduction  bands  along  the  E 
direction  are  Illustrated  in  Figure  3-2.  As  the  caption  explains,  states  with 
a  large  (220]  plane  wave  component  (these  contribute  most  strongly  to  the 
normal  emission)  are  marked,  and  experimentally  determined  values  of  energy 
and  k  for  selected  transitions  observed  in  the  normal  emission  data  from  CdTe 
are  also  shown.  While  the  data  clearly  falls  in  the  region  of  high  [220] 
density,  the  points  occur  at  energies  consistently  higher  than  the  calcula¬ 
tion.  The  effect  of  mixing  of  the  [220]  and  [002]  parabolas  is  especially 
strong  in  the  region  around  k  =  0.7  (2Tfa)  ,  and  it  is  here  that  the  structures 
marked  No.  1  and  No.  3  exhibit  a  nonmonotonic  dependence  of  energy  in  the 
spectra  of  k  in  Figure  3-4. 

The  calculation  was  also  used  in  a  second  way.  The  wave  functions  of 
each  band  were  used  to  examine  the  optical  transition  matrix  elements  for 
Interband  transitions.  These  are  the  central  elements  determining  the  intens¬ 
ity  of  peaks  in  the  spectra  and  contain  the  polarization  dependence  of  the 
photoexcitation.  In  all  cases  that  were  examined,  the  calculated  polarization 
dependence  and  the  measured  behavior  agreed.  This  analysis  requires  knowledge 
of  which  bands  participate  in  the  interband  transition;  the  agreement  obtained 
is  therefore  confirmation  of  the  accuracy  of  the  assignment  of  the  k-space 
origin  of  the  features  in  the  photoemission  data. 


L 


■.■’o -.'V 


RPT 60411 


HUGHES 


SANTA  SANSANA  NCSiAACH  Cf  NTIN 


>  20 
a> 


Figure  3-4. 


k  along  (110)  (2Tr/a) 


Calculated  CdTe  Conduction  Bands  of  Symmetry  and  Identifica¬ 
tion  of  States  With  Greater  Than  30^  (  )  or  Between  %  and  30% 

(a)  [220]  Plane  Wave  Character.  Also  shown  are  portions  of  the 
experimental  final  states  for  peaks  No.  1,  No.  3,  and  No.  7  (0) . 
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3.3  ANGLE -RESOLVED  PHOTOEMISSION  STUDY  OF  THE  CdTe,  HgTe  and  Hr,  ^Cd  Te  BAND 

1  “X  X 

DISPERSIONS  ALONG  MAJOR  SYMMETRY  DIRECTIONS,  AND  THE  ALLOYING  EFFECT  IN 
Hgi.^Cd^Te 


Angle-resolved  photoenlssion  spectroscopy  (ARPES)  study  on  the  band 

structure  of  HSi-^Cd^Te  along  [110]  direction  i^lng  normal  emission  method  has 

been  published  by  Sllberman  et  al.^  This  work  Is  an  extension  of  the  band 

mapping  along  other  major  symmetry  directions  ([111]  and  [100])  using  off- 

normal  emission  method.  By  assuming  direct  transition  model  (DTM)  and  the 

free  electron  final  state,  the  band  mapping  along  [100]  direction  shows  very 

2 

good  agreement  with  the  band  calculations  by  Chelikowsky  et  al.  for  CdTe 

•j 

using  nonlocal  pseudo potential  method,  and  Hass  et  al.  for  HgTe  and 
Hgj_jjCd^Te  using  empirical  tight  binding  method.  Part  of  our  mapping  results 
for  the  lower  bands  along  [100]  for  CdTe  and  HgTe  are  shown  in  Figure  3-5. 

The  assumption  of  free  electron  final  state  and  DTM  in  the  band  mapping  along 
[111]  and  mapping.  Realistic  final  state  calculation  may  be  necessary  for 
further  analysis. 


r  A  X 


Figure  3-5.  The  Band  Mappings  of  CdTe(.)  and  HgTe  (o)  Along  [100]  Direction 
Using  Off-Normal  Emission  ARPES.  The  theoretical  calculations 
for  CdTe  from  Ref.  2  (solid  line)  and  HgTe  from  Ref.  3  (dash 
line)  are  also  shown  for  comparison.  Only  the  lower  bands  are 
presented . 


■  •  '.v 
•V.'A 

C 


72 


RPT 60411 


HUGHES 


SANTA  BAAtANA  AESCAACH  Cf  NTEA 

4  '1  t'. 

Besides  the  band  mapping  on  the  CdTe,  HgTe  and  Ifei_„Cd^Te,  ARPES  is  very 

X  "X  X 

useful  in  investigating  the  alloying  effect  on  the  substitutional  alloy 
Hg^ _jjCd^Te .  In  Hg ^ _jjCdj^Te ,  the  crystal  structure  remains  single  crystal  (l.e. 
the  lattice  is  periodic);  however,  the  potential  is  aperiodic  due  to  the 
random  fluctuation  of  the  chemical  constituents.  In  treating  the  disorder 
problem  in  pseudobinary  alloys,  the  Coherent  Potential  Approximation  (CPA) 
which  uses  the  virtioal  crystal  potential  as  the  unperturbed  Hamiltonian  and 
treats  the  chemical  disorder  as  the  perturbation,  has  been  very  successful  in 
calculating  the  electronic  structures  of  the  substitutional  alloys.  In  1982 
Spicer  et  al.^  published  an  experimental  result  showing  the  breakdown  of 
virtual  crystal  approximation  (VCA)  in  they  used  the  angle- 

integrated  photoemlsslon  spectroscopy  to  show  that  the  density  of  states  as  a 
fimctlon  of  energy  agree  with  the  CPA  rather  than  VCA  calculation  result. 
Besides  the  prediction  in  the  density  of  states,  the  CPA  also  shows  that  the 
pure  Bloch  waves  in  VCA  picture  are  now  subjected  to  the  alloy  scattering 
resulting  from  the  aperiodic  potential  term.  This  will  result  in  the  broaden¬ 
ing  and  the  peak  shifting  of  the  spectral  density  function  deviating  from  the 
VCA  Bloch  state.  Because  the  disorder  in  Wg^_j^Cd^Te  manifests  itself  primar¬ 
ily  in  the  difference  of  cation  s-level  energy  associated  with  relativistic 
shifts,  only  the  self-energy  associated  with  the  cation  s-states  enters  the 
spectral  density  functions  A  (n  k,  E): 

t<  y  1* 

A(n  k,  E)  -  -ir"^  f  Z  { [E  -  E(n  k)  -  f  Z  ]^  +  (f  Z  )^}"^ 
s  s  s 

S  S  S 

where  E(n  k)  is  the  energy  of  the  VCA  eigen  state  |  n  k  >,  fg^(n  k)  represents 

the  fraction  of  time  an  electron  in  state  )  n  k  >  spends  on  cation  s  orbitals, 
»•  ? 

Zg  is  the  imaginary  part  of  self  energy,  and  Zg  is  the  real  part  of  self 

t 

energy.  As  one  can  see  from  the  expression  of  the  spectral  function,  Zg  and 
Zg  reflect  the  energy  shift  and  the  damping  of  a  particular  VCA  eigen 
state.  Using  ARPES,  we  can  restrict  ourself  into  the  optical  transition 
involving  a  particular  initial  state.  Then  by  comparing  the  ARPES  peaks  of 
HgTe,  CdTe,  and  Hgj_j^Cd^Te  of  the  same  k  value,  we  are  able  to  deduce  the 
information  of  the  peak  shifts  and  broadenings.  The  alloy  scattering  is  the 
strongest  in  the  lowerband  along  [100]  direction  because  it  contains  Che 
largest  cation  s  orbital  contribution.  In  Figure  3-6  we  show  the  experimental 
result  of  the  alloy  broadening  and  the  energy  shift  of  the  lower  band  along 
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I  [100]  direction  (the  details  of  the  experimental  deduction  will  be  published 

'  » 

'  elsewhere).  We  also  show  the  theoretical  calculation  of  the  self  energy  Z 

"3 

I  '  and  by  Hass  and  Ehrenreich  in  Figure  3-6  One  should  realize  that  the 

I  ® 

[  correct  theoretical  expression  should  also  contain  f^*^  (n  k)  which  is  not 

j  included  in  Figure  3-6.  The  theoretical  calculation  including  f^*^  (n  k)  Is 

I  underway.  The  experimental  result  shows  strong  variation  of  the  peak  shifts 

I  and  the  broadenings  of  the  peaks  versus  binding  energies  is  very  significant 

!  —  the  positions  of  the  cross  over  of  the  peak  shift  (from  negative  to  posi¬ 

tive)  and  the  maximum  broadening  agree  very  well  with  the  theoretical  predic- 

1  »• 

tion.  The  cross  over  of  Zg  and  the  maximum  of  Zg  happen  at  binding  energy 
I  ~4.8  eV  which  corresponds  the  transition  from  the  s  orbital  contribution  to 

the  Cd  s  contribution  (see  the  density  of  states  projected  on  cation  s  orbi¬ 
tals,  Figure  3-7). 
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Comparison  of  the  Theoretical  Calculation  of  (E)  for 
HgQ  iCcIq  oTe  and  the  Experimental  Result  of  the  ARPES  Peak 
Shift  From  the  VGA  Value  of  <100>  Lower  Band* 


dll^DlA/6r  (^V) 

Figure  3-6(b).  Comparison  of  2Eg  (E)  for  HgQ^yCdQ^jTe  and  the  Experimental 
Result  of  ARPES  Peak  Broadening.  The  broadening  data  was 
obtained  by  comparing  ARPES  data  of  HgTe  and  HgQ^yCdQ^jTe  of 
the  same  photon  energy  hv,  and  analyzer  angle  0y^. 
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Section  4 
ALLOY  BONDING 


Results  obtained  on  the  alloy  bonding  task  since  July  1985  are  contained 
In  the  following  five  papers.  The  first  four  are  reprints  of  papers  which 
were  published  In  this  period  while  the  last  one  Is  a  preprint  accepted  for 
publication  In  the  Journal  of  Vacuun  Science  and  Technology  A. 
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Dislocation  energies  and  hardness  of  semiconductors 
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The  dislocation  energies  and  hardness  of  semiconducton  are  calculated  by  an  extension  of 
Harrison’s  method.  It  is  demonstrated  in  agreement  with  experiment  that  dislocation  energies  per 
unit  length  are  proportional  tod  ~*-d  ~’,wheredis  the  bond  length  and  hardness  is  proportional 
tod  d  ~  The  hardness  is  related  to  the  interaction  energies  among  dislocations.  It  is  argued 
that  dislocation  densities  of  semiconductors  will  be  reduced  if  they  are  alloyed  with  a  second 
constituent  that  has  a  shorter  bond  length.  Experimental  evidence  supporting  this  strategy  is 
noted. 


Dislocations  in  semiconductors  are  detrimental  to  de¬ 
vice  function;  they  serve  as  channels  for  impurity  migration 
and  trapping,  which  cause  nonuniform  doping  and  degrades 
p-n  jimctions.'  They  also  decrease  the  material’s  resistance 
to  plastic  deformation.  The  aim  of  this  letter  is  to  provide 
insights  into  the  underlying  physical  mechanisms  control¬ 
ling  dislocations  and  semiconductor  hardness,  and  then  to 
suggest  strategies  for  decreasing  dislocation  densities.  It  is 
well  esublished  that  the  hardness  of  tetrahedrally  coordi¬ 
nated  semiconductor  materials — groups  IV,  III-V,  and  II- 
VI  compounds— exhibits  a  sharp  variation  with  their  near¬ 
neighbor  distance  d,  approximately  proportional  tod~*  for 
one  group  of  seven  compounds.^  Thus,  semiconductors  with 
small  lattice  constants  tend  to  be  harder  materials.  These 
same  materials  have  larger  stiffness  coefficients^  and  have 
fewer  dislocations  in  as-grown  crystals.^ 

The  shear  coefficients  (combinations  of  C«4.and 
^1 1  ~  ^12  in  the  Schoenflies  noution)’  depend  on  crystal 
orientation  and  (in  Harrison’s  notation*)  are  proportional  to 
Kj/d’(Kj  +  where  Kiccd'^isthecovaientand  K, 

is  the  ionic  energy.  The  metallic  interaction  modifies  the 
functional  dependence  of  the  shear  coefficient  on  K,  and 
but  introduces  no  explicit  dependence  on  the  hopping  inte¬ 
grals,  denoted  K,  by  Harrison.*  *  In  a  pure  covalent  material, 
the  bond  energy  is  proportional  to  K,  (or  d~%  and  the  bond 
volume  h  ecd\  hence,  in  this  case,  the  shear  coefficient  var¬ 
ies  as  d  In  the  limit,  K,>  K,,  C, ,  —  C,2  ee  d  ~  For  most 
polar  semiconductors,  d  is  a  good  approximation. 

Hardness  is  determined  by  applying  a  known  force  F  to 
a  probe  of  a  prescribed  shape  driving  it  into  the  surface  of  the 
sample.*  The  aread  ol  uie  resulting  indentation  is  measured, 
and  the  hardness  is  the  force  per  unit  indented  area.  Many 
disloeations  must  be  formed  to  allow  the  probe  to  indent  the 
semiconductor.  If  the  indenter  is  a  rectangular  pyramid, 
then  the  hardness  is^sF/d  ^  fh  /Ah  m  ej./AM,’wbener 
is  the  work  required  to  cause  the  indenter  to  penetrate  to  a 
depth  A.  A  side  view  of  the  indentation  in  a  cut  through  its 
center  is  illustrated  schematically  in  Fig.  1.  The  top  of  the 
indentation  has  side  length  W-,  thus,  A  »  W^.  The  Burger’s 
vector  has  magnitude  A,  proportional  to  the  bond  length  d. 
The  number  of  dislocations  A'’*  required  to  accommodate  an 


indentation  to  depth  A  is  1V„  «  A /A  IF  cos  ^ /A,  where  ^ 
is  the  angle  between  the  normal  to  the  tip  of  the  indenter  and 
aside.  Figure  1  also  shows  a  model  of  one  possible  configura¬ 
tion  of  the  dislocations.  The  edges  of  the  extra  atom  planes 
that  are  driven  from  the  indented  volume  into  the  bulk  of  the 
semiconductor  are  shown  as  lines  terminated  by  dots.  The 
dotted  ends  of  these  lines  are  the  positions  of  the  disloca¬ 
tions,  which  are  perpendicular  to  the  plane  of  the  figure.  The 
planes  driven  to  the  sides  each  have  a  finite  extent  and  a 
trapezoidal  shape.  The  planes  driven  down  under  the  indent¬ 
er  have  a  square  shape. 

Much  of  the  work  done  on  the  mdenter  goes  into  the 
energy  to  form  the  indicated  dislocation  configuration,  al¬ 
though  some  certainly  goes  into  heat.  There  are  two  major 
contributions  to  this  formation  energy.  The  first  is  the  enegy 
needed  to  generate  each  dislocation  as  an  isolated  entity,  and 
the  second  is  the  interaction  energy  among  these  disloca¬ 
tions.  Because  the  interaction  term  dominates  H,  approxi¬ 
mations  made  to  simplify  the  first  term  are  relatively  unim¬ 
portant.  The  extra  planes  driven  to  the  sides  of  the 
indenUtion  have  a  finite  extent;  accordingly  (in  this  idealized 
picture),  there  are  both  edge  dislocations  at  their  base  and 
screw  dislocations  associated  with  their  termination.  The 
square  planes  driven  below  the  indentation  have  edge  dislo¬ 
cations  around  the  sides  and  screw  dislocations  at  the 
comers  to  make  the  turns.  Moreover,  there  are  interactions 


1 


I  . 
t 


HG.  I .  Schematic  represenution  of  an  idealized  minimum-eneriy  disloea- 
lion  eonfiguraiion  produced  by  a  square  cross-section  indenter  in  a  hard¬ 
ness  measurement.  H'is  the  side  length  and  A  it  the  depth  of  the  indenution. 
In  this  ideal  case,  half  the  material  from  the  indented  region  is  displaced 
along  the  glide  planes  (indicated  by  the  dashed  arroersi  to  the  sides  and  half 
is  displaced  below  the  indenter. 
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among  the  dislocations,  which  can  produce  a  minimum-en* 
ergy  configuration.  For  the  arrangement  depicted  in  Fig.  1. 
it  always  cosu  energy  to  position  a  second  dislocation  on  a 
parallel  glide  plane  to  one  already  present.  However,  the 
magnitude  of  this  extra  energy  can  be  minimized  and,  for 
proper  configurations  of  the  dislocations,  there  are  attractive 
forces  along  the  glide  planes  that  will  tend  to  position  the 
dislocations  into  the  minimum-energy  configuration.  The 
minimum  configuration  arises  when  half  of  the  atoms  from 
the  indented  volume  go  respectively  to  the  side  and  below  the 
indenter.  Then,  in  both  regions,  the  maximum  angle  made 
between  successive  close-spaced  dislocation  lines  and  their 
glide  planes  =  ir/2  — /2,  as  shown  in  the  figure.  This  is 
the  minimum  realistic  energy  configuration.  If  the  disloca¬ 
tions  are  separated  more  than  shown  in  Fig.  1,  then  there  is 
more  volume  of  strained  material  and  the  interaction  energy 
would  be  larger  still. 

An  approximate  expression  for  the  energy  required  to 
indent  the  material  is^ 


tant  However,  care  most  be  taken  with  the  interaction 
terms:  Eq.  (1)  neglects  a  number  of  secondary  interactions, 
some  positive  and  others  negative;  these  will  be  added  later. 
The  principal  neglected  terms  are  the  interactions  between 
the  dislocations  in  the  different  regions  on  each  side  (posi¬ 
tive),  and  the  interactions  between  adjacent  side  and  bottom 
regions  (positive),  and  the  interr^on  between  the  opposite 
sides  in  the  bottom  regions  (negative).  Comparison  of  the 
results  with  experiment  will  indicate  how  important  these 
neglected  terms  are  likely  to  be.  The  length  L,  'tsW[N  —  i)/N 
for  i  from  1  to  N.  This  is  the  largest  length  of  the  side-inserted 
planes  and  its  choice  partially  accounts  for  interacitons 
between  the  otherwise  neglected  screw  dislocations.  The  dis¬ 
tance  =  r,  —  tj  is  given  by  2>/2F  |/  —  J\  for  /  and  j  ranging 
from  1  to  N  for  the  minimum-energy  configuration  uid  a 
tetrahedrally  bonded  semiconductor.  Finally,  in  the  indicat¬ 
ed  configuration,  ^  =  ir/2  —  d  /2.  Inserting  these  ex¬ 
pressions  into  Eq.  (1)  and  retaining  only  terms  of  order 
yields 


er  =  2a\  f  E,L,  +  ^  f  [min(/.,Z,)]  1  (1) 

Lib  1  2  1  J 

In  the  first  term,  using  an  isotropic  medium  approximation 
and  neglecting  core  terms,  the  energy  per  unit  length  to  form 
an  edge  dislocation  is^ 


Gb^  ,  R 

- - In — , 

4it(1— v)  ro 


(2) 


the  shear  coefiicient  is  G,  the  Burger’s  vector  d  /3  for  an 
indentation  along  a  <  100)  axis,  the  range  of  the  elastic  defor¬ 
mation  of  a  dislocation  R  is  taken  equal  to  fV  (for  want  of  a 
better  approximation),  ro—d  is  the  dislocation  core  radius, 
the  Poisson  ratio  is  v— 0.2  for  most  semiconductors,  and  L, 
is  the  length  of  the  /th  dislocation.  In  the  second  term,  is 
the  interaction  energy  per  unit  length  between  dislocations  t 
and  j.  Assuming  they  have  parallel  glide  planes  and  their 
Burger’s  vectors  have  the  same  sign,  is  given  by^ 


where  is  the  separation,  and  is  the  angle  that  a  line 
perpendicular  to  and  joining  the  dislocations  makes  with  the 
glide  plane,  as  shown  in  Fig.  1.  Because  the  various  disloca¬ 
tions  in  a  region  have  different  lengths,  the  net  interaction 
energy  is  approximated  by  multiplying  the  energy  per  unit 
length  by  the  length  of  the  shorter  one.  The  upper  limits  on 
the  sums  Naie  the  number  of  dislocations  in  one  region  (side 
or  bottom)  associated  with  one  edge.  For  the  minimum-ener¬ 
gy  configuration,  S  =  N„/2.  The  four  that  multiplies  the 
bracket  accounts  for  the  four  sides,  and  the  two  for  the  two 
regions  for  each  side. 

We  now  encounter  our  first  surprise.  As  we  can  see 
from  Eq.  (2)  and  (3),  and  £,  have  comparable  magnitudes. 
Because  there  are  approximately  N  ^  terms  in  the  interaction 
energy  sum,  only  N  terms  in  the  formation  energy  sum,  and 
Sisi  large  number  (A>  1 ),  the  interaction  energy  completely 
dominates  the  hardness.  In  fact,  N  is  typically  of  the  order 
10*.  Hence,  terms  owing  to  screw  dislocations,  core  energies, 
heat  dissipation  as  the  dislocations  propagate  to  their  places, 
and  other  effects  associated  with  the  first  term  are  unimpor- 
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One  can  also  get  a  number  for  the  hardness  of  a  disloca¬ 
tion  in  which  all  the  material  is  pushed  along  the  same  glide 
plane,  e.g.,  to  the  side,  to  the  bottom,  or  normal  to  the  face  of 
the  indenter  (a  possibility  not  depicted  in  Fig.  1 ).  In  this  case, 
the  factor  of  2  in  front  of  Eq.  1 1 1  is  removed,  N  =  and 
d  =  ir/2  — 17.  Then  a  higher  nonequilibrium  hardness  in  the 
context  of  this  model  (denoted  N,},  is  obtained. 

„  G  cot  i?  r  ,  ( cot  rf\  4  ,  .  2  a  1  ,r. 

"  1  " 
The  proper  answer  for  most  materials,  and  depending  on 
crystal  orienution,  probably  lies  somewhere  between 
and  H,.  For  an  indenter  with  ^  =  jr/4,  we  have 
=  0.0969(7/(1  -V)  and  f/, ///„,„  =2.39.  Harrison'  has 
shown  that  one  contribution  to  the  shear  coefficient  (actually 
C,,  —  Cjj)  is  C  =  2.38  fFmal/md^,  where  m  is  the  free- 
electron  mass,  o,  is  the  covalence,  =  Fj/(F|  +  V\]''^, 
and  d  is  the  bond  length.  We  will  approximate  G  by  this 
expression.  Using  this  G  and  v  =  0.2,  and  changing  the  di¬ 
mensions  to  those  in  terms  of  which  experimental  hardness 
numbers  are  customarily  quoted  gives  =  2.38  X 10* 
[a]/d  ’)  kgm/mm^,  where  d  is  in  angstroms.  (Calculated  val¬ 
ues  of  and  are  plotted  against  experimental  results 
in  Fig.  2  for  a  number  of  semiconductors. 

Figure  2  has  the  theoretical  /f„,„  and  /f,  values  con¬ 
nected  by  arrows  from  to  /f ,  for  each  compound,  plot¬ 
ted  as  a  function  of  the  corresponding  experimental  val¬ 
ues.*’  If  the  theory  were  perl  wt  and  the  experimental  values 
were  accurate,  the  points  would  fall  on  the  indicated  unity 
slope  line.  Several  conclusions  can  be  drawn.  First./,  the 
order  of  magnitude  of  the  predicted  and  measured  values  are 
the  same,  a  result  obtained  with  no  adjustable  parameters  in 
the  theory.  Secondly,  the  trends  from  one  compound  to  an¬ 
other  are  properly  given  by  the  theory.  Although  the 
values  are  generally  too  small,  they  fit  the  soft  materials  bet¬ 
ter,  and  the  /f ,  values  fit  the  harder  materials  better.  Thirdly, 
from  Eqs.  (1)  and  (3),  H  is  given  in  a  rough  but  revealing 
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FIG.  2.  ThcoreticaJ  vs  expeiimcntai  hardness  of  several  semkonduetors  in 
Refs.  2  and  9.  The  two  theoretical  values  for  each  semiconductor  are  con¬ 
nected  by  an  arrow  from  /f„,  to  /f,.  Perfect  agreement  would  correspond 
ic  points  being  on  the  unity-slope  solid  line. 

approximation  by  AN^E,jL  /Ah,  where  the  averge  disloca¬ 
tion  length  L  =  It' /2  and  £,^  is  the  average  dislocation  pair 
interaction  energy  per  unit  length.  Notice  that  N  (or  IT)  can¬ 
cels  from  this  expression;  thus,  H  is  independent  of  If'  (or  £ ) 
and  therefore,  H  is  truly  a  measure  of  the  properties  of  the 
material.  This  result  would  not  be  found  if  the  dislocation 
energies  ( «  )  were  to  dominate  H  rather  than  the  pair  inter¬ 
action  energy  (acJV’|.  Finally,  the  Berger's  vector  cancels 
from  the  leading  term  and  appears  only  in  the  argument  of 
the  logarithm  in  Eq.  (3).  Thus,  the  answers  are  also  insensi¬ 
tive  to  its  choice. 

Dislocations  are  often  found  in  materials  as  they  are 
grown.  Their  density  is  determined  by  the  thermal  and  me¬ 
chanical  stresses  to  which  they  are  subjected  in  the  growth 
process.  A  dislocation  constitutes  a  metastable  excitation 
relative  to  the  perfect  crystal  ground  state.  At  the  elevated 
growth  temperatures  and  temperature  gradient  behind  the 
growth  front,  the  number  of  dislocations  present  is  con¬ 
trolled  by  the  relative  rate  at  which  vacancies  anneal  or  con¬ 
dense  into  dislocations. The  dislocation  formation  rate  will 
be  slower  in  a  material  grown  at  the  same  temperature  if  £,  is 
higher.  If  an  alloy  is  formed  from  a  material  of  interest  and  a 
second  constituent  with  a  shorter  bond  length,  one  expects 
the  average  bond  energy  (and  thus  both  the  melt  temperature 
and  vacancy  formation  energy)  to  increase  proportional  to  a 
low  inverse  power  of  the  average  bond  length.*  Hence,  the 
equilibrium  vacancy  density  just  below  the  melting  point 
tends  to  be  the  same  in  lowest  order  for  all  materials,  inde¬ 
pendent  of  the  bond  lengths  of  the  constituents.  However, 
the  shear  coelBcient  and  dislocation  energy  per  unit  length 
will  increase  with  much  higher  inverse  powers  of  the  bond 
length.  Omsequently,  dislocation  densities  should  be  re¬ 
duced  in  such  alloys  relative  to  those  found  in  the  longer 
bond  length  pure  constituent.  This  expectation  is  confirmed 


in  the  recent  work  reported  on  Zni  _  „  Cd.Te  bulk  material.* 
The  best  CdTe  that  has  ever  been  grown  has  dislocation  den¬ 
sities  in  excess  of  3  X  lO’  cm  ~  The  addition  of  only  4%  Zn 
reduced  the  dislocation  count  to  less  than  5x10*  cm~*.  The 
ZnTe  bond  length  is  2.643  A,  while  the  CdTe  bond  length  is 
2.80S  A.  a  6%  difference.  This  6%  difference  in  bond  length 
translates  into  a  2%  difference  in  the  dislocation  energy  per 
unit  length  for  1  —  x  =  0.04.  Dislocation  energies  per  unit 
length  are  typicaUy  10  eV  per  lattice  spacing;  accordingly,  a 
2%  increase  can  be  expected  to  slow  their  formation  rate 
considerably. 

The  argument  just  presented  naturally  leads  to  a  strate¬ 
gy  for  decreasing  dislocations  other  semiconductors.  If  an 
alloy  is  made  of  the  material  of  interest  with  another  com¬ 
pound  with  a  shorter  bond  length,  then  the  dislocation  den¬ 
sity  should  be  reduced.  For  example,  this  suggests  that  the 
addition  of  a  small  amount  of  GaP  (d  =  2.359  A)  may  signifi¬ 
cantly  reduce  the  dislocation  density  of  bulk  grown  GaAs 
(d  s  2.448  A).  It  has  been  demonstrated  that  the  addition  of 
approximately  1%  GaN  (d  =  1.946  A)*-‘  or  of  a  10“  cm“* 
BAs(d  =  2.069  A)  concentration"  to  GaAs  can  yidd  a  large 
volume  of  dislocation-free  material.  An  InAs  additive  with 
its  longer  bond  length  (d  =  2.623  A)  serves  the  same  fimction 
indirectly,  by  causing  GaAs  bonds  in  its  neighborhood  to  be 
compressed.  This  indirect  mechanism  should  be  less  effec¬ 
tive  than  substituting  short  bond  length  additives. 

We  have  demonstrated  that  the  dislocation  energies 
and  hardness  of  tetrahedrally  bonded  semiconductors  are 
rapid  functions  of  the  reciprocal  of  the  bond  length.  This 
rapid  d  dependence  of  dislocation  energies  provides  a  ratio¬ 
nale  for  the  dramatic  decrease  of  the  dislocation  density  in 
bulk  grown  Znoa«  Cdo.M'Te  material  relative  to  that  found  in 
CdTe,  and  suggest  means  for  accomplishing  the  same  ends  in 
other  materials. 
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The  bond  length  and  energy  changes  of  the  constituents  of  alloys  relative  to  their  pure  crystal 
values  are  calculated  from  an  extension  of  Harrison’s  method.  It  is  demonstrated  that  the  already 
weak  HgTe  bonds  are  destabilized  by  adjacent  CdTe,  HgS,  or  HgSe,  but  are  stabilized  by  ZnTe.  It 
is  also  argued  that  dislocation  energies  and  the  hardness  of  semiconductors  vary  as  a  high  inverse 
power  of  the  bond  length  of  the  constituents.  Hence,  the  shorter  ZnTe  bond  as  ap  additive  should 
improve  the  structural  properties  of  HgTe  and  CdTe.  Experiments  that  support  these  predictions 
are  noted.  The  electronic  transport  properties  of  0. 1  eV  band  gap  HgZnTe  are  about  the  same  as 
those  of  HgCdTe,  and  the  structural  properties  of  the  Zn  compound  are  superior;  thus,  we 
conclude  that  HgZnTe  is  likely  to  be  the  better  material  for  IR  devices. 


I.  INTRODUCTION 

The  objective  of  this  work  is  to  understand  the  microscopic 
mechanisms  that  govern  the  stability  of  Hg,  .^Cd,Te  al¬ 
loys,  and  then  to  suggest  changes  in  the  material  that  im¬ 
prove  its  strength  without  adversely  affecting  its  electronic 
behavior.  In  pursuit  of  this  goal  we  have  extended  Harrison's 
bonding  theory,'  which  is  applicable  to  all  tetrahedral  struc¬ 
tured  semiconductors,  to  calculate  bond  length  and  energy 
changes  in  an  alloy — including  charge  shift  and  reconstruc¬ 
tion  effects — relative  to  their  pure  crystal  values,  vacancy 
formation  energies,  dislocation  energies,  and  hardness. 

In  this  paper,  we  review  the  experimental  situation  in  sev¬ 
eral  of  these  areas  and  compare  some  of  the  results  with 
theory  for  all  the  group  IV,  III-V  compound,  and  II- VI 
compound  semiconductors  and  their  alloys.  For  example, 
the  theory  properly  predicts  the  observed  inverse  ninth  pow¬ 
er  bond  length  (</  ~’)  dependence  of  the  hardness^  of  semi¬ 
conductors  and,  with  no  adjustable  parameters,  their  correct 
magnitude.  We  demonstrate  that  the  weak  HgTe  bond  is 
destabilized  by  alloying  it  with  CdTe,  HgS,  or  HgSe;  how¬ 
ever,  the  bond  is  stabilized  by  ZnTe.  Moreover,  because  the 
bond  length  of  ZnTe  (2.406  A)  is  14%  shorter  than  that  of 
HgTe  (2. 797  A)  or  CdTe(2.804  A),  the  dislocation  energy  per 
unit  length  and  hardness  of  the  alloys  Hg,  .^Zn„Te  and 
Cd,  _„Zn,Te  are  predicted  in  agreement  with  some  experi¬ 
ments  to  be  significantly  higher  than  those  of  the  compounds 
with  X  as  0.  Measurements  indicate  that  the  electron  an  ' 
hole  mobilities  of  Hg,  _,Zn,Te  with  x-~0. 16,  correspond¬ 
ing  to  a  0.1  eV  band  gap,  are  comparable  to  those  of 
Hg,.  ^CdxTe  for  x  =  0.2.^  This  is  to  be  expected,  because 
the  electron  effective  mass  in  a  narrow-gap  material  is  small; 
hence,  the  electron  wave  functions  at  r,c  are  distributed 
over  many  atoms,  with  the  result  that  alloy  scattering  rates 
are  small  at  the  band  edge.  The  hole  mobilities  are  large  for  a 
somewhat  different  reason:  The  valence  band  edge  is  mostly 
composed  of  Te  p  states;  however,  there  is  also  some  cation  p 
state  contribution.  Because  the p  state  energies  of  Hg  and  Zn 


are  nearly  the  same,  the  alloy  scattering  strength  at  the  va¬ 
lence  band  edge  is  small.  We  conclude  that,  because  of  its 
structural  advantage,  HgZnTe  may  be  superior  to  HgCdTe 
for  infrared  devices 

II.  ALLOY  MODIFICATION  OF  THE  BOND  ENERGY 
AND  LENGTH 

A.  Pure  compounds 

The  first  task  is  to  recast  Harrison’s  bond  energy  and  bond 
length  formalism  into  a  structure  suitable  for  generalization 
to  alloys.  Focus  attention  on  the  nth  bond  of  a  pure  zinc 
blende  structured  compound  semiconductor;  for  the  time 
being,  suppress  any  notation  identifying  it.  Then,  in  a  Il-VI 
compound,  the  bond  energy  is': 


+  +  yji 


+  2  +c.)  +  :iT- 

(-0  a 

where  the  first  term  is  the  energy  per  bond  needed  to  transfer 
two  electrons  from  ap  state  on  the  group  VI  atom  (the  anioni 
to  a  p  state  on  the  group  II  atom  (the  cation),  so  that  both 
start  with  four  electrons  (ultimately  in  the  final  bonding  ar¬ 
rangement  there  is  a  net  electron  transfer  from  the  cation  to 
the  anion),  the  second  term  is  the  promotion  energy  per  bond 
to  forr*  tp^  hybrids  on  both  atoms,  the  third  term  is  the  bond 
formation  energy  owing  to  the  covalent  and  ionic  terms,  the 
fourth  term  is  Harrison’s  metallization  energy,  and  the  fifth 
term  is  the  repulsive  overlap  energy  that  prevents  the  bonds 
from  collapsing.  The  various  symbols  are  defined  below; 

e,i  I  is  the  {“','0°"}  p  state  correlated  term  values,' 
i  I  =(f;>i  1 1  ~  1 1 1/“*  hopping  integrals 

between  two  adjacent  bonds  coupled  through  a  }• 
,  I  are  the  s  state  correlated  term  values, 

1^2=  —  24.5/</*(eV)  is  the  covalent  energy, 
d  is  the  bond  length. 
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K,»(6*-6_)/2. 

f  I  ±  I  *(f.i  ±  I  +  J I  )/4  are  the  sp'  hybrid  energies. 

.  1  *  i  (Uii  * ,  Ui,  ,  ,  ^ ,  )/(e?  -  ei) 

./-I 

is  the  contribution  to  the  reduction  of  the  bonding  state  ener¬ 
gy  calculated  by  second-order  perturbation  theory  of  the 
bond  in  question  (denoted  by  superscript  oj,  owing  to  its  in¬ 
teraction  with  the  three  neighboring  (i  =1-3)  antibonding 
states  sharing  the  same 

=VH-“r)/2  =  C/.-. 

u,.  = 

are  the  probability  amplitudes  of  finding  the  cation  and  an¬ 
ion  rp'-hybrid  contributions  to  the  bonding  and  antibonding 
states,  e.g.,  ^  is  the  probability  amplitude  of  finding  the 

cation'  hybrid  in  the  bonding  state,  etc., 

a^^Vf/yjVl  +  VI  is  the  polarity. 
a^^\y,\/yjvl  -(-  is  the  covalency, 

f,:, 


+  f- 
2 


are  the  bonding  and  antibonding  one  electron  energy  states 
owing  only  to  covalent  and  ionic  interactions, 

.  .  ..  ,,, 

- - — - - .  I -1.2.3 

are  the  contributions  to  the  bond  in  question  arising  from  its 
antibonding  state  interactions  with  a  neighboring  bonding 
sute  /,  and  y  is  an  adjusuble  coefficient  chosen  so  the  mea¬ 
sured  bond  lengths  are  reproduced. 

Harrison  has  written  Eq.  (1)  without  explicit  reference 
to  the  surrounding  bonds.  Equation  ( 1 1  in  its  present  form  is 
rather  easily  used  for  alloy  calculations  by  replacing  host 
atom  terms  by  one  appropriate  to  the  alloy. 

The  parameter  y  for  each  compound  is  fixed  by  the 
condition  that  there  is  no  net  force  on  atoms  at  their  equilib¬ 
rium  lattice  positions.  This  condition  is  satisfied  if  the  gradi¬ 
ent  of  the  total  energy  with  respect  to  all  the  bond  lengths 
vanishes: 

where£r  =  X^*- 

tl 

The  condition  is  equivalent  to  the  requirement  that  the  net 
bond  tension  T„  +  T„^  for  the  nth  bond  equals  the  repul¬ 
sive  force 

=  (2) 

a  - 


T„»^(-Uv{^y\)  =  ^,4y\^y\/d.  o) 

da.  ,-0 


Table  1.  Correlated  term  values  and  bond  energies  for  three  II-VI  com¬ 
pounds.  The  experimental  bond  lengths,  polarity,  and  metallic  contribution 
to  the  bond  energy  are  also  plotted 
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Correlated  term  values  |eV| 

Atom 

5 

P 

sp' 

Zn 

-  10.224 

-  4,920 

-  6.246 

Cd 

-  9.611 

-  4.784 

-  3,991 

Hg 

-  10,946 

-  4.872 

-  6.391 

Te 

-  19,620 

-  9.824 

-  12.273 

d(Ar 

Bond  energies 

Compound 

“r 

^metallic 

£»(eV) 

ZnTe 

2,643 

0,632 

-  0.625 

-  0.985 

CdTe 

2.g0S 

0,710 

-0.466 

-0.944 

HgTe 

2,797 

0,683 

-  0.614 

-0.480 

'Reference  12. 

The  last  identity  in  Eq.  (4)  defines  the  indicated  bond  ten¬ 
sions.  Collecting  all  the  results  from  the  earlier  paragraphs, 
one  finds  (dropping  the  bond  designation  n): 

Tl-  =  --j[(l+o,)a, -t-^i^^^]<;-.(5a) 

L  ^,-e'a  \ 

Tl.  =  Tl.  =  Tl. 

=  ,5b, 

r;-=.i  C-.  (5c)  Q 

■yl  _  yJ  _  yJ 

®  L  ^  • 

Finally  combining  Eqs.  (l)-(5),  the  bond  energy  is 
£,=  _l(l+aJ)(F’-)-K5)''^ 

8  '  (vl  +  y\)''^ 

+  (6) 

Equation  (6)  was  first  derived  by  Harrison  using  a  much 
simpler  argument. '  The  bond  energies  quoted  in  Table  I  and 
on  Fig.  S  are  calculated  from  this  expression,  using  as  input 
the  correlated  term  values  from  Table  I  and  the  experimental 
bond  lengths. 

B.  Alloy  calculations 

There  are  several  layers  of  sophistication  that  can  be  used 
to  compute  the  bond  energy  and  length  changes  between 
pure  compounds  and  their  alloys.  However,  all  the  calcula¬ 
tions  are  motivated  by  the  EXAFS  experimentol  result 
found  first  for  Ga,  _  In,  As  and  later  in  other  compounds.^ 

It  was  found  that,  while  the  average  bond  length  in  the  alloy 
follows  Vegard’s  rule  and  varies  linearly  between  those  of  the 
pure  constituents,  the  individual  bond  lengths  of  GaAs  and 
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InAs  arc  iwiriy  the  same  as  those  of  the  pure  materials — 
changug  by  only  about  *%  The  shorter  GaAs  bond  in- 
creaiai  its  length  when  it  is  an  unpurity  in  an  InAs  host  by 
about  *%.  and  the  longer  InAs  bond  shrinks  by  about  4%  in 
the  opposite  extreme  where  it  is  the  impurity  in  GaAs.  The 
lattice  fits  together  by  having  the  metal  sublattice  nearly  re¬ 
tain  its  fee  structure  with  a  lattice  constant  following  Ve- 
gard’s  rule,  and  each  As  atom  adjusts  its  local  position  and 
accommodates  to  which  metal  atoms  occupy  its  four  neigh¬ 
boring  sites.  Thus  there  are  five  local  arrangements;  in  the 
first  an  As  atom  is  surrounded  by  four  Ga  atoms  [a 
Ga(4)In(o)  configuration],  in  the  second  there  are  three  Ga 
atoms  and  one  In  [a  Ga|3)In(l)  configuration],  and  the  last 
three  configurations  are  Ga|2)In|2),  Gall)In(3),  and 
Galo)In|4).  In  a  general  cation  A,  alloy  or  anion 

CA,  alloy,  the  generalizations  are  obvious,  and  one 
can  examine  A|4  —  rf)^ri)\7j  »  1,2,3,4  configurations. 

I 

1.  Method  1 

The  lowest  level  of  approximation  that  takes  account  of 
these  findings  calculates  [from  Eqs.  (1)  and  (2)]  a  given 
bond’s  energy  and  length  shift  when  its  surrounding  bonds 
have  different  constituents.  The  calculation  is  done  assum¬ 
ing  that  the  surrounding  bonds  retain  their  pure  crystal  bond 
lengths.  This  model  is  somewhat  unrealistic,  in  that  it  ig¬ 
nores  the  strains  that  must  be  present  to  allow  the  lattice  to 
fit  together.  However,  it  is  nonetheless  useful  because  it 
treats  the  largest  terms,  viz.  those  arising  from  charge  trans- 
fer  between  the  substituted  species  in  the  alloy,  and  results  in 
analytic  expressions,  while  the  successively  more  complicat¬ 
ed  treatments  must  be  done  numerically.  It  is  also  useful 
later  because,  taken  together  with  the  more  complete  calcu¬ 
lations,  it  allows  one  to  judge  the  relative  contribution  from 
different  physical  effects.  From  Eq.  1 1 1  in  a  cation-substitut¬ 
ed  alloy  one  can  deduce  an  expression  for  the  energy  shift  of 
an  AC  bond  when  t;  surrounding  bonds  arc  replaced  by  BC. 
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where  the  partial  derivative  is  taken  with  respect  to  the  cen¬ 
tral  bond  length  |a  term  designated  with  a  superscript  ol  and 
in  this  simple  approximation 

=  77  =  0.1,2,3.  (8) 
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FiC-  1.  Schematic  view  of  vanous  bond  configurations,  lai  Ai3)Bi  1 1  configu> 
ration,  ibi  Al2)Bl2l  configuration,  and  ic)  A(3lBl  M  configuration  embedded 
in  a  virtual  crysuj. 


The  bond  length  shifts  are  found  from  Eq.  (2).  The  modifica¬ 
tion  of  these  expressions  for  an  anion-substituted  alloy  fol¬ 
lows  from  symmetry. 

2  Method  2 

A  more  realistic  calculation  includes  strain  and  permits 
the  lattice  to  fit  together.  The  simplest  non-self-consistent, 
near-neighbor  version  of  this  calculation  proceeds  as  fol¬ 
lows.  In  a  cation-substituted  alloy  A,.,B.C 
(e.g..  Hg,_  ,Cd,Te  or  Zn,_,Hg,Te),  the  cations  are  as¬ 
sumed  to  occupy  their  regular  fee  sublattice  sites,  while  the 
anions  accommodate  to  their  local  configuration.  The  aver¬ 
age  lattice  constant  is  taken  to  follow  Vegard’s  rule: 

5  =  ( 1  -  xkf 0  (AC)  -(-  xdo(BC).  19) 

where  d,,  and  dglBC)  are  the  pure  crystal  lattice  constants. 
Then  in  an  AI3)BM)  configuration  [see  Fig.  Hal],  the  anion 
will  be  displaced  along  the  BC  bona.  For  definiteness,  sup¬ 
posed  (AC)  >  d  (BC);  then  the  anion  will  be  shifted  away  from 
the  tetrahedral  site  toward  the  B  atom.  Because 
dolAC)>d>dn(BC),  thed  (BC)  bond  will  be  elongated  and  the 
d  (AC)  bonds  will  also  be  locally  elongated  in  a  strained  con¬ 
figuration.  We  can  for  this  case  [see  Fig.  l(al]  wnte 

d(BC)  =  (l -(5)d  (101 

and  for  the  three  AC  bonds 


d(AC)  =  . 
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1  -h  —  <5  -  d. 
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In  an  A(2)B(2)  configuration  [see  Fig.  1(b)],  we  have  for  both 
BC  bonds 


(12) 

and  for  both  AC  bonds 

d(AQ^  ll+jS  +-jS^d. 

(13) 

a 

The  values  of  5  are  determined  by  the  condition  dEr/ 
dSxsQ,  or 

,““0  1  "*  dl  )  dS 

(14) 

Equation  ( 14)  is  solved  by 

( 1 )  Assigning  the  s  their  pure  crystal  values. 

(2)  Supposing  that  the  four  tetrahedrally  configured  an¬ 
ions  are  embedded  in  the  virtual  crystal  [see  Fig.  1(c)]. 

(3)  For  each  x,  5  is  varied  and  the  bond  tensions  are 
calculated  and  inserted  into  Eq.  (14)  until  a  value  of  £  is 
found  that  satisfies  this  equation. 

The  computer  outputs  are  S  and  the  AC  and  BC  bond  ener¬ 
gies.  This  method  can  be  extended  to  take  explicit  account  of 
the  bond  length  variations  of  larger  clusters,  and  d  can  be 
calculated  explicitly  (instead  of  using  Vegard’s  rule)  so  the 
results  can  be  made  to  be  self-consistent.  Also,  effects  arising 
from  bond  angle  distortions  can  be  included.  The  major  ef¬ 
fects  that  are  currently  included  in  the  calculation  arise  from 
the  bond  length  dependence  of  Kj,  and  charge  shifts  driven 


Fig.  3.  Bond  energy  as  a  Tunction  of  concentration  x  for  the  Hg,  .,Cd,Te 
alloy  in  different  configurations.  The  short  dashed  lines  are  for  the  CdTe 
bonds  and  the  long-short  dashed  lines  represent  the  behavior  of  the  HgTe 
bonds.  The  solid  lines  are  the  concentrauon  weighted  averages  of  the  respec¬ 
tive  bond  energies. 


Fic.  2.  Bond  lengths  as  a  fiincnon  of  concentration  Jt  for  the  Hg, . ,  Cd,  Te 
alloy  m  different  configurations.  The  short  dashed  lines  are  for  the  CdTe 
bonds  and  the  long-shon  dashed  lines  represent  the  behavior  of  the  HgTe 
bonds. 


by  hybrid  energy  differences  and  coupled  through  the  metal¬ 
lization  terms.  However,  K;  is  also  sensitive  to  bond  angle 
distortions,  and  these  effects  are  not  included  in  the  present 
results.  We  expect  their  inclusion  will  modify  the  quantita¬ 
tive  results  by  10%-30%,  but  not  the  trends. 

C.  Results 

The  correlated  atomic  term  values*  for  the  elements  in  the 
more  common  II- VI  compounds  are  listed  in  Table  I,  along 
with  the  bond  energies  of  their  nine  compounds.  The  valence 
s  state  energies  for  Zn  are  large,  become  smaller  as  for  Cd, 
and  then  (in  an  unusual  occurrence)  the  trend  reverses  and 
the  Hgr  levels  are  deeper  again.  This  is  caused  by  the  relativ¬ 
istic  terms  that  become  important  in  Hg.  The  p  state  levels 
have  a  similar  trend;  however,  the  differences  among  the 
atoms  are  much  smaller.  The  levels  of  the  anions  are  deeper 
than  those  of  the  cations.  Thus,  there  is  a  net  electron  trans¬ 
fer  from  the  cations  to  the  anions,  which  is  responsible  for 
the  ionic  contribution  to  the  bonding.  The  polarity,  also  list¬ 
ed  in  Table  I,  reflects  the  relative  contribution  of  ionic  and 
covalent  character  to  the  bonds.  Notice  that  is  small  for 
ZnTe  larger  for  CdTe,  and  smaller  again  for  HgTe.  By  con¬ 
trast,  the  metallic  contribution  E„  is  large  for  ZnTe  and 
HgTe  and  small  for  CdTe.  In  fact,  HgTe  would  not  be  bound 
if  it  were  not  for  the  metallic  contribution  to  the  bond,  not  a 
terribly  surprising  result,  because  HgTe  is  a  semimetal.  The 
bond  lengths  of  CdTe  and  HgTe  are  very  nearly  the  same, 
but  this  is  an  accident  resulting  from  the  balance  of  contribu- 
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tions  to  their  bonds;  CdTe  is  dominated  by  the  covalent  and 
ionic  contributions  to  its  bond,  while  HgTe  has  bonds  that 
are  mostly  covalent  and  metallic.  We  should  therefore  not  be 
too  surprised  if  their  alloy  displays  some  unusual  bond 
length  variations. 

Figures!  and  3  illustrate  the  Hg,  ^.Cd^Te  alloy  concen¬ 
tration  variations  of  the  length  and  energy  for  the  five  config¬ 
urations.  These  results  were  calculated  using  the  second 
method  discussed  in  Sec.  II  B.  In  Fig.  2  we  see  the  unexpect¬ 
ed  result  that  the  shorter  HgTe  bond  becomes  still  shorter  in 
the  alloy,  while  the  CdTe  bond  length  increases.  Moreover, 
the  changes  are  large  compared  to  the  original  length  differ¬ 
ence.  Thus,  one  can  expect  to  find  local  microstrains,  even  in 
this  nearly  lattice-matched  material.  The  bond  energy  varia¬ 
tions  shown  in  Fig.  3  also  display  an  important  result.  The 
already  weak  HgTe  bond  is  destabilized  by  the  presence  of 
Cd.  The  flat  concentration  variation  of  the  bond  lengths  and 
energies  is  peculiar  to  HgCdTe;  much  more  structured  beha¬ 
viors  are  exhibited  by  other  compounds. 

Several  considerations  to  keep  in  mind  when  dealing  with 
alloys  are  depicted  schematically  in  Fig.  4  for  a  Hg(3)Cd(l) 
configuration.  Due  to  the  charge  shift  alone,  the  HgTe  bond 
shrinks  by  0.030  A,  the  minimum  in  the  configuration  dia¬ 
gram  moves  to  a  smaller  bond  length  and  the  depth  of  the 
well  is  decreased.  Then  the  bond  is  stretched  by  0.0 11  A  to  its 
final  length.  Thus,  the  bond  is  not  centered  at  a  local  mini- 


FiC.  S.  AC  bond  energy  shifts  for  the  II-VI  compounds  m  All  IBI3I  configu¬ 
rations  expressed  as  a  percenuge  of  the  pure  crystal  bond  energy.  The  B 
atom  IS  designated  in  the  line  that  shows  the  shift.  The  pure  crystal  AC  bond 
energies  and  lengths  are  given  at  the  top  of  the  figure.  These  numbers  were 
calculated  by  method  I  in  Sec.  II  B 
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Fic.  4  Schematic  picture  of  a  potential  configuration  diagram  for  the  HgTe 
and  CdTe  bonds  in  an  Hgi3ICdlli  configuration  For  the  HgTe  bond  the 
energy  and  position  of  the  bond  center  are  calculated  by  method  I  of  Sec. 
II  B  and  the  final  stretched  position  and  energy  by  method  2.  The  results 
quoted  for  the  CdTe  bond  are  those  determined  by  method  2. 


mum.  For  example,  this  implies  that  the  force  constant  of 
local  phonon  modes  will  be  related  to  the  second  derivative 
of  this  configuration  diagram  evaluated  at  the  strained  posi¬ 
tion  away  from  the  well  minimum.  This  effect  is  even  more 
pronounced  for  materials  with  a  larger  lattice  mismatch. 
Notice  that  the  CdTe  bond  is  also  stretched  relative  to  its 
local  minimum  in  this  configuration. 

The  result  for  bond  energy  variations  for  all  the  II-Vl 
compound  alloys  is  gathered  in  Fig.  S.  Here  the  energy  shifts 
of  an  AC  bond  in  a  A|  1  )B|  3)  configuration  ^Ei  =  ME  ^  are 
presented  as  a  percentage  of  the  pure  material  bond  energy 
for  cation  substitute  alloys  along  with  equivalent  results  for 
anion  substantial  alloys,  e.g.,  HgS,  _ .  Te, .  The  energy  shifts 
are  calculated  following  the  first  method  presented  in  Sec. 
II  B.  We  have  not  had  the  opportunity  to  run  all  the  com¬ 
pounds  following  the  more  sophisticated  calculational  pro¬ 
cedure.  The  bond  at  the  top  of  a  column  is  the  one  whose 
energy  shift  is  calculated,  and  the  other  constituent  is  identi¬ 
fied  in  the  line  that  indicates  the  energy  shift.  A  negative 
ME^/E^  is  destabilizing  and  a  positive  one  stabilizing.  Fo¬ 
cus  on  HgTe.  It  is  destabilized  by  CdTe,  HgS,  and  HgSe,  but 
is  stabilized  by  ZnTe.  The  CdTe  destabilizes  the  HgTe  bond 
by  18%  according  to  the  simple  calculation  and  by  13% 
according  to  the  better  one.  The  difference  between  Zn  and 
Cd  occurs  because  in  the  Cd  case  there  is  a  net  electron  trans¬ 
fer  from  the  Cd  to  the  Hg,  while  in  the  case  of  Zn  there  is 
practically  no  such  transfer.  These  extra  transferred  elec¬ 
trons  must  occupy  antibonding  states,  because  the  bonding 
states  on  the  HgTe  are  full,  so  the  net  effect  i  including  other 
energy  state  shifts)  is  to  destabilize  the  HgTe  bond.  This  pre¬ 
dicted  destabilization  of  the  HgTe  bond  in  HgCdTe  alloys 
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FlC.  6.  Theoretical  VI  expcnmeatal  hardaem.  The  two  valuci  at  the  top  ind 
bottom  of  the  vertieal  liaes  rtrei|iiiatcd  for  each  compound  were 
for  two  different  poeable  ditloeatioo  arrangementt.  The  lower  one  oorre- 
qioodt  to  the  minhnum  energy  configuration  and  the  upper  one  to  a  likely 
higher  energy  configuration.  The  eaperimental  numben  were  taken  from 
Ref.  2. 

agrees  with  the  d  state  energy  shift  reported  by  Spicer  et  al.  at 
last  year’s  workshop,*  and  by  the  oxygen  uptake  experi¬ 
ments  presented  by  the  same  group  at  this  conference.'' 

III.  DISLOCATION  ENERGIES  AND  HARDNESS 

The  dislocation  energies  and  hardness  of  semiconductor 
are  calculated  by  an  extension  of  Harrison’s  method.'  It  is 
demonstrated,  in  agreement  with  experiment,  that  disloca¬ 
tion  energies  per  unit  length,  depending  on  the  covalency, 
are  proportional  tod  where d is  the  bond  length  and  hard¬ 
ness  is  proportional  to  d  to  d  ~  The  low  powers  are  for 
pure  covalent  materiak.  The  hardness  is  relat^  to  the  inter¬ 
action  energies  among  dislocations.  The  detailed  theory  will 
be  published  elsewhere.* 

Figure  6  shows  the  theoretical  hardness  of  several  semi¬ 
conductors  plotted  against  experimental  values.*  The  theory 
is  calculated  for  two  different  dislocation  configurations 
caused  by  the  extra  material  pushed  by  the  indentor  into  the 
crystal.  The  lower  end  of  each  vertical  line  is  a  realistic  hard¬ 
ness  for  the  minimum  energy  configuration,  and  the  upper 
ends  of  the  lines  are  the  values  found  for  a  higher-energy 
configurati'n.  In  most  cases,  one  would  expect  real  configu¬ 
rations  to  have  hardnesses  lying  somewhere  between  these 
two  values,  once  the  effect  of  heat  generated  by  the  indention 
process  is  subtracted  from  the  experimental  numbers.  This 
theory  has  no  adjusuble  parameters,  yet  it  predicu  the  right 
trends  and  magnitudes. 

IV.  NATIVE  IMPERFECTIONS  GROWN  INTO  ALLOYS  . 

There  are  two  kinds  of  imperfections  we  shall  discuss:  va¬ 
cancies  and  dislocations.  At  present,  the  conclusions  are  ten- 
utive,  because  the  detailed  theory  is  still  incomplete.  Va- 

JL  Vac.  M.  TsennoL  A  VoL  3,  No.  1,  Jan/ffeb  INS 


cancy  formations  energies  and  the  melting  temperature  are 
both  proportional  to  the  bond  energies  of  the  cnn«titiientf 
Thus,  one  would  expect  the  vacancy  density  just  below  the 
melting  point  to  be  about  the  same  for  all  materials  with  a 
similar  coordination  number;  in  our  semiconductor  case,  the 
coordination  number  is  four.  Dislocation  energies  are  most¬ 
ly  owing  to  long-range  bond  angle  distortions*  and  accord¬ 
ingly  are  sensitive  to  the  sheer  coefficients,  which  vary 
roughly  asd  for  many  compounds.'  Based  on  this  general 
argument,  one  would  expect  to  find  fewer  dislocations  in 
materials  with  short  bond  lengths. 

Dislocations  appear  in  bulk  grown  through  a 

sequence  of  steps.  "*  We  shall  only  rii«eii<t.<e  the  simplest 
where  there  are  no  mechanical  stresses  on  the  growing  mate¬ 
rial.  In  the  temperature  gradient  behind  the  growth  front, 
there  is  usually  a  supersaturated  vacancy  content  These  ex¬ 
cess  vacancies  may  diffuse  to  surfaces  where  they  do  no 
harm,  or  they  may  condense  into  vacancy  clusters  rh«r  can 
then  organize  into  dislocation  loops.  If  a  loop  is  not  paralM 
to  the  growth  front  it  can  climb  toward  the  front  by  absorb¬ 
ing  more  vacancies.  If  it  climbs  Cuter  th«n  the  growth  pro¬ 
gresses,  then  the  loop  will  reach  the  front  It  can  then  spread 
and  grow  along  with  the  crystal.  Observed  from  the  growth 
surface,  such  a  loop  will  appear  to  be  two  edge  dislocations. 
If  this  dislocation  formation  process  can  be  inhibited  in  any 
of  its  steps,  the  resulting  material  will  have  fewer  disloca¬ 
tions. 

Zn  in  Zn,  _  .Cd^Te  or  Zn,  _  ,Hg;,Te  wiU  serve  this  func¬ 
tion  in  two  ways.  The  first  is  that  Za  stabilizes  the  weaker 
bond  in  these  alloys;  accordingly,  vacancy  formation  ener¬ 
gies  in  the  alloy  should  be  larger  than  those  in  the  pure  con¬ 
stituent  CdTe  or  HgTe.  The  melting  temperature  «ill  also  be 
raised  somewhat;  which  effect  “wins”  is  a  matter  still  to  be 
decided.  The  second  effect  is  clear  cut;  The  shorter  bond 
lengths  of  the  ZnTe  will  introduce  stiff  struts  into  the  system 
that  will  increase  the  dislocation  energies  and  thus  should 
inhibit  dislocation  loop  formation  and  subsequent  climb. 
The  large  dislocation  density  reduction  found  by  Bell  and 
Sen  in  Zn,  _.Cd.Te  with  only  4%  Zn  supports  this  specula¬ 
tion." 

V.  CONCLUSIONS 

We  have  argued  that,  in  an  alloy,  the  bond  energies  and 
bond  lengths  arc  modified  from  their  pure  crystal  values. 
Bond  energies  can  be  strengthened  or  weakened  in  an  alloy 
depending  on  the  relative  state  energies  and  the  conse¬ 
quent  net  electron  shift  between  the  constituents.  For  most 
compounds,  the  average  bond  length  modifications  for  each 
type  of  bond  are  genetally  small  compared  to  those  suggest¬ 
ed  by  Vegard’s  rule,  in  agreement  with  experiment.  How¬ 
ever,  HgCdTe  is  an  exceptional  case,  because  the  bond 
lengths  for  HgTe  and  CdTe  are  nearly  the  same  (by  acci¬ 
dent).  The  bond  length  and  en^gy  shifts  also  depend  sensi¬ 
tively  on  the  local  configuration  and  the  bonds  are  locally 
strained.  These  effects  must  be  taken  into  consideration  in  a 
proper  theory  of  alloy  vacancy  formation  energies,  phonon 
frequencies,  dislocation  energies,  etc. 

Hie  conclusions  most  relevant  to  the  initial  question 
posed  in  the  introduction  are  that  for  the  narrow  gap  allows; 
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Cd  destabilizes  the  already  weak  HgTe  bond  and  has 
little  effect  on  the  dislocation  energy. 

Zn  both  stabilizes  the  HgTe  bond  and  increases  the  dis¬ 
location  energy. 

Zn  also  decreases  the  dislocation  density  of  CdTe  to 
improve  substrates. 

HgZnTe  electron  and  hole  mobilities  are  about  the 
same  as  those  of  comparable  band  gap  HgCdTe. 
HgZnTe  is  likely  to  be  a  superior  material  to  HgCdTe 
for  detectors  operating  in  the  8-12  /zm  regime. 


"'niis  work  was  supponed  in  part  by  AFOSR  Contract  No.  F496280-81-K- 
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,  Bond'length  relaxation  in  pseudobinary  alloys 
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Arden  Sher 
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(Received  17  September  1984) 

The  bond-length  relaxation  in  paeudobinary  alloys  can  be  predicted  by  a  simple  radial  force  model.  In 
tetrahedral  structure  alloys  the  bond-length  deviation  of  the  solute  in  the  dilute  solution  is  a  quarter  of  the 
bond-length  difference  between  the  two  components.  This  result  agrees  with  the  experimental  work  done 
on  the  Gai.^ln^As  system  performed  by  Mikkelsen  and  Boyce. 


In  calculating  pseudobinary  alloy  properties,  the  nearest- 
neighbor  disunce  is  a  very  important  parameter.  For  many 
years  people  have  used  the  virtual-crysul  approximation 
(VCA)  which  assumes  that  all  atoms  occupy  the  average  lat¬ 
tice  positions.  In  1982  Mikkelson  and  Boyce  used  extended 
x-ray-absorption  fine  structure  (EXAFS)  to  study  the  atom¬ 
ic  scale  structure  of  the  Gai-j,In.As  system.'  They  found 
that  the  Ga-As  and  In-As  bond  lengths  did  not  follow  the 
VCA;  instead,  they  only  slightly  deviated  from  their  natural 
bond  lengths  even  in  the  dilute  solutions,  and  the  magni¬ 
tude  of  the  bond-length  deviation  at  the  dilute  limit  was 
about  a  quarter  of  the  difference  of  the  natural  bond  lengths 
of  the  two  components  (GaAs  and  In  As).  These  results  can 
be  explained  by  a  simple  radial  force  constant  model  in 
-hich  the  weaker  angular  forces  are  neglected. 

Consider  a  compound  GaAs  in  which  one  of  the  Ga 
atoms  is  replaced  by  an  In  atom.  This  configuration 
represents  the  dilute  limit  of  In-As  in  a  GaAs  host  (see  Fig. 
1).  In  this  configuration,  four  In-As  bonds  are  identical 
and  the  tetrahedral  symmetry  is  preserved  at  the  In  atom. 
We  regard  any  difference  between  Ga  and  In  as  a  first-order 
difference.  Let  us  call  dcuu  the  natural  bond  length  of  the 
host  and  K  the  stretching  force  constant  for  the  host.  Ad¬ 
ding  In  will  change  K  by  an  amount  we  regard  as  first  order, 
but  that  alone  will  not  give  the  distonion.  However,  be¬ 
cause  the  natural  bond  length  of  InAs,  is  larger  than 
^otAa  by  a  first-order  amount  the  As  atoms  around  the  In 
impurity  move  outward  by  a  first-order  distance;  we  hold 
the  second  neighbors  fixed  and  to  first  order  can  take  K  to 
be  the  value  for  the  host.  Let  us  call  u  the  outward  relaxa¬ 
tion  displacement  of  the  As.  Then  the  amount  that 

the  final  In-As  bond  length  deviates  from  is 

dotAS  ~  ^inA«  7- u.  The  bond  lengths  of  the  neighboring  Ga¬ 
As  are  decreased  by  —  u  '^sB  equal  to  —  u/3  'o  first  order. 
The  energy  associated  with  the  distortion  is 

E-4t(k/2)(do.A.-dinA.+«)’  +  3(k/2)(u/3)n  . 

By  minimizng  the  energy  with  respect  to  u,  we  obtain 

r  the  deviation  of  the  In-As  bond  length 

^tf|«Ai  “  T  ^^foaAi  “  d|nA») 


Similarly,  the  bond-length  deviation  of  dilute  Ga-As  in  an 
InAs  host  is 

AdG»A»  "  T  (  d  1ba»  ~  do»Ai )  • 

This  predicts  exactly  Mikkelsen  and  Boyce’s  result.'  Note 
that  this  displacement  leaves  all  metallic  atoms  in  their  ori¬ 
ginal  sites;  only  the  As  atoms  move.  To  a  good  approxima¬ 
tion  this  is  true  in  the  concentrated  alloy,  where  the  metal 
atoms  form  a  face-centered-cubic  lattice  of  the  weighted 
average  lattice  parameter  and  the  As  atoms  move  to  fit  local 
bond  lengths. 

One  can  apply  this  principle  to  the  rocksalt-structure  sys¬ 
tems,  e.g.,  Ki.^RbxBr,  and  predict  that  the  bond-length  de¬ 
viation  of  the  solute  in  the  dilute  solution  is  half  of  the 


H  u  h 


(b) 


FIG.  I.  (a)  An  indium  atom  substitutes  for  gallium  in  GaAs 
the  tetrahedral  symmetry  is  preserved  at  the  In  atom  (b)  The 
displacement  of  As  atoms,  u.  results  in  the  Ga-As  bonds  being 
decreased  by  u  cost).  The  In-As  bond  lengths  are  varied  by  </c,ai 
+  "-d|„A. 
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difference  between  the  two  components: 

Adubii— ■ 

However,  because  of  the  lattice  geometry,  the  approxima¬ 
tion  that  all  the  cations  stay  in  their  original  sites  might  not 
be  as  good.  In  rocksalt  structure,  the  restoring  force  on  the 
displaced  first  neighbor  results  from  the  atom  right  behind 


■J.  C.  Mikkelsen.  Jr.  and  J.  B.  Boyce,  Phys  Rev  B  2g.  7130  (19g3) 


the  first  neighbor;  one  would  expect  that  the  relaxation  ex¬ 
tends  more  than  one  interatomic  distance  and  the  value  of 
the  bond-length  deviation  would  be  somewhat  smaller  than 
what  is  predicted  by  this  model.  This  agrees  with  the  expe' 
imental  result.’ ' 

We  thank  Dr.  Alex  Zunger  for  valuable  discussions  This 
work  was  supported  by  the  Defense  Advanced  Research 
Projects  Agency  under  Contract  No.  800989-B8 


^J.  B  Boyce  and  J  C.  Mikkelsen,  Jr  .  Bull  Am  Phys  Soc  29.  371 
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Binding  Energy  and  Spectral  Width  of  Si  2pCore  Excitons  in  SixGei.,  Alloys 
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f  and 


A.-B.  Chen‘*‘ 
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(Received  24  January  I98S) 

A  calculation  is  presented  to  explain  the  anomalous  experimental  behavior  of  the  Si  2p  core- 
exciton  binding  energy  and  linewidth  in  Si^Gei-j,  alloys.  The  observed  minimum  in  the  iinewidth 
near  x^0.15  can  be  explained  as  the  result  of  a  competition  between  intrinsic  broadening  due  to 
screening  and  extrinsic  alloy  broadening.  For  pure  Si.  the  binding  energy  is  estimated  to  be 
0.15  ±0.05  eV  and  the  width  is  shown  to  be  smaller  than  that  observed  at  x  =»  0.15. 


PACS  numbers:  71.5S.Fr.  71.35.  -H.  78.70.Dm 

Until  1984,  the  Si  2pcore  exciton  was  believed  to 
have  an  anomalously  large  binding  energy. Later, 
Newman  and  Dow*'  proposed  a  radically  dilTerent  pic¬ 
ture  in  which  the  Si  2p  core  exciton  is  in  fact  a  reso¬ 
nance  with  a  negative  binding  energy.  They  further 
predicted  that  the  exciton  binding  energy  remains  neg¬ 
ative  throughout  most  of  the  Si^Ge).,  alloy  composi¬ 
tion  range,  except  near  x  —  0.20  where  it  becomes 
positive.  In  a  recent  experiment,'^  Bunker  et  at.  found 
an  anomalous  sharpening  of  the  exciton  spectra  near 
x-O.IS;  the  data  were  interpreted  to  support  the 
Newman-Dow  point  of  view.  Yet  the  most  recent  ex¬ 
periment'^  still  suggests  a  positive  value  for  the  bind- 
^g  energy  £$  in  silicon. 

In  this  Letter,  we  present  a  calculation  that  offers  a 
plausible  resolution  to  the  above  problem.  In  our 
theory,  the  calculated  Si  2pcore-exciton  binding  ener¬ 
gy  E^ix)  and  the  Iinewidth  A(x)  in  Si^Gei-j,  alloys 
are  sensitively  dependent  on  the  parameter  £«.  A 
comparison  of  the  calculated  A(x)  with  the  experi¬ 
ment'^  suggests  a  positive  value  0.15  ±0.05  for  £§. 
The  anomalous  experimental  spectrum'^  near  x  ^  0. 1 5 
is  explained  as  a  result  of  a  competition  between  an  in¬ 
trinsic  broadening  A/  due  to  screening  and  an  extrinsic 
alloy  broadening  In  the  present  theory,  there  is  no 
need  to  suppose  that  the  exciton  suddenly  changes  its 
character  from  an  extended  effective-mass-like  state  to 
a  deep  localized  sute. 

We  need  to  calculate  £*  and  A-Ay-l-A^  as  a  func¬ 
tion  of  alloy  concentration  x.  The  calculations  are 
based  on  a  quantitative  coherent-potential-approx¬ 
imation  (CPA)  band  structure.  Details  of  the  CPA 
calculations  will  be  presented  elsewhere.  Below,  we 
briefly  discuss  a  Green's-function  method  for  calculat¬ 
ing  £»  and  A^. 

The  one-particle  effective  Green's  function  in  CPA 
takes  the  form 

C(£)-l£- W-l(£)l“'.  (I) 

^here  //  is  the  virtual-crystal  approximation  Hamil- 
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tonian  and  I(£)  is  the  self-energy.  The  site-diagonal 
Green’s  function  is  denoted  as 

£.(f)-(0alG(£)l6.>.  (2) 

where  <ba  is  a  localized  orbital  of  specified  symmetry. 
Here  we  only  need  to  consider  a  -  s  for  A\  symmetry. 
The  corresponding  function  in  pure  Si  is  denoted  as 
Ff(E).  Following  the  theoretical  treatment  of  deep 
substitutional-impurity  levels,^  the  core-exciton  level 
for  pure  Si  is  determined  by 

£«(£)-(  K-£»)-'.  (3) 

where  £/'  is  the  site  potential  seen  by  an  s  electron  in 
bulk  silicon,  and  F  is  a  central-potential  parameter. 
For  a  chosen  value  of  P*,  Eq.  (3)  can  be  solved  for  £. 
and  vice  versa.  Then  £°  “  £°  -*  £.  where  £°  is  the 
conduction-band  edge  in  pure  silicon.  Because  of  the 
uncertainties  in  the  value  of  experimental  £°  and 
theoretical  K  we  treat  £®  (or  P)  as  a  parameter.  The 
binding  energy  £»  in  a  Si^^Gei-^  alloy  can  be  calculat¬ 
ed  by  solving 

£,(£)-(P-£,-I,(£)l-'.  (4) 

where 

£,-x£/'-l-(l-x)£°*.  (5) 

Then  £*  is  given  by 

£»-£,-£.  (6) 

The  calculated  val '  s  of  the  conduction-band  edge 

and  the  exciton  level  measured  relative  to  the  top  of 

the  valence  band  are  plotted  in  Fig.  1.  The  band  gap 
increases  with  x  with  a  slope  discontinuity  at  =015. 
The  dashed  lines  a.  b,  and  f  represent  exciton  levels 
obtained  with  £°-0  1.  0.15,  and  0.30  eV.  respective¬ 
ly.  The  binding  energy  £*  is  also  an  increasing  func¬ 
tion  of  X,  with  a  slope  discontinuity  near  .v  =  0  15 
The  CPA  introduces  a  slight  bowing  m  £,  and  £» 
Stnnati"'  has  calculated  the  variation  of  A/  wnh  /;» 
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FIG.  1.  Variation  of  the  band  gap  (solid  line)  and  the  Si 
2p  core-exciton  level  with  x  in  Si,Gei_x  alloys.  The  energy 
is  measured  from  the  top  of  the  valence  band.  The  dashed 
curves  a,  b.  and  c  represent  exciton  levels  calculated  with 
£*  -  0. 1 ,  0. 1 5,  and  0.3  eV,  respectively. 


by  replacing  the  short-range  Coulomb  potential  with  a 
spherical  square  well  of  variable  depth  and  a  screened 
Coulomb  tail.  Strinati’s  results  can  be  used  to  estimate 
corresponding  to  the  calculated  £».  A/  decreases 
rapidly  with  £4,  then  saturates  for  larger  £4. 

The  contribution  to  the  natural  linewidth  from  the 
alloy  broadening  is  calculated  by  a  consideration  of  the 
electron  part  of  the  exciton  wave  function,  il/,.  The  <1/, 
is  expanded  in  a  linear  combination  of  the  spart  of  the 
conduction-band  wave  functions  <b^: 

«li,(k)-2;C„<6i(k).  (7) 

We  found  that  alloy  scattering  is  only  moderate  and  s 
scattering  is  dominant,  thus,  the  alloy  broadening 
( £)  is  well  approximated  by 

A,,(£)  =  jr(l-r)6/lm£,(£).  (8) 

where  5,  is  the  difference  between  £,^'  and  £,^* 
Hence,  the  alloy-broadening  contribution  to  A  is  relat¬ 
ed  to  the  alloy  broadening  of  the  band  states. 
A,(k.E): 

^  k 

--c22C,^A,(k.£) 

^  k  n 

=  /p,(£)A, (£)</£ 

=  X  (  1 -  v)ft;TrJp-(  £)</£.  i9) 

The  integral  in  Fq  (9)  is  evaluated  numerically 


FIG.  2.  Variation  of  A  (solid  lines)  and  A/  (dashed  lines) 
with  jrfor  three  f*  values 


The  calculated  A,  which  is  the  sum  of  A/  and  A<  .  is 
plotted  against  xin  Fig.  2  for  three  values  of  £*.  In  all 
three  panels,  the  dashed  curve  represents  A/  and  the 
solid  line  represents  A.  It  is  seen  from  Fig.  1  that  the 
exciton  level  follows  the  -V  edge  of  the  conduction 
band.  Hence  the  binding  energy  £4,  relative  to  the 
conduction  band  edge,  remains  almost  constant  (for  a 
given  £4)  until  the  minimum  switches  from  the  X 
edge  to  the  L  edge.  Because  of  the  change  in  the  slope 
of  £,.  £4  decreases  rapidly  when  L  becomes  the 
minimum.  Correspondingly,  A/  varies  slowly  until  the 
X  to  L  crossover  and  then  increases  rapidly  This 
feature  is  clearly  seen  in  Fig.  2. 

For  £4°  =  0  15,  the  A/  and  are  comparable  near 
V  =  0  50,  and  A/  dominates  for  all  small  v  and  large  v 
These  two  competing  mechanisms  give  a  relative 
minimum  near  x=OI5,  a  broader  maximum  near 
X  =  0  50,  and  a  smaller  minimum  lor  pure  silicon  .^s 
IS  decreased,  the  relative  minimum  is  shifted  10 
larger  x.  e  g  .  the  minimum  shifts  10  v  =  0  rd  for 
£4^  =  0  10  eV  For  £4’  =  0  15  cV  .  the  position  of  ihe 
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relulive  minimum  is  in  agreement  with  the  experi¬ 
ment.'^  (By  measuring  the  relative  width  at  .taiO.lS 
to  that  jf  -  1,  one  can  make  a  better  estimate  of  £°.) 
To  correlate  the  theory  with  experiment,  the  calculated 
1/A^  is  compared  with  the  measured'^ 
dE)nt,^  in  Fig  3.  where  Am  is  the  edge  step  and 
(dfi/dE)^^  are  the  maximum  values  of  the  derivative 
of  absorption  spectra  with  respect  to  photon  energy. 
Because  the  experimental  values  are  given  in  arbitrary 
units,  the  values  are  normalized  to  agree  at  t-0.5 
The  observed  anomalous  behavior  near  x  —  0.15  and 
the  qualitative  x  dependence  in  that  region  is  clearly 
replicated  by  the  theory.  However,  the  calculation 
predicts  a  larger  maximum  at  xr  —  1.  It  would  be  in¬ 
teresting  to  have  experiments  that  cover  the  entire 
range  of  xto  further  test  this  prediction. 

For  larger  values  of  E^,  the  calculated  £»  is  also 
large  and  hence  A/  decreases  slowly  with  x  Because 
the  broadening  is  determined  mainly  by  A^,  the 
linewidth  is  expected  to  be  small  for  jr-0  and  x-  1 
only;  this  occurs  for  £j  —  0.30  eV.  For  negative  values 
of  £°.  £4  remains  negative  for  all  values  of  x  Accord¬ 
ingly.  the  linewidth  is  broad  for  all  x.  and  there  would 
be  no  such  anomaly  as  in  Fig.  3. 

The  calculations  presented  in  this  Letter  are  slightly 
different  from  alchemy  approximations.*  We  treat  the 
central-cell  potential  as  a  parameter  and  narrow  its 
range  from  other  considerations.  We  examine  values 
of  F-£,*'  of  -8.49.  -7.09,  and  -6.56  eV.  corre¬ 
sponding  to  Ei  values  of  0.30.  0.15,  and  0.10  eV, 
respectively.  If  the  strict  alchemy  approximation  were 
taken,  the  value  of  F-  £/'  would  be  £/-  £,*'-■  -4.59 
eV  in  the  tight-binding  approximation,  and  a  negative 


KIG  3  The  tdituljicd  I/A’  vjIuc'  (solid  line)  compured 
wiih  the  experimenul  rcsulis  (nurks)  Irom  Rcl  12  The 
s.iluc  ol  e.ikul.iicd  AiSi>  is  (H  2'^  eV  The  expcnmeni  is 
iiirm;ili/ed  to  the  ihcnrs  .it  \  =  0  5 


£4  (  —  -0.10  eV)  would  be  obtained  When  long- 
range  interactions  are  included,  however,  the  above 
resonance  state  becomes  a  shallow  donor  level,  which 
is  the  experimental  situation  for  a  P  impurity  in  Si 
Our  results  suggest  that  P  for  core  excitons  in  Si  is 
deeper  than  those  implied  by  alchemy  approximations. 
However,  if  we  use  the  alchemy  approximation  as  a 
means  of  scaling,  the  value  of  y  for  Ge  3p  core  exci¬ 
tons  should  be  deeper  than  that  for  Si  2pcore  excitons 
Hence,  the  curve  corresponding  to  £4*- 0.30  in  Fig.  2 
should  be  a  reasonable  estimate  for  Ge  3pcore-exciton 
binding  energy  in  alloys  Therefore,  we  do  not  expect 
to  see  an  anomalous  behavior  of  A  in  alloys  for  this 
case. 

In  summary,  the  present  calculations  of  the  Si  2p 
core-exciton  binding  energy  and  linewidth  suggest  that 
the  exciton  level  is  about  0.15  ±0.05  eV  below  the 
conduction-band  edge  for  pure  Si.  It  follows  the  X 
edge  for  X  >  0.15  in  the  Si^^Gei.,  alloys,  and  £4  may 
eventually  reach  zero  in  the  dilute  limit  x— 0.  Our 
value  for  E°  represents  the  lower  end  of  the  previous 
measurements, but  is  in  good  agreement  with  a  re¬ 
cent  experimental'^  value  of  0.120  ±0.03  eV.  By  con¬ 
sidering  the  intrinsic  linewidth  and  the  alloy  broaden¬ 
ing.  we  can  explain  the  observed  relative  minimum  in 
the  linewidth  near  x  =  0.15.  without  requiring  a  sud¬ 
den  change  of  the  exciton  character.  On  the  basis  of 
this  calculation,  we  expect  the  corresponding  width  in 
pure  Si  to  be  even  smaller  than  that  observed  near 
x-0.15.  We  further  argue  that  the  binding  energy  of 
Ge  3pcore  excitons  should  be  larger  than  that  of  Si  2p 
core  excitons  and  there  should  be  no  anomaly  in  the 
Ge  3p  linewidth  in  alloys 
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Abstract 

Correlations  between  atoms  of  Cd  and  Zn  in  Cdi_,Zn,Te  and  Hg  and  Zn 
in  Hgi_,Zn,Te  are  calculated  within  an  extended  quasichemical  approximation. 
A  tetrahedron  of  alloy  atoms  is  used  for  the  microcluster  unit.  The  statistical 
mechanics  of  the  correlations  is  outlined  for  the  present  problem,  an  approximate 
method  of  solution  is  indicated.  Microcluster  energies  are  obtained  from  an  ap¬ 
proximation  of  a  Keating  model  due  to  Martin,  plus  a  correction  due  to  Harrison 
to  account  for  the  coupling  between  the  tetrahedral  bonds. 

It  is  found  that  the  correlations  are  significant  at  typical  growth  temperatures 
(300-800K).  Hgi_,  Zn,Te  is  predicted  to  undergo  spinoidal  decomposition  between 
300K  and  600K. 
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Introduction 


We  present  here  a  preliminary  investigation  of  the  statistical  mechanics  of  mul* 
tiatom  correlations  (microclusters)  in  pseudobinary  alloys  Ai_,B,C  of  tetrahedral 
semiconductors.  Thermodynamic  equilibrium  is  a  balance  between  entropic  forces, 
which  favors  random  site  occupation  of  the  pseudobinary  elements,  and  more  ener¬ 
getically  favorable  ordered  configurations.  Below  a  critical  temperature,  the  equi¬ 
librium  state  is  a  generally  an  ordered  state,  spinoidal  decomposition  for  positive 
values  of  the  mixing  enthalpy,  or  compound  formation  for  negative  values.  For 
growth  temperature  above  the  critical  temperature  the  alloy  may  still  exhibit  cor¬ 
relations,  which  vanish  in  the  infinite  temperature  limit. 

With  suitable  approximations  this  alloy  problem  is  isomorphic  to  the  well- 
known  Ising  model,  a  notoriously  difficult  problem  to  solve.  Since  a  rigorous 
treatment  of  correlations  between  units  is  intractable,  we  make  here  a  quasichemical 
approximation  and  treat  the  alloy  as  composed  of  microclusters  of  several  atoms 
embedded  in  an  effective  medium,  which  is  tantamount  to  neglecting  correlations 
among  microclusters.  The  problem  reduces  to  the  original  one  in  the  limit  of  large 
microcluster  size;  in  general  the  errors  introduced  by  neglecting  coupling  between 
microclusters  should  not  be  of  great  importance  well  above  the  phase  transition 
temperature.  This  model  is  clearly  inappropriate  at  or  below  a  critical  temperature, 
but  should  give  a  reasonable  indication  of  the  short  range  correlations  provided  the 
correlation  lengths  do  not  exceed  the  microcluster  size. 


Microcluster  Hamiltonian 


Tetrahedron  as  a  Microcluster  Unit 


If  the  microcluster  has  only  pair  interactions,  the  free  energy  must  be 


symmetric  in  x  and  1— z.  For  that  reason  we  will  work  with  a  larger  unit  and  choose 
the  tetrahedral  unit  shown  in  Fig.  1  for  the  microcluster.  It  is  the  natural  unit 
for  calculating  corrections  to  pair  interactions,  and  is  also  compact.  In  tetrahedral 
solids  the  strains  associated  with  an  impurity  tend  to  be  short  ranged;  thus  the 
tetrahedron  will  emphasize  the  most  important  differences  between  the  various 
configurations.  The  microcluster  is  taken  to  be  embedded  in  an  effective  (virtual 
crystal)  medium. 

A  microcluster  of  two  alloy  atoms  corresponds  to  the  Bethe-Pierls  ap* 
proximation  [1]  in  the  Ising  model  analogue.  The  quasichemical  approximation 
used  here  is  comparable  to  it,  though  somewhat  better  (as  larger  microcluster  sizes 
are  taken)  to  that  approximation. 

Force  Constant  Model 

Harrison  ( 2  ]  has  shown  that  the  cohesive  energy  in  a  tetrahedral  crystal 
derives  mainly  from  nearly  independent  bonds  that  only  weakly  couple  to  their 
environment,  and  that  the  variation  in  the  bond  energy  depends  mainly  on  the  bond 
length,  with  rather  weak  dependences  on  the  constituent  atoms  or  the  environment 
of  the  bond.  Metallization,  the  relatively  small  coupling  of  the  bond  to  neighboring 
bonds  Harrison  calculates  in  perturbation  theory.  This  implies  that  the  electronic 
structure  contributing  to  the  bonding  is  quite  short-ranged.  Long  range  interactions 
arising  from  the  electronic  structure  (apart  from  the  Madelung  terms)  can  be 
estimated  [3];  they  decay  exponentially  with  internuclear  distance  and  are  too 
small  to  be  of  significance  here. 

We  use  these  findings  as  a  justification  for  constructing  a  force  constant 
model  in  which  the  bonds  are  treated  nearly  independently  and  the  cohesive  energy 
in  the  alloy  differs  from  the  virtual  crystal  approximation  (VGA)  partly  through 
strains  of  bonds  stretched  or  bond  angles  bent  from  their  “natural"  positions. 
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Indeed,  such  a  model  can  be  derived  from  from  a  tight-binding  theory  [  4  ]  in  that 
theory  in  the  only  additional  corrections  are  the  changes  in  metallization  in  the 
microclusters  relative  to  those  in  the  VGA.  These  corrections,  called  previously  by  us 
[5]  the  “chemical”  energy,  are  also  included  here.  It  is  this  coupling  that  gives  rise 
to  departures  from  the  pair  interactions  and  asymmetries  in  the  dependence  of  the 
free  energy  on  composition.  In  the  heavier  elements  the  chemical  energy  can  be  quite 
significant,  and  can  approach  the  order  of  the  strain  energies.  Following  our  earlier 
work  we  will  use  a  Keating  model  (6),  neglecting  the  Madelung  contributions  to 
the  force  contents.  This  model  is  essentially  similar  to  that  of  Sokel  [4],  except 
that  we  fit  the  force  constants  to  experimental  data. 


Strain  Energies  in  an  Effective  Medium 

A  microcluster  differs  in  size  and  shape  from  the  effective  medium  sur¬ 
rounding  it,  and  so  displaces  of  all  of  the  surrounding  atoms  so  as  to  minimize  the 
total  strain  energy.  (The  chemical  energy  does  not  depend  strongly  on  the  strains 
and  can  be  ignored  when  calculating  the  strain  energy.)  In  the  approximations  made 
here  the  entire  effective  medium  is  characterized  in  terms  of  the  bonds  connecting 
the  microcluster  to  its  twelve  surrounding  second  neighbors.  Each  second  neighbor 
couples  to  the  adjacent  microcluster  atom  with  effective  force  constants  AT,  and 
Kt  for  strains  parallel  and  perpendicular  to  the  second-neighbor  bond,  respectively. 
Further  reduction  in  the  strain  energy  through  displacement  of  more  distant  neigh¬ 
bors  enters  here  as  a  modification  of  K,  and  K(.  So  far  this  is  merely  a  recasting  of 
the  Keating  model,  since  the  interaction  with  the  second  neighbors  can  obviously 
be  represented  in  terms  of  effective  force  constants;  the  approximations  lie  in  our 
choice  for  K,  and  Kt-  It  is  implicitly  assumed  that  the  24  force  constants  K  are 
independent  of  each  other.  Because  displacements  of  more  distant  neighbors  fall  off 
rapidly  away  from  an  impurity,  the  coupling  between  the  different  force  constants  is 
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indeed  quite  small,  and  errors  introduced  by  this  assumption  should  not  be  serious 
in  comparison  to  errors  inherent  in  the  effective  medium  approximation. 

To  guide  us  to  a  suitable  model  for  K,  and  Kt,  we  use  for  the  present 
purposes  the  Shih  approximation  [  7  ]  for  the  strain  energy  outside  the  microcluster 
itself.  Using  a  force  constant  model,  Shih  et  al  calculated  the  bond  length  of  a 
single  impurity  in  a  host  crystal,  holding  the  second  neighbors  fixed  and  neglecting 
the  smaller  angular  forces.  One  might  intuitively  expect  that  the  distortion  about 
an  impurity  is  quite  well  localised  in  tetrahedral  solids,  since  the  impurity-second 
neighbor  bond  angles  are  quite  large.  One  would  expect  this  to  be  far  less  the  case 
if  atoms  were  arranged  collinearly  as  they  are  in  the  alkali  halides,  for  example. 
In  tetrahedral  solids  a  displacement  along  the  impurity-nearest  neighbor  axis~biia 
only  a  small  component  along  the  neighbor-second  neighbor  axis,  and  the  distortion 
tends  to  damp  out  rapidly  as  the  angular  forces  are  weak.  The  simple  model  of 
Shih  el  al  predicts  that  the  nearest  neighbor  displaces  so  as  to  accommodate  3/4  of 
the  strain  in  the  impurity  bond.  This  is  in  excellent  agreement  with  experimental 
data  [8],  and  provides  poet  facto  justification  for  the  argument  that  distortions  are 
reasonably  well  localized.  Subsequent  analysis  by  ourselves  [  5  ]  has  shown  that  the 
weak  angular  forces  and  the  relaxation  of  more  distant  neighbors  make  corrections 
to  this  that  nearly  cancel.  The  Shih  approximation  abo  predicts  reasonable  mixing 
enthalpies.  Thus,  for  the  present  analysis  we  use  the  Shih  approximation  to  calculate 
the  effective  force  constants  Kr  and  Kt  ■ 

The  Shih  approximation  predicts 


where  or  is  a  Keating  force  constant  to  be  defined  below,  and 

a:,  =0 

because  it  neglects  angular  forces.  One  might  worry  that  the  latter  approximation  is 
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unrealistic,  but  a  detailed  analysis  [9]  shows  that  it  only  makes  a  small  contribution 
(~  3%)  to  the  total  strain  energy. 

Microcluster  Energies  in  an  Effective  Medium 

Following  our  earlier  work  [5],  strain  energies  are  calculated  from  the 
Keating  model 

^strain  =  -^  a  ^-r,  r.-  +  ^  ^  ■  r,  .  (1) 

above  signifies  the  square  of  the  difference  in  the  dot  product  relative  to  the 
unstrained  value.  The  a  terms  are  twobody  (radial)  forces  and  the  ,3  terms  are 
three-body  (angular)  forces,  and  are  typically  5-lOx  smaller  than  the  a  terms. 

The  total  strain  energy  is  thus  formulated  as  the  sum  of  strains  internal 
to  the  microcluster,  which  we  calculate  from  the  Keating  model,  the  three-body  3 
strain  energies  connecting  the  microcluster  to  the  adjacent  effective  medium  atoms, 
and  the  strain  energy  from  the  effective  medium, 

(2) 

»7 

where  the  index  »  signifies  each  of  the  three  external  bonds  adjacent  to  microcluster 
atom  j. 

The  strain  energy  must  be  minimized  with  respect  to  the  various  possible 
distortions  of  the  microcluster.  To  this  is  added  the  chemical  terms  [  5  ]  for  the  total 
energy.  The  two  can  be  reasonably  treated  separately  since  the  chemical  terms 
do  not  depend  strongly  on  bond  lengths,  and  thus  do  not  modify  the  distortions 
significantly. 

From  symmetry  considerations  the  and  B4  microclusters  have  only  the 
breathing  mode  shown  in  Fig.  2a,  in  which  every  atom  displaces  radially  by  an 
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equal  amount.  The  resulting  strain  energies  for  the  A4  and  B4  microclusters  can  be 
written 


•EB4“^o^trwn  “  where 


6A—Q-dA)R  and 

Sg  a=  (3  —  dg)/"! 


(4) 


are  the  fractional  differences  in  bond  lengths  of  the  pure  crystab  A  and  B  relative  to 
that  of  the  virtual  crystal.  The  A2B2  microcluster  has  five  independent  dbtortions 
rtj,  where  the  uj  represent  some  linear  combination  of  the  dbplacements  of  the 
microcluster  atoms.  With  judicious  selection  of  the  linear  combinations,  the  stress 
matrix  V,,.  Egtraio  can  be  made  block  diagonal  with  two  of  the  dbtortions 
corresponding  to  the  same  breathing  mode  as  in  A(,  and  the  remaining  three 
corresponding  to  an  internal  strain  shown  in  Fig.  26.  The  microcluster  strain  energy 
takes  the  form 


The  A3B  and  AB3  microclusters  have  four  independent  distortions,  but  again  a 
breathing  mode  can  be  block  diagonalized  with  the  strain  energy  taking  the  form 


with  a  similar  expression  for  £3.  C4  and  Co  depend  on  the  choice  of  model  taken 
for  the  Kr  and  Kt.  We  use  the  Shih  approximation  described  earlier  for  their 
calculation. 


The  total  energy  for  microcluster  j  b  then 


—  ^/.strain  4"  ^y.chem  i* ) 

The  chemical  energy  b  calculated  as  in  our  previous  work  [5].  Since  only  the 
differences  in  microcluster  energies  are  relevant,  it  b  convenient  to  define  linear 
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(8) 


combiaation  of  the  Eji 

By  definition  Ao  and  A4  are  0. 

Statistics  of  Correlations 

The  entropy  is  obtained  from  the  number  of  ways  to  configure  the  alloy  for  a 
given  energy.  Since  the  microclusters  are  assumed  to  be  independent,  the  energy 
of  a  configuration  is  specified  by  the  number  of  each  kind  of  microcluster.  For 
the  present  work  it  is  assumed  that  the  surrounding  medium  has  cubic  symmetry, 
so  that  the  microcluster  energy  depends  only  on  the  number  of  A  (or  B)  atoms  it 
contains.  We  distinguish  each  kind  of  microcluster  with  an  index  which  can  be 
taken  in  this  symmetric  case  as  the  number  of  A  atoms  in  the  microcluster,  so  that  j 
ranges  from  0  to  ns=4.  We  define  is  the  number  of  distinguishable  ways  of  rotating 
microcluster  j  on  the  same  site,  i.e.  =  1,4, 6, 4  and  1  for  ;  =  0, 1,2,3  and  4.  If 
an  external  shear  b  imposed  these  these  energies  are  not  degenerate,  in  which  case 
it  would  be  necessary  to  treat  each  orientation  separately. 

The  entropy,  the  logarithm  of  the  number  of  ways  to  configure  a  system  with 
jV/y  microclusters  of  type  j,  b  difficult  to  calculate  precisely.  We  use  here  the 
approximate  formula  [  10  ] 


where  and  Ng  are  the  number  of  A  and  B  atoms,  respectively,  with  the  number 
of  unit  celb  N  =  Ny^  +  Ng  and  the  number  of  clusters  M  =  Mj,  which  for 
thb  problem,  .V/  =  N.  The  products  over  j  range  from  0  to  n=4.  Finally 

x?  =  gyr'(l (10) 
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is  the  probability  of  finding  microcluster  ;  in  a  random  alloy.  The  expression 
in  brackets  in  Eqn.  0  is  the  probability  of  finding  Mo,  Mi,...M^  microclusters 
in  the  random  alloy  out  of  all  the  possible  configurations  of  microclusters.  The 
probabilities  i®  of  findic*  microcluster  j  are  assumed  to  be  independent. 


The  microcluster  energies  Ay,  Eqns..  8,  and  the  entropy  Eqn.  9  completely 
specify  the  equilibrium  statistical  mechanics  of  the  alloy.  With  Stirling’s  approxima¬ 
tion,  the  mixing  free  energy  is 

M  (l-i)(E4~^)  +  *(£b-£o) 

■+■  NkT[x]3ix  +  (1  —  x)ln(l  —  z)]  -I-  M^^XjAj  -I-  ArTzy  ln(zy/zy)  (11) 

i 

The  meaning  of  £4  and  £§  are  the  cluster  energies  in  a  pure  A  and  B  materials 
respectively:  (1  —  z)£4  -H  zEy  is  the  VCA  alloy  energy.  The  equilibrium  energy 
is  the  minimum  with  respect  to  the  microcluster  concentrations  Zy ,  leading  to  five 
simultaneous  equations  subject  to  the  constraints 

(1  -  z)n  =  53  j  *;• 

l  =  i;xy 

} 

Using  the  method  of  steepest  descents,  these  equations  may  be  collapsed  [9]  into 
a  single  equation  for  an  effective  activity  coefficient  tf.  Once  q  is  obtained,  the 
microcluster  concentrations  are  calculated  from 


where  Z  is  the  partition  function 


(13) 

(14) 

The  mixing  free  energy  oecomes 


SF=  N 


(1  -  x}(E4  -  £4)  +  z(£o  -  £^) 


NitT((l  —  n)(z  Inz  +  (1  —  z)ln(l  —  z))  +  nx  Inr/  —  InZ] 


(15) 
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Results 


Figs.  3  show  the  microcluster  energies  as  a  function  of  composition  for  the 
alloys  Cdi_,Zn,Te  and  Hgi_,Zn,Te.  The  asymmetries  in  those  curves  result  from 
the  chemical  terms,  and  the  asymmetry  at  the  extremes  of  alloy  composition  give  a 
rough  measure  of  their  importance  relative  to  the  strain  energy.  Near  the  extremes 
of  the  composition  range,  the  microcluster  energies  increase  roughly  in  proportion 
to  the  number  of  minority  species  present.  Near  the  5095  composition,  however,  the 
CdoZno  composition  has  very  little  excess  energy,  reflecting  the  fact  that  it  nearly 
fits  into  the  effective  medium,  while  the  Cd4  and  Zn^  microclusters  must  contract 
and  expand  respectively,  to  fit  into  the  room  available.  For  the  same  reason,  mutatis 
mutandis,  ZnCda  and  ZnaCd  show  the  smallest  energies  near  the  25^  and  75*^ 
compositions,  respectively. 

Figs  4  and  3  show  the  difference  in  expectation  values  of  microcluster  popula¬ 
tions  relative  to  their  random-alloy  values  at  T=300K  and  1500K,  respectively. 
These  fipires  show  strong  deviations  from  randomness  at  low  temperature,  and 
significant  ones  at  high  temperature.  Microclusters  having  least  strain  energy  show 
greatest  gains  in  population,  following  closely  the  rule  mentioned  above  for  the 
microcluster  energies. 

« 

The  free  energies  of  mixing  are  shown  in  Figs  6.  At  high  temperature,  the 
mixing  energy  is  nearly  symmetric  and  rather  strongly  negative,  reflecting  the 
dominance  of  entropy  terms.  At  300K  some  finer  structure  in  the  Cdi_,Zn,Te  is 
in  clear  evidence,  reflecting  a  tendency  toward  spinoidal  decomposition.  The  same 
tendencies  are  present,  and  more  pronounced  in  Hgi_,Zn,Te.  The  free  energy  is 
less  at  the  extremes  of  alloy  composition  in  Figs.  6  than  throughout  most  all  of 
the  composition  range,  thus  predicting  that  a  phase  transition  towards  spinoidal 
decomposition  occurs  in  Hgi_,Zn,Te  between  300K  and  600K. 
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Synopsis 

The  calculation  of  correlations  in  pseudobinary  alloys  is  outlined  within  a 
quasichemical  approximation,  using  tetrahedra  of  alloy  atoms  for  the  microcluster 
unit.  Using  a  Keating  model  with  corrections  for  coupling  between  the  bonds,  the 
energies  of  the  various  tetrahedral  microclusters  are  calculated  for  Cdi-,Zn,Te 
and  Hgi-aZn,Te.  The  correlations  are  found  to  be  significant  in  both  systems, 
particularly  in  Hgi-aZn,Te.  A  phase  transition  to  spinoidal  decomposition  is 
predicted  for  that  system  between  300K  and  600K. 
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Figure  Captions 

Fig.  1.  (a)  Tetrahedral  unit  of  four  alloy  atoms  (•)  surrounding  the  central 
atom  (o).  The  unit  b  embedded  in  an  effective  medium,  coupled  to  twelve  second 
neighbors  (o)  with  effective  force  constants  Kr  and  Kt  as  described  in  the  text. 
The  four  alloy  atoms  can  be  one  of  A4,  AB3,  A2B2,  A3B  or  B4.  (b)  the  external 
unit  connecting  to  the  cluster  atom  (•)  in  the  present  model.  Using  the  Shih 
approximation  the  bottom  three  atoms  are  assumed  not  to  move. 

Fig.  2.  Distortions  of  a  tetrahedron,  (a)  breathing  mode  distortion  of  A4 
and  also  in  the  other  microclusters  as  described  in  the  text,  (b)  distortion  in  A2B2 
corresponding  to  bond  length  differences  6y^  ~  Sg.  There  are  three  independent 
distortions:  a  radial  distortion  along  the  bond  axes,  a  distortion  transverse  to  the 
bond  axes,  and  the  displacement  of  the  central  C  atom,  (e)  distortion  in  A3B  or 
AB3  corresponding  to  bond  length  differences  S^—6g.  There  are  three  independent 
distortions  here  as  well:  a  uniform  translation  of  all  atoms  in  the  microcluster  (not 
shown),  the  displacement  of  the  central  atom,  and  the  displacement  of  the  three  A 
(or  B)  atoms  transverse  to  the  bond  axis. 

Fig.  3.  Microcluster  energies  as  a  function  of  composition  in  (a)  Cdi_,  Zn,Te 
and  Hgi_,Zn,Te. 

Fig.  4.  Microcluster  populations  relative  to  the  random  alloy  value  in 
Cdi_,Zn,Te  at  (a)  300K  and  (b)  1500K. 

Fig.  6,  Microcluster  populations  relative  to  the  random  alloy  value  in 
Hgi-,Zn,Te  at  (a)  300K  and  (b)  1500K. 

Fig.  6.  Mixing  free  energies  as  a  function  of  composition  in  (a)  Cdi_,  Zn,Te 
and  ('6y  Hgi_,Zn,Te. 
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